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PREFACE . 

My object in the following treatise is to apply to the Kinetic 
Theory of Gases a method of analysis different from that 
generally employed. It has been treated always on a certain 
fundamental assumption, namely, that the molecules of a gas 
are, as regards their relative motion, independent of one 
another. As a consequence, we may say as the expression, 
of that independence, the law of distribution of momenta 
assumes the exponential form e"'^^, and, so far as concerns 
translation velocities, 

m being mass, and w, v, w component velocities. From this 
independence and from this form of Q are deduced Boltzmann's 
theorems, namely the H theorem, and that of the equality 
of mean kinetic energy for each degree of freedom. 

I propose to give to Q the more general form of a quadratic 
function, namely 

Q = %m (w" + v^-\' w^) + SS6 (uu' + vv' + ww). 

Here 6 is a negative function of the distance r at the instant 
considered between the two molecules whose velocities are 
Uy u\ etc., which function is inappreciable except for very small 
values of r. I shall endeavour to prove in Chapters IV., V. 
that without the h coefficients the motion cannot be stationary. 
It has been proved abundantly that, assuming the independence, 
the motion is stationary with the usual form of Q. I question 
the axiom, not the demonstration. 

The consequence of attributing to Q the proposed new form 
is, that molecules near to each other have on average a motion 
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in the same direction. They tend to form streams. That 
result, if it can be established, is worth investigation. 

For ordinary gases under ordinary conditions the b co- 
efficients are probably very small, and their effect negligible 
in such investigations as those of Tait and Boltzmann con- 
cerning diffusion, viscosity, etc. But I think that the law 
^-hQ j^ j^g altered form will express the state of the system 
without restriction as to density, except as follows. A physical 
limit there must be, when the gas liquefies under pressure, 
if not before. For it will not be contended that the dis- 
tribution of momenta among the molecules of the liquefied 
gas is represented by the same exponential form as in the 
gaseous condition. An analytical limit there is, when Q in 
its altered form ceases to be necessarily positive, that is when 
the determinant of the coefficients ceases to be positive. It 
can be shown that this determinant does diminish as density 
increases, or temperature diminishes. But I have not cal- 
culated its value. It is therefore no more than a conjecture, 
though perhaps a plausible conjecture, that the vanishing of 
the determinant may coincide with the physical change in 
the substance. 

It appears to me that the law of equality of mean kinetic 
energy for each degree of freedom cannot be reconciled with 
my proposed form of Q; that in fact the law holds only for 
the limiting case of a very rare gas. 

It is no light thing to question a conclusion maintained 
by Boltzmann, if indeed he does maintain this conclusion for 
all substances, or for all gases irrespective of density. I can 
but state the objections to this theorem, and to a certain aspect 
of the H theorem, as they appear to me. The reader will 
judge what weight is to be attributed to them. 

S. H. BURBURY. 

10 Ma^, 1899. 



TABLE OF CONTENTS. 



Chapter I. Arts. 1—4. Statement of the theory. 5—10. The defini- 
tion of homogeneity, density, and stream motion. 11 — 13. Of the 
pressure of a gas. 14. Of the temperature. 15. Intermolecular 
forces defined. 16 — 20. Two fundamental assumptions A and B 
discussed pp. 1 — 13 

Chapter II. Arts. 21—26. Clausius' theorem. 27 — 28. Value of 
22i?r for elastic spheres. 29—30. A vertical column of. elastic 
sphere molecules pp. 14 — 23 

Chapter III. 31 — 37. Law of stationary motion for elastic spheres. 
38. The H theorem. 39—42. Discussion of the ff theorem. 
43 — 45. Application of Boltzmann's general method, making the 
density proportional to e~^^^, where x I'elates to external forces. 
46. Proposed inclusion of yj/ the potential of intermolecular forces. 
47 — 52. Boltzmann and Watson's generalization. Note on the 
ff theorem . . .pp. 24 — 45 

Chapter IV. Arts. 53 — 54. An extension of Clausius' theorem. 
55. Case of no intermolecular forces. 56 — 57. Correlation of 
velocities proved when intermolecular forces exist. 58. Possible 
case of attractive forces pp. 46 — 54 

Chapter V. Arts. 60—62. Definitions of /, ^, etc. and M. 63—68. 

Treatment of -r- in the case of infinitely small non-colliding mole- 
cules. 69 — 78. Treatment of -^ for elastic spheres having finite 
diameters pp. 55—78 

Chapter VI. Art. 79. Proposed form of solution for elastic spheres. 
80. Completion of the function M. 81— 82. Cases in illustration. 
83 — 84. The form of the coefficient b not determinate, but 2b is 
determinate. 85. A suggested form of (b) ; relation of m'2uu' + &c. 
to the Virial of intermolecular forces .... pp. 79 — 88 



Vlll TABLE OF CONTENTS. 

Chapter VII. Arts. 86 — 92. Treatment of molecules under finite inter- 
molecular forces, the encounters being supposed binary. 93. A sug- 
gested extension to non- binary encounters . . . pp. 89 — 97 

Chapter VIII. Arts. 94—96. The general law of distribution of 
momenta. 97. The fundamental property of the distribution. 

98 — 100. Boltzmann's law of equal kinetic energy. 101 — 102. Ap- 
plication of Boltzmann's general method giving for density ^""'^(x+w. 
103. The vertical column again considered. 104. Boltzmann's 

minimum function. 105 — 110. General character of the motion 

pp. 98—112 

Chapter IX. Arts. Ill — 112. Determination of mean free path. 
113 — 118. The fundamental theorem for Diffusion, &c. pp. 113—121 

Chapter X. On the relation between temperature and kinetic energy. 
Arts. 119—123. Natanson's theorem. 124. Method em- 
ployed by Bryan and Boltzmann. 125 — 126. By Professor 
J. J. Thomson. 127 — 131. The Second Law of Thermodynamics. 
132-133. Original proof. 134. Modified proof. 135—142. Other 
proofs of the same law pp. 122 — 142 

Appendix, (a) — (p) Certain propositions relating to determinants. 
(q) The condition that the quadratic function Q shall be necessarily 
positive. (r) On the value (referring to Chapters V., VI.) of 

l///"^S'^-^«'^^^ pp. 143-157 



CHAPTER I. 



OUTLINE OF THE THEORY. 



1. A GAS according to the Kinetic Theory consists of a 
great number of molecules in rapid motion. And the object of 
the theory is to explain on this hypothesis certain of the 
physical properties of gases. 

Any quantity of gas which can be isolated for the purpose 
of experiments is to be regarded as containing a number 
of molecules practically infinite. It is not possible to control 
or to observe the motions of individual molecules. But the 
theory assumes that such motion is subject to the usual 
dynamical laws. Also that if the gas, as an aggregate of 
molecules, be at rest, no dissipation of energy takes place. 

2. A molecule may consist of one or more than one atom 
according to the chemical constitution of the substance to 
which it belongs. It may be that hereafter we shall be able to 
explain on djniamical principles the chemical relations of atoms 
as constituent parts of a molecule, and of molecules inter se. 
And some progress has been made in this direction. At 
present the theory is concerned not with the chemical pro- 
perties, but with those properties of gases which may change 
without any change taking place in the chemical composition 
of the gas: for instance, density, pressure, and temperature. 
And as depending on these latter, it is concerned with the 
phenomena of viscosity, diffusion, and conduction of heat or 
electricity. 
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3. A molecule moving as such, one and indivisible, has 
three degrees of freedom in respect of its motion of translation 
in space. And it may conceivably have three other degrees of 
freedom, namely, motion of rotation about its principal axes. 
In respect of its internal constitution it may have many degrees 
of freedom, but we know nothing of the internal constitution 
of molecules except what chemists tell us. 

4. Any finite quantity of gas may, if protected from dis- 
turbance from without, retain for an indefinite time un- 
impaired its physical properties, for instance, its pressure and 
temperature, and also its chemical constitution. If therefore 
these properties are not changing, and are to be explained 
on dynamical principles as resulting from the motion of the 
molecules, such motion must be stationary motion. 



OF HOMOGENEITY. 

5. The gas, or system of molecules, may be in a certain 
sense, i.e. if we regard all the molecules at any instant within 
a finite space, homogeneous. For instance, a cubic centi- 
metre of air has sensibly the same properties in all respects 
whether taken from one part of a horizontal table or from 
another part. It is true that it will generally not have in all 
respects the same properties if taken from near the floor, as if 
taken from near the ceiling. It will, namely, be heavier, i.e. 
will contain more molecules, in the former than in the latter 
case. That is because the air in the room is in a field of 
uniform force, gravity. But this force will not necessarily 
affect any other property of the gas, at all events under 
ordinary conditions to which it is subject. 

If we assume the infinite divisibility of matter, or fchat any 
space, however infinitely small, contains an infinite number of 
molecules, the condition of homogeneity may apply to the 
molecules contained in any such infinitely small space. We 
may in that case say without ambiguity that the gas at a 
point P has the same properties as at another point Q. 
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6. Modern physicists however teach us that the molecules 
have finite dimensions, or that there is a finite distance, how- 
ever small, within which two molecules cannot approach each 
other. That being the case, the homogeneity can be asserted 
only in the following form. If we consider a quantity of gas 
within a sphere of radius r described about P as centre, 
and a quantity within an equal sphere described about Q as 
centre, either sphere containing a number of molecules practi- 
cally infinite, then r being within certain limits, the first 
quantity has the same physical properties as the second. 



OF THE DENSITY OF A GAS. 

7. The same doctrine of the finite size of molecules causes 
considerable difficulty in defining certain other properties of 
the system, for instance the density. 

The density of the gas is usually defined to be the sum of 
the masses of all the molecules in unit of volume. But if the 
molecules have finite dimensions, that may make the density 
at any point relative to that at other points depend on the unit 
arbitrarily chosen. It is not therefore a complete definition. 
We can however say, and it is generally sufficient for our 
purpose, 

(a) If there be N molecules within a definite volume S, 
the average density at any point within 8 is N/S. 

(b) Again, we can define the density at a point which 
answers some general description, or as we may call it the relative 
density. For instance, suppose / to be a continuous function 
of the position of a point, then we may define the density 
with reference to /, as for instance at points where /= 0, 
provided that there exists an infinite number of such points. 
For consider any n of them, and about each a small sphere 
of volume g). Let nco be the unit of volume. Let N be 
the number of molecules supposed of equal mass which, or the 
centres of inertia of which, are within some one or other of 
these small spheres. Then the density where y= is If/na) in 
the limit, as n is increased indefinitely, noy remaining constant. 

1-^2 
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To define density at a point, as a function of the coordinates 
of that point, is with finite molecules a more difficult matter. 

Let r be a line of any arbitrary fixed length. As a first 
step to our definition, let us say the density at a point P, or 
X, y, z is the number of molecules, supposed of equal mass, con- 
tained within a sphere of radius r described about P as centre, 
divided by the volume of that sphere. If n be that number, p 
the density, then our definition is 

But that makes p generally a function of r as well as of 
X, y, z. How are we to get rid of r ? Boltzmann and other 
writers after or before him would complete the definition thus, 

p at P = ~. — - when r becomes infinitely small, any space how- 

ever infinitely small being assumed, for this purpose only, to 
contain an infinitely great number of molecules. This defi- 
nition does not materially differ from that which we should 
give if our gas, instead of consisting of discrete molecules, were 
a continuous fluid, any portion of which may be conceived as 
infinitely subdivided. 

But it is inconsistent with the teaching of these same 
writers, who maintain that molecules have finite dimensions, 
that in fact we can ascertain limits between which those dimen- 
sions must lie. It may be admitted that no ill consequences 
have yet arisen from the use of this inconsistent definition, 
and if it be either necessary or useful, the inconsistency need 
not be regarded as a serious objection. Inconsistency may 
indeed be an advantage in mathematics as it often is in politics. 
It may be admitted also that this definition occurs to us by 
instinct when we are asked to define density at a point. That 
may be because we have passed through a former state of 
existence in which matter was supposed to be continuous, and 
therefore the definition strictly accurate. And our instinctive 
resort to it is due to unconscious memory of that former state. 

8. We may avoid the difficulty by making a further 
assumption with regard to the properties of our system. Let 
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Vq be the radius of a sphere described about P as centre which 
contains a very great number of molecules. Let Vi be greater 
than Vq in a finite ratio. Let as before n be the number of 
molecules within a sphere of radius r. Now let us assume, as 
a property of our system, that for all values of r between Tq and 
Viy the ratio n:r^ is constant. Then we may define /o at P as 

r— ^ using this constant value of n/r^. 

This method implies either that the state of our gas does 
not change from point to point in space, or that it changes 
so gradually as to be sensibly constant throughout the space 
between the spheres of radii Tq and Vi, 

9. Assuming that the state of our gas changes so rapidly 
in space that we cannot make the above assumption, and 
assuming that we scruple to base our definition on an assertion 
which we elsewhere contradict, it seems to me that the ex- 
pression "density at a point" or "number of molecules per 
unit of volume at a point " has of itself no meaning whatever 
when the molecules are of finite dimensions, but a meaning may 
be given it in an arbitrary way. One way would be to take an 
arbitrarily chosen volume for the unit. I would suggest the 
following as the best form of definition. Let /be a continuous 
positive function of ?^ which is equal to unity for all values of r 

less than a certain distance a, and for which ~ is always 

negative when r > a, and such that I rYdr is a convergent 

series. Then define the density at P as %mf, r being measured 
from P, m being the mass of a molecule and the summation 
including all molecules in space. For with molecules of finite 
dimensions, density at P must depend upon the distances of 
many molecules from P, but cannot be affected by those whose 
distances from P are very great. 

THE STREAM VELOCITY. 

10. As with the density, so precisely the same difiiculty 
presents itself if we try to define the stream velocity, or in 
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Professor Tait's language, mass motion, of the gas in given 
direction at any point P. M. Ladislas Natanson in a paper, 
to which I shall refer later, makes use of this function 
without giving any definition of it except by calling it 
" vitesse moyenne ou apparente." He is however, as I under- 
stand, willing to assume that any infinitely small space contains 
an infinite number of molecules. As the same difficulty arises, 
so the same solution mutatis mutandis may be resorted to, as in 
the case of density. 



OB^ THE PRESSURK OF A GAS. 

11. The motion which we attribute to our molecules cannot, 
assuming the gas to be at rest relatively to the surrounding 
space, be related to any particular direction in space. That is, 
at or near auy point P there are on the average of time as 
many molecules moving with given velocity in any one direction 
as in any other. For if otherwise the pressure of the gas would, 
as shown by the next paragraph, be unequal in different 
directions, which is known not to be the case if the gas be at 
rest. 

12. Assuming the property mentioned in the last paragraph, 
and assuming the distribution of momentum among the mole- 
cules to be known, we can calculate the quantity of momentum 
transferred across an imaginary plane within the gas per unit 
of area and time. Let m be the mass of a molecule, p the 
number of molecules per unit of volume, to be defined as above 
explained, and pf(V) dV the number per unit of volume whose 
velocity of translation lies between V and V + dV. Then by 
our assumption the number of these last whose direction of 
motion makes with the normal to the plane angles between 
6 and + dd is 

y,f(V)dVsmedd. 

Their velocity normal to the plane is Fcos 0. Therefore the 
quantity of momentum normal to the plane transferred across 
the plane in one direction per unit of area and time by 
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molecules having velocity F... V+dVis 



^mpV^ I sin cos^ 0d0 = ^7np V\ 
Jo 



And the quantity of momentum so transferred by all the 
molecules is 



Jo 



= ^mpV^ suppose, 
V^ being the mean value of V'^ for all the molecules. 

13. If further the molecules exert on one another no forces, 
twice the quantity of momentum so transferred through unit of 
area of the plane is the normal pressure of the gas per unit of 
area of the plane. See Watson s Kinetic Theory of Oases, 
second edition, pp. 56 — 58. 

For we may suppose the gas confined in a vertical cylinder 
under a heavy moveable piston from which molecules striking 
rebound as elastic bodies. The momentum normal to the base 
of the piston, i.e. vertically upwards, which if it were an 
imaginary plane would be transferred through it, is in fact, 
the mass of the piston being very great, reversed, so that twice 
that quantity is transferred to the piston, and tends to make 
the piston rise in the cylinder. That is, the piston receives 
per unit of area and time from the impacts of the molecules 
momentum equal to ^mp V^ vertically upwards. It also acquires 
by the action of gravity per unit of area and time momentum 
Mg vertically downwards, where g is the force of gravity, and 
M bears to the whole mass of the piston the same ratio which 
unit of area bears to the base of the piston. For equilibrium 

Mg = ^mp V^ = py the pressure per unit of area. 

If the impacts of the molecules be sufficiently numerous 
per unit of time the piston will remain sensibly at rest. 
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OF THE TEMPERATURE OF A GAS. 

14. We have then in the absence of intermolecular forces 

p varies as p F^. 

If the molecules do exert mutual forces on one another, the 
above value of the pressure is not accurate. But there are 
sufficient reasons for believing that for ordinary gases under 
ordinary conditions the pressure is not materially affected by 
the intermolecular forces. For such gases therefore, accurately 
or approximately, 

p varies as p F^ 

But for such gases, accurately or approximately, the laws of 
Boyle and Charles hold good, whence p varies as p0y 6 denoting 
the absolute temperature. It is thence inferred that, accurately 
or approximately, 6 varies as F^ And that statement we 
must accept provisionally. Temperature measured from the 
absolute zero is represented by the mean square of the velocity 
of translation. On this theory temperature is a quality which 
cannot be attributed to a molecule, but only to an aggregate of 
molecules. 

OF INTERMOLECULAR FORCES. 

15. It is however certain that for sufficiently dense gases 
intermolecular forces do exist, and we must now make a further 
assumption regarding them. We shall assume, namely, that the 
intermolecular force R between two molecules 

(1) Acts in the line joining their centres of inertia, or 
within very small limits of error may be regarded as so acting. 

(2) Is a continuous function of the distance r between 
those centres of inertia, which becomes evanescent for all values 
of r exceeding a very small finite distance ro, which distance 
again is negligible compared with the dimensions of any mass 
of gas upon which we can make experiments. 

(3) The force must for sufficiently small values of r become 
repulsive, and should become infinite as r is indefinitely dimin- 
ished — but save for these restrictions it may be that the force 
IS attractive for some values and repulsive for other values of r. 
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R shall be taken as positive when repulsive. 

If c be the smallest radius beyond which the force exerted 
by a molecule at the centre on another molecule is negligible, 
c is defined to be the radius of action of the central molecule, 
and a sphere of radius c described about the molecule is its 
sphere of action, or its effective volume. 

THE FUNDAMENTAL ASSUMPTION. 

16. As above stated, the number of molecules of our gas 
per unit of volume may be different at different points of the 
system considered. Representing by p that number, we may 
have to treat p as a function of x, y, z. At present it is con- 
venient to treat only the case in which p is constant throughout 
the system, and the distribution of velocities among the 
molecules is the same throughout the system. That being the 
case we may define as follows. 

At any instant let pf^ (w) du be the number per unit of 
volume of molecules whose component velocities in the direc- 
tion taken for x lie between u and u -h du. Similarly let 
pfy {v) dv be the number per unit of volume whose component 
velocities in direction y lie between v and v + dv, and pf^ (w) dw 
the number whose component velocities in direction z lie 
between w and w + dw. 

We might say, and it would be precisely equivalent to what 
we have already said, that f^. (a) du is the chance that a given 
molecule shall have velocity in x between u and u + du, and 
fy (v) dVyfz (w) dw are the corresponding chances for y and z, 

17. These definitions, however apparently unambiguous, 
may be based upon either of two diametrically opposite 
assumptions. 

We may make namely, assumption A, 

A. The chance of any molecule having velocity in x between 
u and u + du is independent, not only of its position in space, 
but also of the velocities v, w, which it has in directions y and z, 
and further except in the case mentioned below, it is inde- 
pendent of the positions and velocities at the instant of all the 
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other molecules of the system. The excepted case is when the 
two molecules are so placed that they are, or very recently have 
been, within one another's sphere of action. The force of this 
exception, and the necessity for it, will appear in the considera- 
tion of the H theorem, Chapter ill. 

On that assumption it follows that the chance of any given 
molecule having at any given instant velocities in Xy y, z 
respectively between xi and ii -h dii, v and v -\-dv,w and w + dw, 

is /^ {u) fy (v)f, (w) dudv dw. 

Further, if we consider any pair of molecules, and denote 
by the suffix 1 quantities relating to the first, by the suffix 2 
quantities relating to the second, then the chance that they 
shall at any instant respectively have velocities 

in X between u^ and u^ + du^, 

U2 and U2 + du2 ; 

in y between Vi and v^+dvi, 

V2 and V2 + dv2 ; 

and in z between Wi dbXid w^ + dwi, 

W2 and W2 + dw2, 

is /i^ (ui)fiy {vi)fiz {Wi)f2x {^2} . . . du^dvidwAh >>.dw2. 

This or its equivalent is the assumption on which the kinetic 
theory of gases has hitherto generally been treated*. I shall 
refer to the state of things assumed as condition A. Without 
doubt it may be assumed legitimately with respect to the 
limiting case of an infinitely rare gas, that is one in which 

* Dr Watson (Kinetic Theory of Gases, 2nd Edition, p. 3, equation B) 
assumes the independence of the chances for two spheres approaching collision. 
And that is sufficient for his purpose. Professor Tait ("Foundations of the 
Kinetic Theory of Gases," Transactions R. S. Edinburgh, 1888, Art. 21) makes 
an equivalent assumption. Dr Boltzmann [Vorlesungen iiher Gas Theone, 
Part I, p. 23) assumes that the motion is, and ever continues to be, 
'* molecular- ungeordnet. " I shall consider later (Chapter in.) what may be the 
effect of that assumption. In the meantime the weapon that Dr Boltzmann 
draws from that armoury is precisely the same as the one used by Watson and 
Tait, namely, that for all pairs of spheres or molecules approaching collision the 
chances are independent, as stated in my condition A. These writers do not 
define the exception, the occasion for definition not having arisen, but they do 
not assert the independence in the excepted case. 
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the dimensions of any molecule, or if it be a centre of force, 
its radius of action, are negligible compared with the mean 
distance between any molecule and its nearest neighbour. 
It is probable also that all the known gases under ordinary 
conditions of pressure and temperature, approximate closely to 
the limiting condition of infinite rarity. And therefore the 
calculations made on the assumption of condition A or its 
equivalent by Tait, Boltzmann, and others, of the rate of 
diffusion, conduction of heat or electricity, &c., for such gases 
remain unaffected. 

18. Or we may make assumption B. /^.. 

B. The chance of a given molecule having at any instant 
assigned velocities is n ot independ ent of the positions and 
velocities of all the other molecules at the instant. On this 
assumption B, instead of deducing the chance of the members 
of a group of n molecules having respectively at any instant the 
velocities 

u^,.,u^-\- du^ &c. 

from the assumed chances for individual molecules, we must 
reverse the process. Let the chance that at the given instant 
the coordinates of the n molecules shall lie respectively between 

Xi and Xi + dxi, 

Vi and yi + dy^, 

-01 and Zi + dzi] 
for the first molecule, 

X2 and X2 + dx2, 
&c. 
for the second, and so on, 

and their component velocities between 

tti and Ui + duiy 

Vi and Vi + dvi, 

Wi and Wi-\-dwi\ 



for the first molecule. 



U2 and U2 + du2 

&c., 
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for the second and so on, be denoted by 

F(xiyi . . . ZnUiVi . . . Wn) dx^dyi . . . dwn- 

Then the chance that the x velocity of the first molecule 
shall lie between 

Ui and i^i + dwi, 

whatever be the positions and velocities of the other n — 1 
molecules, is 



K 



(dxidx2 . . . dzndv^dwidu^_ . . . dw^^ F{xiyi ...Ui ... w^. 



19. By way of illustration let us suppose, under condition A, 
f{u) = Ge-^-\ 
where C is the usual constant, chosen so that 



I'OO 
J —CO 



6-'"''dM = l. 

J —CO 

Then on assumption A the chance that the members of the 
group of n molecules supposed to be of equal masses, shall 
respectively have their coordinates and velocities between the 
above limits is of the form 

Ce~^Mxi . . . dzndui . . . dwn , 
where Q = X(u^ + v^-{- w^), and h is constant ; 

and in dealing with assumption B, let us assume, the molecules 
being all of the same mass and the same structure, 

F(x,y,...u,... Wn) = C€-^^, 
where Q = u^^ + v^^ + w^^ + u^^ -\- ... 4- w„^; 

+ 6i2 {1I1U2 + V1V2 + W1W2) 

Here the b coefficients shall be functions of the distance r, at 
the instant between the pair of molecules to which the suffixes 
relate, which functions become evanescent when r exceeds a 
certain very small distance. In that case omitting the constant 



/(%)=/// 



+00 

€~^^ dxi . . . dzndvidwidu2 . . . dWn 



= 6 2^11 ' ' , as shown in the Appendix (g). 



2 


612 


bn .. 


tia 


2 


t^, .. 


6is 


b^ 


2 .. 
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Here J9 is the determinant of 71^ constituents 



D = 



and Dn is its first coaxial minor. 

The form oi f{u) under assumption B is the same, except 
for the factor of h^ as under assumption A, but the physical 
state of the system may be very different. For instance, under 
assumption B, two molecules at a distance r apart, so small 
that the corresponding h is not negligible, are, if h be negative, 
on average moving in the same direction. According to Boltz- 
mann (Vorlesungen, p. 21) the motion in this case is not 
" molecular ungeordnet.'' 

20. Of these two assumptions A and B, B includes A as a 
particular case. It is therefore better, if and in so far as it 
may turn out to be possible, to establish a proposition on 
assumption B, rather than on assumption A. My object in 
the present treatise is to show that it is generally possible. It 
will appear also that A requires us to assume our gas to be 
infinitely rare, whereas B imposes upon us no stringent con- 
dition with regard to the density. The only condition in fact 
which appears to be required at present is this : the coefficients 
b must be such that for all possible values of the us, vs, and 
w's, Q shall be positive. Mathematically that is expressed by 
saying that the determinant D and all its coaxial minors must 
be positive. 

D will appear to be generally a function of density and 
temperature, and if ever D = the mathematical treatment 
changes, whether or not the physical system changes its state, 
and if so whatever the nature of that change may be. 

It will appear also that if u, v, w be the component velocities 
of a molecule of mass m, u\ v\ w' those of any other molecule, 

\m {uu + vv' -h ww') 
is related to the ViHal defined in Chapter ii. of the intermolecular 
forces acting on m. See Arts. 58, 85. 



CHAPTER IL 



CLAUSIUS THEOREM. 



21. It will be necessary to consider at length later the law 
of distribution of velocities or momenta among the molecules 
for any given value of the mean kinetic energy — that is the 
form of the function which we denoted above, Art. 12, by /(F). 
The following proposition, due originally to Clausius, is inde- 
pendent of that law. 

If m be the mass of any one of a system of molecules in 
stationary motion, ^, y, z the coordinates of its centre of 
inertia at any instant, we have, S denoting summation for all 
the molecules of the system, and 2m^, Smi/, ^mz being zero 






dt 
dt 



or 2m f ^ j + 2m^^ = 0, 



with similar equations for y and z. 
It follows that 



2' 



m 



{(sy-(S)"-(S')=-^(^-''^--* 



X, F, Z being the component forces, whether external or 
intermolecular, acting on a molecule. 
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22. We will now however reserve X, F, Z for the com- 
ponents of external force only, and let R be the intermolecular 
force between two molecules acting in the line joining their 
centres, and let R^yRy, Rz be its components. So that Clausius' 
equation becomes 

-t2{R^-^Ryy^R,z) (1). 

Half the expression on the right-hand side of the last 
equation is called the Virial. And separately 

-\{Xx-¥Yy + Zz) 

is the Virial of the external forces, and — \1X{Rafc + Ryy + R^z) 
is the Virial of the intermolecular forces. 

23. Let one molecule be at P or x, y, z, another at P' or 
X, y\ z'. Let r^ = (^ - xj + {y- yj 4- (^ - zj. Then for the 
force on P due to P\ R being taken as positive when repulsive, 

Rx = Ri Ry = Ri Rz == — "~ R' 

Similarly for the force on P' due to P, 

Rx = Ri Ry = Ri Rz = Ri 

and therefore for two molecules at P and P' 

And so the term in Clausius* equation due to the inter- 
molecular forces is 

^^{R^ + Ryy-\-Rzz)=i;iRr (2), 

each pair of molecules being counted once. 

24. We have next to consider the external forces, whose 
components are X, F, Z. I shall follow Van der Waals' 
method. 

Let the whole system be enclosed by a surface /S, through 
which the molecules cannot pass. Let p be the normal pressure 
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per unit of area on S, directed inwards. At any point P on S 
let X, fjL, V be the direction cosines of the normal s to 5^ at P, 
I, m, n those of a the line joining P with the origin. Then for 
a molecule at P 

— {Xx + Fy + Zz) —p (l\a + m/xa + nva) 

=^ pa cos (s, a) 

da 

if ds be an element of the normal, measured outwards, 

, da" 
= ^P-d^' 

Let us now assume p to be constant over the whole surface 
,Sf. Then 



■// 



(Xx + Yy + Zz) dS 



, dS 
as 



= iP nj ^''cL^dxdydjz 

throughout the space enclosed by S, by Green's Theorem, 

= Spa), 
if ft) be the volume enclosed by S. 
We see then that the Virial 

orby(l) ^p(o = I.^m{u^ + v^ + w^) + iX%Rr (3); 

or, if T be the mean kinetic energy of the system of molecules, 

CO 

T 

and — represents the mean kinetic energy of a molecule. 

If we take the volume m for a generalised coordinate, 
supposing it to change in magnitude but not in shape, the 
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other coordinates of the system may be a?a) , yceTy &c., where 

T dT 
X, y, &c. are numerical. Then I - = 3~ > ^^d if % be the 

potential of the intermolecular forces ^XXEr/co = — -^ , and our 

equation assumes the Lagrangian form p — -i -r^ , p being 

the applied force which maintains a> constant. 

25. If the surface S be an imaginary closed surface, the 
system of molecules in the external space being homogeneous 
with the system within S, the pressure per unit of area of S 
is given by the above equation, 

The whole Virial or ^pco — ^22iJr is thus equal to the 
whole mean kinetic energy of the molecules within S. It 
does not follow that the two parts of the Virial, namely ^pa> 
and — Xl^^RTy are separately proportional to X^m (u^ + t;^ + w^). 

26. If S be the external elastic boundary of the system 
through which molecules cannot pass, we might calculate p by 
the same method as we calculated it in the case of there being 
no forces, namely p = ^mpV'^, where p is the mean density, and 
\mV'^ the mean kinetic energy, of molecules near S, This 
result is not inconsistent with Art. 24, because the mean 
potential energy near S is, in the case now supposed, different 
from the mean potential energy throughout the enclosed space, 
V^ and soT'2 ^ 72 ^^d p' + p. 

The vahie of XZRr for Elastic Spheres. 

27. We have assumed iZ to be a finite force. In the 
limiting ease of so-called elastic bodies, the force which acts in 
collision between two such bodies is, according to the usual 
convention, an infinite force which reverses the velocities of the 
colliding bodies normal to their tangent plane at collision in an 
infinitely short time. As a consequence of the time being 
infinitely short, no change takes place in the coordinates during 

B. 2 
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collision, and no body can ever be in collision with more than 
one other body at the same instant. 

It is of course admitted that no such body exists in nature. 
But the behaviour of gas molecules in regard to their mutual 
action is probably, at least in the case of rare gases, approxi- 
mately the same as it would be if they were elastic spheres of 
this form. And the assumption made in the case of the 
spheres that the kinetic energy of translation of two spheres in 
collision is in the aggregate unchanged by the collision, must 
if molecules are stable be true on average for any form of 
molecule. 

28. It is then possible to calculate X^Rr for a system of 
elastic spheres, regarding the force as finite, and sufficiently 
great to produce the reversal of the normal velocities in a 
finite time r, and treating the coordinates as constant, to proceed 
to the limit when the force becomes infinite, and r zero. 

Let the colliding spheres be each of mass m, and of 
diameter c. Let p be the number of spheres per unit of 
volume. 

Let T be now the mean kinetic energy of a sphere. And 
let w, V, w, u\ v\ w\ be the component velocities of any two 
spheres, V their relative velocity. Then 

V^ = (u- uj + (y - vy -\-(w- wj, 
and mF^ = 2m {u^ + ^^ + vJ^) on average 

Let 6 be the angle between the relative velocity and the 
line of centres at collision. Then for each sphere the normal 
component of relative velocity is \V qos $, and it is reversed in 
time T. The assumed force is then R = mV cos 0/t, and it acts 
during the time r. The number of collisions per unit of 
volume and time is, given F, irc^Vp, And the number for 
which d lies between and 6 -^-dd is 

27rc^Vp sin cos dd0. 
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Also r = c. Hence for all collisions of a pair of spheres with 
given V 

XRr = 2irm&V^p psin ^cos^ OdO 
Jo 

= ^7rmc^VY 

And for all collisions of a pair of spheres with whatever 
value of V, 

iRr = iTTC^pmV^ 

or \tRr = ^TTC^pT. 

But each collision is now counted twice, namely, once from 
the point of view of each of the spheres concerned. Counting 
each collision once, or from the point of view of one sphere only, 

\^Rr = ^iTC^pT (4), 

whence we obtain for each sphere 

and for all spheres 

and therefore \ ^IRr — JfTrc^p . pT. 

I shall write /c for ^7r(fp. If a be the radius of a sphere 
which on average contains one molecule, evidently 

or K is four times the aggregate volume of all the spheres in 
unit of volume. 

It follows from Art. 24 that 

or i> = |(l+Wpy (5). 

On a system of Elastic Spheres in Vertical Column. 

29. The following investigation depends upon the usual 
assumption made with respect to collisions between elastic 

2—2 
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bodies, that a certain momentum is instantaneously transferred 
through a finite space. The responsibility for that assumption 
rests on the inventor, whoever he was, of the theory of rigid 
elastic bodies. Consider a vertical column of gas, or material 
system, whose molecules are equal elastic spheres, under a 
constant vertical force/. 

If p be the pressure per unit of surface, x the height of a 
point in the column above a fixed plane, m the mass of a sphere, 
p the density, p^ the density at the base, we have 









1=-^^^' 




also by (5) 




P 


= f(i + i/.)pr. 




Whence, 


if we 


make 


3 

= constant = 77", 




we find 




-2m/A = ^'7^ + f7rc= 


dx' 



or /o = poe-^^'^-^^e- ^''-*«> . 

Again, consider N spheres crossing the plane ^ = with u 
for vertical component of velocity. Of these some, say N — N\ 
will reach the plane x=dx without collision. N' will undergo 
collision before reaching dx. But for these N' there will in 
stationary motion be substituted, as the result of collisions, N' 
other spheres with the same vertical component u. 

Now, if the impact were direct, i.e., the line of centres at 
collision vertical, the effect of the 
collision would be to substitute for a 
molecule with vertical velocity ii at 
height X a molecule with the same 
vertical velocity u at height x + c, 
where c is the diameter, or, as we 
may express it, the substituted sphere 
would gain a vertical height c, with- 
out losing in respect of that distance 
any kinetic energy to the force /. We might without affecting 
the general motion suppose the two spheres to change places 
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after collision, so that one would lose, and the other gain, the 

height c without change of kinetic eoergy. If, therefore, all the 

N' collisions were direct, the average height of the iV^ spheres 

doc 
or their substitutes at the end of the time — would be, not doo, but 

u ' 

N' , , 

dx-\- ^Cy while their loss of kinetic energy would be Nmfdx, 

But all impacts will not be direct ; we must consider then 
the result of indirect impacts. For this purpose consider two 
classes of collisions, in one of which the sphere A has vertical 
component u before collision, and in the other A' has vertical 
component u after collision. The effect of a pair of collisions one 
from each class, is to substitute A' for A as the sphere with verti- 
cal component u. Now let I (see figure) denote the vector line 
of centres at collision, and cos {ul) the cosine of the angle which 
I makes with the vertical. Then in the first of the pair of 
collisions the centre of A is below the point of contact by 
\c cos {ul). In the second, the centre of A' is above the point 
of contact by \c cos {id). There is no reason why the point of 
contact should be higher or lower in one case than in the other. 
It will be on average at the same height. Therefore on average 
of all pairs of collisions substituting A' for A with vertical 
velocity Uy A' i^ above A by 



c cos {ul) = r, suppose. 

Let q be the relative velocity of the two colliding spheres. 
Then considering q as radius vector of a cone of axis u^ and I as 
radius vector of a cone of axis q, we get 



r = c cos {uq) cos {ql) 



= ^c con {uq)y 

r- / — 

because cos {ql) = j cos^ sin dO P cos sin d0 = ^. 

Let V be the absolute velocity of the sphere whose vertical 
component of velocity is u, so that 

cos iuv) = - . 

V 



J 
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Then r = %g cos (iiq) 



= |c - cos {vq). 

Let i|r be the velocity of the other colliding sphere, E the 

angle between v and x/r. Then 

It V — yfr cos E 




:|C 



V 



We have to multiply this by the 

number of collisions which N 

spheres having velocity v undergo with spheres of velocity 

sjr . . . '\jr -\- dyjr, making with v angles E ... E -{- dE in time dt, or 

dijC 

, and then integrate for all values of ^jr and E. 



u 



Let pf('^)dyfr be the number of spheres in unit volume 
with velocity yfr ...yfr-^- d'^. The result is 

^r o /"° 7 . /-/ I \ f"" 1 • r^ Tri cy U V — 'dr cos E dx 

Nwc^pj^ dffWJ^ i sm EdEqic - Y^ ^ T 

= §7rc^/o . Ndx 
= KNdx. 

Therefore at time dt the average height of the N spheres or 
their successors above the plane x—Q is {1 + ic) dx. 

But the energy which allowing for substitutions they lose 
in the ascent is Nmfdx. The loss takes place only during 
free path. It follows that the loss of energy due to the ascent 
dx is, allowing for substitutions, m/(i^(l—/«:) per sphere. 

Now suppose that at ^ = the number per unit of volume 
of spheres having \mv? .,. d (u^) for energy of vertical velocity 
is 

Ke-^"^""' d (u^) (a), 

where K is constant. Then, by what has been proved above, 
the number which at height dx have ^mu^ ... d (u^) for energy 
of vertical velocity is 
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and the number which at height dx have 

^mu^ ..,d (u^) — mfdx (1 — /c) 
for energy of vertical velocity is 

^—2hmifdx^—9K J^£—h {mu'^—2infdx) (i— «) ) ^y^ 

that is, Ke-^'^'^'du^ (6), 

because, neglecting /c^, &c., 

and therefore 2hmffcdx = — Bk. 

The two groups (a) and (6) are equally numerous, and 
therefore either can by ascending or descending, allowing for 
substitutions, exactly replace the other. Now this is the 
reasoning by which in the ordinary case, when a: = 0, we prove 
p = p^e-'^'^^'^. It now proves p = p^e'^'^^'^e-'', 

I shall refer to this problem later for illustration. 

30. The pressure per unit of surface is increased in the 
ratio 1 : (1 + \ic) as the molecules, from being material points, 
become spheres with finite diameter c. 

But the pressure per unit of surface is the quantity of 
momentum which is carried through unit of surface in unit of 
time, Art. 12. Now, so far as this momentum is carried 
through the surface by molecules during their free path, it is 
not altered in the least by ic acquiring finite value. The 
increase of the transfer of momentum consists in the case of 
elastic bodies in the process above explained, namely, the 
instantaneous transfer of momentum through a finite distance 
which occurs on collision. 

The result of Art. 29 may also be explained thus. The 
mean effect of all collisions on any sphere is equivalent to a 
force k/ acting from the greater to the less density. 



CHAPTER III. 

THE DISTRIBUTION OF MOMENTA AMONG THE MOLECULES, 
ASSUMING CONDITION A. 

31. The simplest form of molecule that we can conceive is 
an elastic sphere, such as described in Arts. 28, 29. 

An elastic sphere, so conceived, may be considered as having 
three degrees of freedom only, namely, motion of translation 
in space. For the spheres being supposed perfectly smooth, 
motion of rotation about an axis, if non-existent, will not be 
produced, and if existing will not be altered, by collisions. 

Let u, V, IV be the component velocities of a molecule, m its 
mass. Then the molecules being elastic spheres of mass m, let 
the number per unit of volume of molecules whose velocities 
in direction w lie between u and u + du be proportional to 
e~f^'>nu^du. This is called Maxwell's distribution. It has been 
proved by several writers that assuming this distribution 
to exist, and assuming condition A (which is necessary), the 
distribution will not be disturbed by collisions between the 
spheres. 

My object in the present chapter is not to give another 
proof of this proposition, but to show in what manner condition 
A affects the truth, as well as the proof, of the proposition. I 
shall therefore follow the general method elaborated by Boltz- 
mann and Watson, as being more convenient for my purpose 
than other known proofs, e,g. Professor Tait's. 

32. Consider two sets of spheres, one set having mass M, 
the other mass m. Let the number per unit of volume of M 
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spheres be p^, of m spheres p^. Let us form a diagram of 
velocities, the velocity of any sphere of either set being repre- 
sented in magnitude and direction by a line drawn on the 
diagram from the origin to any point P or x, y, z. At P 
suppose an element of volume dxdydz. The number per unit 
of volume of spheres M whose velocities are represented by 
lines drawn from to points within that element shall be 

PmF{x, y, z) dxdydz, 
or shortly p^Fdxdydz, 

Call these the class F {dxdydz). 

Similarly the number per unit of volume of spheres m whose 
velocities are represented by lines drawn from to points 
within the element of volume d^drjd^ at the point ^, 77, f shall 

be pmf(^>V, K)d^d'nd};, 

or shortly Prnfd^di^d^, 

and we will call these the class 

We may also express the above definitions as follows. The 
chance that a sphere M shall belong at any instant to the class 
Fdxdydz is Fdxdydz. The chance that a sphere m shall belong 
at any instant to the class fd^drjd^ isfd^drjd^, 

33. Now assume condition A. 

Then the two chances are independent. And therefore the 
chance that a pair of spheres, M and m, shall belong, M to the 
class Fdxdydz, and m to the Glass fd^drjd^, is 

Ffd^drjd^dxdydz, 

On our assumption of condition A this is true even if the 
two spheres are on the point of collision, and therefore very near 
each other. For that case is not within the exception defined 
in Art. 17. If we assume B instead of A, it may be true only 
if the two spheres are at a considerable distance from each 
other. But in this chapter I assume condition A. On this 
assumption the number per unit of volume of pairs of spheres, 
M and m, which have velocities between the limits aforesaid, 
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and for which m lies within the element of volume dco, so situated 
that M and m are approaching collision, is Ffd^dr^d^dxdydzdw, 
If the element day be so situated that M and m are separating 
after collision^ the case is within the exception of Art. 17, and 
therefore condition A does not apply. Let us now define dco. 
Let, namely, E" -=^ {x - l^f ■{■ {y - rjf -^ {z - ^)\ so that R is the 
relative velocity of M and m. 

Let c be the sum of their radii. About the centre of M 
suppose a circular area described with radius c in a plane 
perpendicular to i?, and let 

day = TTC^Rdt. 

Then the number per unit of volume and time of collisions 
between spheres M of the class Fdxdydz, and spheres m of the 
class /dfdlTyd^, is 

pMpmFdxdydzfd^drid^iirc^R. 

And the whole number per unit of time of collisions which 
members of the class Fdxdydz undergo with spheres m of any 
class is 

PuPm^dxdydz 1 1 1 fd^drjd^Trc^R, 

Each of these diminishes the number of the class Fdxdydz 
by one. So the rate of diminution per unit of time of the class 
Fdxdydz due to collisions between its members and spheres m 



IS 



PMPmFdxdydz j ( j fd^d^d^irc'R (6). 



The Effect of Collisions. 

34. We have now to find the number of spheres M which 
by collisions with spheres m pass into the class Fdxdydz per 
unit of volume and time. Let P, p be the centres of M and m. 

On a collision the only thing that changes is the direction 
of jB. And that change depends upon (1) the angle yfr which 
the line of centres Pp at the instant of collision makes with J2, 
(2) the angle ^ which the plane through a line drawn through 
the centre of M parallel to R, say Pn, and the point of contact 
of M and m makes with a fixed plane through Pn, 
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The direction of R after collision is in the same plane 
through Pn, and makes the angle 2i/r with Pn. That is the 
easiest way to picture to ourselves the effect of a collision. 

In order that a collision may take place, it is necessary that 
a line ps drawn through p parallel to R should cut the circular 
area above mentioned. We may divide that circular area into 
elements whose type is 2c^ sin y^ cos yfrdyfrd<f>. Assuming con- 
dition A, the line ps is as likely to cut the circular area in any 
one point as in any other. Therefore the chance that the new 
direction of R shall make the angle 2^...2i/r + d-v^r with the 
original direction of R is proportional to 2 sin yfr cos ylrdsjr, that 
is to sin 2ylrdyfr, It follows that, assuming condition A, for any 
given direction of R before collision, all directions after collision 
are equally probable, and on average the whole kinetic energy 
is divided equally between the spheres. 

35. The direction of R after collision is then a function of 

^, y> z, f , 97, f, -i/r, </>. 

Let the new values of these variables be denoted by the 
corresponding accented letters. Evidently '^, <f> are unaltered, 

or (j)' = (j>, yjr^ = yjr. 

All those pairs of spheres for which before collision the 
variables a; .,. ^ are between 

00 ... 00 -{-dx ,,, ^ ... ^+ d^y 

and yjr, <p between limits 

yjr ... ylr^ dyjr, (jy . . . (f) + d<l> 

will after collision have these variables between limits denoted 
by the corresponding accented letters. 

Now let p j^F' dx dy' ds! be the number per unit of volume of 
M spheres which at any instant belong to the class F'dxdy'dz. 
Similarly after collision m passes into the class f'd^'dr]^d^\ 
whose numbers Siref^d^^di]'d^\ 

36. If after collision the velocities of M and m were 
reversed, without change of their positions, they would by a 
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reverse collision pass into the classes Fdxdydz and fd^drjd^ 
respectively but with reversed velocities. 

It follows that the number of pairs of spheres, M and m, 
which before collision belong to the classes F'dxdydz' and 
f'd^d'i)d^ respectively, and by collision with given values of 
^ and </> pass into the classes Fdxdydz and fd^drjd^ respectively, 
is, again assuming condition A, 

Pj^p^F'dxdy'dz'fd^'d'n'd^ 2 sin y\r cos fdylrd(f>c'R 

37. Now x\ y, z, ^\ 77', f ' are functions of x, y, z, ^, 9;, f, i/r, 
<^, and of those quantities only. Also the system of two spheres 
passes with unchanged kinetic energy from one state to the 
other. 

Whence it follows by the general proposition proved later, 
Art. 50, and can be proved independently, see Watson s Kinetic 
Theory of Oases, second edition, p. 9, that 

dxdy'dz'd^'d'q'd^' = dxdydzd^drjd^ (7). 

Therefore the number of pairs of spheres, M and m, which 
per unit of volume and time pass out of the classes F'dxdy'dz' 
and f'd^'drj'd^' into the classes Fdxdydz and fd^drjd^ respec- 
tively is with given yfr and <}> 

PmpmF' f'dxdydz d^dr]d^ 2 sin -^ cos yfrdyjrdcliC^R. 

And if now ^, 77, f, yfr, (f> vary, we find for the number of 
M spheres by which, owing to collisions between M and m, the 
class Fdxdydz is increased per unit of time the expression 

PMpmdxdydz [[( ^ d^drjd^ (^T 2 sin yfr cos y^dy^d<\>Fyd'R 

J J J -oc J oJ 

= pMpmdxdydz 1 1 1 d^drjd^FfiTC^R. 

But the number by which the same class is diminished per 
unit of time by collisions with the m's we found to be 



pmpmFdxdydz 1 1 j fd^dyd^irc^R. 
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dF 
If therefore -^ dxdydz denote the change per unit of time 

of Fdxdydz, due to collisions of M with m, 

^ dxdydz = dxdydz jjj d^dr^di; {F'f - Ff) ttc'R . . .(8). 
Similarly 
^^d^dvd^=d^dvd^ jjj dxdydz {F'f ~ Ff) ttc'R (8a). 

38. If Ff = Ff in all cases in which a pair of spheres M 
and m can by collision pass from the classes F and / to the 

classes F' and/' respectively, or vice versa, — = 0. Similarly 

4-.—0. And this being true for all values of F and /, the 

motion is stationary. Assuming then condition A, Ff = Ff is 
a sufficient condition for stationary motion. 

The proof that it is a necessary condition is due to Boltz- 
mann, and is as follows. 

Let ^"^IIJ ^ d^dydzF (log F-1) 

-^jjj^'y^dvd^/aogf-n 

^jjjd^dvd^^^logf 
And substituting for -^ and -j- their values above found, 

'^jt = 1111 ^^y^' wfy^^'^^^w^^'^' ~ ^^^ ^^""^ ^ + ^""^^^ ''^'^ 

= ^^^ dxdydz jjj d^dvd^ {F'f - Ff) log {Ff) ^c-R. . .(9). 
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But in the course of this integration x, y, z and f , t], ^ assume 
all possible positions on the diagram of velocities, and therefore 
assume the positions corresponding to any F' and /', and there- 
fore F and / among other values assume the values F' and /'. 
Whence 



dH 

dt 



= jjj dxdydz jjj d^drjd^ (Ff - F'f) log F'firc'R, . .(10). 



and adding together (9) and (10) 

^ = i jjj dxdydz jjj d^dvd^(F-f'-Ff)\og^,-K<fR... (11). 

And this is necessarily negative, if not zero, and then only 

zero when F'f = Ff for every case in which two spheres can 

pass by collision from the classes F and / to the classes F^ and 

dfT 
f respectively, or vice versa. But in stationary motion — y- 

must be zero. Therefore in stationary motion, given condition 
A, Fy = Ff, is a necessary condition. A distribution of 
velocities in which F'f = Ff for all cases in which a pair of 
spheres can pass by collision from the state F and / to the 
state F' and/' or vice versa shall be called the normal distri- 
bution. Any other distribution shall be called an abnormal 
distribution. The H theorem proves that when the distribution 

is abnormal --^- is negative. 

This is known as Boltzmann's minimum theorem, or the 
H theorem, and H is Boltzmanns minimum function, or, as he 
sometimes calls it, the Entropy function. The theorem is as 
here given, founded on, and depends for its existence upon, 
condition A, without which we cannot use the product Ff to 
express the frequency, or chance, of the simultaneous occurrence 
in collisions of the two states denoted by /and F, 

071 the H theorem. An objection considered. 

39. If when the entire system of elastic spheres has, with 
H diminishing, reached a certain state, all the velocities were 
simultaneously reversed without change of the coordinates, the 
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system would exactly retrace its course with H increasing. It 
is therefore possible, when a system is set in motion, for H to 
increase in that motion. This was at one time considered an 
objection to the theorem, because the theorem was supposed to 
prove that imder no possible circumstances could H increase. 
It however does not profess to prove that H always diminishes, 
but that, assuming a certain condition H must diminish on 
average. For the assumption underlies the whole proof as 
hitherto given. 

Now to meet this objection, consider two spheres, M 
having velocities w, y, z, and m having f, rj, f. Construct 
a spherical surface of radius c about the centre of M, and 
consider the element of volume 

2c^ sin yfr cosyjr dy^ d<j> Rdt = dco. 




If m be within dco, the two spheres are either (1) approaching 
collision, or (2) separating after collision, according as dco is on 
one or the other side of M, Condition A asserts that the 
number per unit of volume of pairs of spheres, M and m, 
belonging respectively to the classes F and /, which are ap- 
proaching a collision of that kind — i.e. a collision with the given 
yjr and <^ — is Ffdx,,,d^d(Oy see Art. 17 and note. But con- 
dition A does not assert that the number per unit of volume of 
such pairs which are separating after a collision of that kind is 
Ffdx,,,d^d(o, because that case is within the exception of 
Art. 17. In fact if x, y, z and ^, rj, f denote the precollision 
velocities of the separating pairs, x\ y\ &c. their velocities after 
collision, the separating pairs are the identical pairs which dt 
seconds ago had respectively velocities 

x,,,x-\- dx ?'•..§'+ c?f, 
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and were approaching a collision of the kind in question. The 
number per unit of volume of the separating pairs is therefore 
by our assumption of condition A 

Ff, doo...d^, dco, 

or which is the same thing, 

Ffdx\..,d^'dco, 
and not F'f'dw...d^'do). (We are assuming at present that 
Ff^F'f.) Now in the reversed course the pairs which in 
the original course were separating have become the approach- 
ing pairs, and vice versa. Therefore in the reversed course the 
number per unit of volume of pairs, M and m, which, belonging 
respectively to the classes F' and/', are approaching a collision 
of the kind in question, is Ffdx'...dl^'da),dind not F'f'dx'... 
d^' d(o, as according to condition A it should be. We see then 
that condition A is not satisfied in the reversed motion, and 
why not. 

The fact then that in the reversed motion H increases with 
the time is no objection to the H theorem as a mathematical 
proposition, because the reversed motion does not satisfy the 
condition on which the E theorem is based. 

40. We now see the necessity for the exception referred to 
in Art. 17. For if condition A continues to exist throughout 
the direct motion without that exception, it must also continue 
to exist throughout the reversed motion. Therefore by the H 

theorem -^ is negative or zero in both motions. But if it be 

negative, and differ from zero, in one of those motions, it must 

JTT 

be positive in the other. Therefore -^ must, if we assume 
condition A without exception, be zero in both motions. 

In other words, if we assume condition A with the ex- 
ception, the mathematical consequence of our assumption is 
that if is either minimum or diminishing, as the theorem asserts. 
If we assume condition A without the exception, the mathe- 
matical consequence of our assumption is that H is minimum 
already. This condition is satisfied by Maxwell's distribution. 
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41. We must now consider Boltzmann's own assumption 
that the motion is " molecular-ungeordnet." It being assumed 
that in the direct course the motion is molecular-ungeordnet, is 
it molecular-ungeordnet in the reverse course, or not ? I think 
Boltzmann's answer to this question would be in the negative, 
on the ground that we, or Maxwell's corps of demons, have 
" ordered " (geordnet) the motion by the very act of reversal. 
If this be so, then " molecular-ungeordnet " has, as applied to 
this theorem, precisely the same properties as my condition A. 
See Vorlesungen, p. 42. 

Let us endeavour to construct synthetically a system which 
shall without doubt be molecular-ungeordnet. The molecules 
being distinguished by numbers, I ask (say) Dr Watson to 
assign velocities to them according to any law he pleases. 
Then I, in complete ignorance of those assigned velocities, 
scatter the molecules at haphazard through space, and they 
shall start from the positions which I so give them with the 
velocities so assigned by Dr Watson. That is, prima facie, a 
molecular-ungeordnet system — in fact it is as near an approach 
to chaos as is possible in an imperfect world. 

Clearly also in this system condition A is satisfied, and the 
system will, or in all probability will, move from its initial 
position with H diminishing. Equally clear is it that if in the 
initial position all the velocities were reversed, condition A 
would be satisfied, and the system would, or in all probability 
would, move off in the opposite direction with H diminishing*. 
It follows then that if the selected distribution of velocities be 
abnormal, the state in which we have placed the system is one 
from which the system can and probably will move with H 
diminishing, but into which it must move with H increasing. 

* This statement is not inconsistent with that of Arts. 39 and 40 that — 

changes sign on reversal of the velocities. For in those articles the pairs of 
molecules having velocities x', ^', &c. are separating after an actual collision, 
before which they had velocities x, ^, &c. Corresponding to the classes x, ^, 
their number is Ff, In the synthetic system no collision has taken place, and 
the number is Fy. 

3 
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dH 

dt 



If the system were to move through that state, — ^ would be 



discontinuous. 

It may indeed, if the distribution of velocities happens to be 
the normal one, be in Maxwell's distribution. It cannot be one 
of a series of states through which the system passes in its 
approach to Maxwell's distribution, as the H theorem requires. 
The system then which I have supposed to be constructed is 
not, as appears to me, in Boltzmann's sense molecular-unge- 
ordnet. It should seem that the molecular ungeordnet state 
must cease to be such on reversal of the velocities. 

What the H theorem proves then is this, that the distri- 
bution of velocities expressed by the equation F'f = Ff is the 
only distribution which can be permanent consistently with the 
existence, and the continued existence, of condition A or its 
equivalent. It is assumed without proof, and as I hope to 
show in the succeeding chapters is not generally true, that 
condition A can and does continue to exist. 

The rate at which H diminishes has been calculated by 
Tait, Watson, and others for certain very interesting cases. See 
Watson's work above referred to. Art. 15. 

42. Boltzmann (see Berlin Sitzungsberichte, February 
1897) considers that we have in this behaviour of the system a 
true irreversible process. That term is generally applied only 
to processes involving dissipation of energy, that is conversion 
of some other form of energy into heat. If however heat, or 
temperature is to be explained as the kinetic energy of mole- 
cules or atoms moving according to usual dynamical laws, 
every process, even if it does involve conversion of other forms 
of energy into heat, is theoretically reversible ; would in fact be 
reversed if at any instant all the velocities were reversed, as we 
assumed them to be in discussing the H theorem. If that be 
so, by asserting that any process is irreversible we assert only 
that by no means within our power can we reverse it, because, 
that is, we cannot control individual molecules. It is in this 
sense that, as I understand, the diminution of H in accordance 
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with the theorem would be called by Boltzmann an irreversible 
process. 

Boltzmann also shows (VorlesungeUy pp. 58 — 60) that — H 
differs by a constant only from the logarithm of the chance of 
the system being in its actual state, and therefore as H 
diminishes in absolute value, the system passes from a less to 
a more probable state. See note at end of chapter. 

43. A solution of the equation 
F'f' = Ff, 
for the two sets of elastic spheres is 






/-p 



(12). 



That is found by treating the kinetic energy as the only 
thing which remains unchanged by the collision. It follows 
from this solution that for the same sphere x^ = y^ = z^, and 
that for different spheres Mx^ = m^^ &c., or the mean kinetic 
energy is the same for each sphere w^hatever its mass. 

44. The complete solution of the equation Ff= F'f\ as 
given by Boltzmann, VorlesungeUy p. 131, is 

where J. is a function of the coordinates, and 

Q = rn{(a--uy-\-(/3-vy + (y-wy\ (13), 

u, Vy w being constant velocities ; F has the corresponding value. 

In stationary motion ^. — ^> that is 



( da? * dt da dt du 



[Ae'^^^0. 



Suppose now the system to be in a field of external force 
whose potential is ^, e.g. the column of Art. 29. 

3—2 
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If no other forces act on the system than those derived from 
%, we have 

-(s-S)=-i <»). 

and therefore the equation becomes 

{cL-u)(^+2hA^e~^^ + Sic, = (15), 

the solution of which is A = e~^^^ (16). 

For instance, if the molecules be elastic spheres of infinitely 
small diameter, or material points, between which no collisions 
occur, there are no forces except those derived from ;)^, and we 
deduce by the above method, A = e~^^^. 

The expression now asserts that the chance of a molecule 
having velocities a . . . a + c?a &c., and being in a position where 
the potential of the external force is %...% + rf^, is 

It follows that the distribution of the velocities is the same for 
all values of %, that is in all parts of the system, and that the 
number of molecules per unit of volume is proportional to 
g-2^x. This agrees with the result of Art. 29, because in the 
case now supposed /c = 0. 

45. If the spheres have finite diameter c, so that collisions 
do occur, there are other forces acting besides those derived 
from 'x.y namely, the infinite forces assumed to act for infinitely 
short times during collisions. And it is no longer true that 

d . . dy 

and we cannot obtain the result 

except by ignoring the infinite forces altogether. It is true 
indeed that they act in equal and opposite pairs, but no 
solution can be satisfactory which takes no account of them. 
See further as to this problem post, Art. 103, 
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46. If the intermolecular forces are finite, and have a 
potential yjr, we have 

A{ _ \ — _ ^% ^'^ 
dt ' dx dx ' 

And our formula leads to the conclusion that 

A = €-^^^+^' (17), 

expressing now the fact that as before the distribution of 
velocities is on an average the same for any molecule, in 
whatever position it happens to be, not only with respect to 
the external field, but also with respect to other molecules — but 
that the number of molecules per unit of volume in a position 
in which the whole potential is % + i/r is on average e'^^^+'^l 
I see no reason why we should hesitate to accept this result, or 
how we can avoid accepting it, except that Dr Watson (in the 
work above referred to, pp. 70, 71) hesitates to accept it, 
though I do not understand him to deny its validity. He 
puts the case that, with finite intermolecular forces, a molecule 
P may have a sphere of action, within which may be many 
other molecules. Any individual molecule outside of the 
sphere of action exerts on P a force which is negligible. But 
the infinite number of molecules outside of the sphere of action 
may exert on P a sensible force. And he shows that this may 
be included in the fixed centre forces. But he refrains from 
drawing any conclusion for the molecules within the sphere of 
action. 



Boltzmann and Watsons generalisation, 

47. Our molecules regarded as elastic spheres have each 
only three degrees of freedom, that namely of motion of trans- 
lation. 

Let us now attribute to our molecules the most general 
form possible. 

Let a molecule m be a material system, whose position and 
state at any instant are defined by n generalised coordinates 
qi ... qm and the corresponding momenta^! -"Pn- 
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Let pm, be the number of molecules of the kind m per unit 
of volume. 

Let Pmfiqi '^'qnPi*"Pn)dqi,., dpn , 

or shortly pmfdqi . . . dpn , 

denote the number per unit of volume of molecules m for which 
the variables lie between the limits 

q, ... qi + dq, 

qn ... qn+dqn . . 

Pj ... pi + dpi ^^^' 



Pii ...pn+dpn^ 

or, as we shall express it, are in the state c. 

We may express this in a different way that the chance of 
an m molecule being in the state c is/dqi ... dpn* 

It is assumed that the molecule is, so far as its own 
internal forces are concerned, stable. That is that if it be now 
in the state c, and if no forces except its own internal forces 
act upon it, it will after the lapse of some time, greater or less, 
find itself in the state (c) again, except so far as the coordinates 
defining its position in space, which are included in qi ... qn , 
are changed by its motion of translation. It is also assumed 
that / will not be altered by reversing the sign of pi... pn 
simultaneously. And therefore / cannot contain odd powers or 
products of the ps. 

In the same region let there be another set of mole- 
cules M, each having r degrees of freedom, whose coordinates 
and momenta are denoted respectively by Qi ... Qr and Pi... Pr. 
And about these make the same assumption as above stated. 

Let pjij[ be the number of M molecules per unit of volume. 

Let p^F(Q, . . . QrPi ...Pr)dQ,... dP,., 

or shortly p^FdQi . . . dP^, 
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denote the number per unit of volume of molecules M for 
which the same variables lie between the limits 

Qr... Qr + dQr I /px 



or, as we may express it, are in the state C. 

48. Let (f){qi.., qnQi-" Qr) be a function of the coordinates 
of two molecules m and M, such that when (j> is positive there 
is no mutual action between the two molecules. 

Let us now make an assumption equivalent, with respect to 
the molecules now under consideration, to condition A. Let 

us assume namely that so long as (j) is positive and -~ negative, 

the two chances F and / are independent. And therefore the 
number per unit of volume of pairs of molecules, one belonging 
to the m's and the other to the Jl/'s, for which at any instant 
the variables lie, for m within the limits c, and for M within 
the limits C, is 

pmpM^fdqi ... dpndQ^,.. dPr. 

When <f> passes through zero the two molecules pass into 
their sphere of mutual action, and an encounter commences 
between them. The chances F and / generally cease to be 
independent. 

It may be the case that </>, having become negative, would 
if the two molecules are not influenced from without, remain 
negative for all time. The two molecules, that is, may, each 
remaining stable, form a permanent union revolving round 
one another in a closed orbit. I shall not consider that case, 
but assume that after a certain time r, ^ will again pass 
through zero and become positive. The encounter ceases when 
</) again becomes zero. All those pairs of molecules for which 

at a given instant <\> lies between zero and -^ dty -^ being 
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positive, will within the time dt after that instant commence, 
and as I am now assuming, will in due time afterwards com- 
plete, an encounter. The number per unit of volume and time 

of such encounters must therefore contain the factor —^ . 

dt 

49. But the two molecules have between them only n -f r 
degrees of freedom. We must therefore use (^ for one of the 
generalised coordinates defining the position of the double 

system. 

Let us substitute </> for qn. Then all those pairs of mole- 
cules, m and M, for which at any instant the variables other 
than qn lie within the limits c and C respectively, and for 

which at the same instant j> lies between zero and -~ dt, will 

(tv 

in the interval of time dt after that instant encounter one 
another. The number of such encounters which commence in 
time dt is therefore 

pmpM^fdqi . . . dqn-idp, . . . dp^dQ, . . . dPr -^ dt. 

As a consequence of the encounter, after the time r the 
variables will be found between the limits 

qi qi +dqi' 



for m, and 



g'n-i ... q'n-i + dq'n^i 
Pi Pi+dpi 

Pn Pn +dpn 



,(cO 



jL ^ . . . -t ^« "T' Ct-L'j' 



.(C) 



for if, and cf) between zero and -^ dt, and the encounter 

ceases as (f) again passes through zero. The pair of molecules 
pass out of their sphere of mutual action. 
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According to our notation the number per unit of volume 
of pairs of molecules, m and M, for which the variables lie 

between these limits c', C, with -^ positive, is 

prnpuFTMi ••• dq'n-idpi ... dpndQi ...dPr -^ dt 

Conversely if at the beginning of the time r the variables 

were in the limits c', C, d>\ —^ dt with reversed velocities, 

at 

they would at the end of the same time r be found within the 

limits c, C with reversed velocities. We may call these reversed 

encounters. The number of reversed encounters less the number 

of direct encounters per unit of volume and time is 

PmPuFTdq^ ... dq\.,dp^ ... dpr:dQ^ ... dP;^ dt 

- pmPM^fdqi . . . dqn-idpi . . . dpndQi . . . dPr -^ dt. 



50. It can now be proved (the system being conservative) 
that 

dqi ... dqn-idpi ... dpn'dQi ... dP/ -^ 

= dqi... dqn-idp^ . . . dpndQi ...dPr-^. 

For the accented variables g/, &c. are all functions of the un- 
accented variables and of the time t. And as t is the same for 
every system within the initial limits, we have 

dS , dS 
dq^r^'' d^r^''^"'^ 

where S is Hamilton's Principal function for the motion in 
question. Similarly 

dS dS 
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Now writing d<f> for -^ dt and dcf)' for -f dt, 

dqi . . . dq'n^^d^'dp^ . . . dpndQ^ , . . dP/ 

= dq; ... dq'^,d(l>'dQ,\.. dQ; x dp/ ... dpn'dP,',,. dP; 
= dq^\,. dq'n.id(l>'dQ,' ... dQr'dq, ... dqn-id(f>dQ,, ,,,dQr 

X 2 4 ^' #/ c?P/ 
"" rfg'i rf^a *" dQr ' 

where the last factor denotes the functional determinant 

dp^ dpi dpi 

dqi dqz "* dcf} ' 

dqi 

and since 

,_dS dp^ _ d^S . 
P' ~d^r d^i~dq,dqr ' 

the functional determinant may be written 

d'S d'S d'^S 



dqidqi dq^dq^, '" dq^d^' 

d'S 
dq^dq^ 

In the same way 
dqi,., dqn-id<}>dpi . . . dpndQi . . . dP,. 

= dqi ,.. dqn-id<j>dQi ... dQr x dq^ ... dqn-idcjy'dQi ... dQr 

dqi dq^'*'' 

and the functional determinant has, but for its sign, the same 
value as before. Hence the continued products of the differ- 
entials are numerically equal, that is 

dqi . . . dq'n-id^'dpi . . . dpn'dQi . . . dPr 

= dqi... dqn-id(^dpi . . . dpndQi . . . dPr (18), 
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and therefore the excess of the reverse over the direct encounters 
per unit of volume and time is 

pMpm {F'f -Ff)dq^,,. dqn-id(l>dp^ ...dpnX dQi ... dPr, 
or 

JfTjl 

-^ = pMpm dQ^.,. dPr JJKF'f -Ff)dqi,.. dqn-id<f>dpi ,..dpn; 
and is zero if in all cases F'f^ = Ff. 

A corresponding equation holds for -^ . See further for this 

proposition Watson's Kinetic Theory of Gases, second edition, 
Arts. 8, 10, 12, 13. 

51. As this property is true for all values of jP and /with 
the corresponding F' and/', we see that if Ff = Ff in all cases 
where a pair of molecules can by encounter pass out of the 
classes F, f into the classes F\ f respectively, or vice versa, 
we have a sufficient condition for stationary motion. 

The H theorem of Art. 38 can easily be adapted to prove 
that it is also a necessary condition. Further we may take 
for solution of the equation Ff = Ff, 

where T is the kinetic energy of the molecule M, and T' that of 
m ; and if T, T can be reduced to the sum of squares of the 
velocities, we get a result analogous to that obtained for elastic 
spheres in Art. 43, namely that, as a consequence of condition 
A, the mean kinetic energy is the same for each degree of 
freedom. 

52. It is necessary here to point out what further 
assumptions, besides condition A, we have made concerning 
our molecules in the above demonstration. We expressly 
assumed in Art. 47, that each molecule left to itself is stable, 
that is, that so far as the mutual forces between its parts are 
concerned, the coordinates and momenta, whatever values they 
have at this instant, will or may, at some future time pass 
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simultaneously through the same values or infinitely near to 
them, the coordinates x, y, z defining the position in space of 
the centre of inertia being excepted. We also expressly assumed 
that if two molecules pass into an encounter, they pass out of 
it again, that is do not form a permanent union. We have 
further assumed tacitly that each molecule is stable, not only 
under the influence of its own internal forces, but also under 
the influence of the new forces which may be developed during 
its encounter with another molecule. Suppose for instance for 
an m molecule /= er^^, where 5 is a quadratic function of the 
momenta pi.,.pn with coefficient functions of q^...qn' A 
necessary condition of stability is then that s shall be positive 
for all possible values of the p's. And that requires that the 
determinant of the coefficients in 5, which we will call d, and all 
its coaxial minors, shall be positive. See Appendix (g). 



If 



' — /If ^t_ 



^+&i2 2?lP2+&C., 



then 



d = 



dl ^12 ^13 
&12 0^2 hz 



Similarly in order that an M molecule shall be stable if left 
to itself, we have if 

F^e-^^ and 5f = ^ ^' + £12^1^2 + &c.. 



D. 



Ai -D12 -^13 
-O12 -"-2 -^23 



and D must be positive. And so long as the two molecules 
do not influence each other, the condition is that Dd shall be 
positive, which is necessarily true if it is true for d and D 
separately. But when the two molecules are within the sphere 
of their mutual action, the corresponding criterion of stability 
is that s-\-8 shall be positive whatever values the velocities may 
have. If A be the determinant of the coefficients in S -{- s, 
A, and all its coaxial minors must be positive. 
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But if we form the determinant A for the coefficients in 
s -{- S, it will generally during an encounter or when <^ < 
contain constituents of the form pP. The condition may fail 
for A, although it held when </> > for d and D separately. 
So that the two molecules, though stable separately, may be 
unstable together. The encounter may be the death of the 
molecule. 

Note on the R Theorem, Art. 42. 

The IT Theorem proves strictly that, assuming the independence of 

the chances (Art. 17), ~t~ is generally negative. It is claimed for it, as a 

physical theorem, that the diminution of JI is analogous to the increase of 
Entropy in irreversible processes. To make that analogy complete, the 

rJTT 

chance of -^ being positive ought to diminish indefinitely as time in- 
creases. I will state the difficulties that this view presents to me, without 
saying that they are insuperable. 

A finite number of molecules in a finite space are set in motion at an 
initial epoch at random, in the manner described in Art. 41, and then left 
for infinite time, undisturbed from without, under conservative forces. 

dH 
Initially -^ is almost certainly negative. After time t the coordinates 

and momenta of every molecule are determinate functions of t and of all 
the initial coordinates and momenta. If (7, C be two successive states of 
the system, each in its turn a consequence of the same initial state, in 
what sense is C less probable than C 1 It may be that the initial state, 
formed as we have formed it, is less likely to be C than C. But the 
chance of state C at time t is not proportional to the chance of the initial 
state being C. It is then improbable that the initial state should be one 

dH 
in which -^r is positive. It is not proved to be improbable that it should 

dH 

be one, a necessary consequence of which is that —r- will be positive at 

time t. 

Again, the continued independence of the chances, if conceded, 
attributes to every subsequent state of the system the same properties 
as the initial state has, which would be right if the system were con- 
tinually receiving external disturbances. Can it be conceded for the 
isolated system ? Let S be the Principal function, x, x' type coordinates 

of two molecules at time t. li x^af nearly, are the momenta -r-, -z-^, 

correlated or independent? I think it can be proved that they are 

d P 
correlated. If so, the formula 8 of Art. 37, expressing -^ , is inexact. 

dH 
And the theorem fails to prove -j- negative, except when the factors 

F'f -Ffd^re great, that is when ^deviates widely from its normal value. 



CHAPTER IV. 

CORRELATION OF VELOCITIES. FINITE FORCES. 

53. I PROPOSE to prove in this and the next chapter that 
in a system consisting of molecules of finite dimensions in 
stationary motion, it is not true for molecules very near to one 
another, that the chances of their having velocities between 
assigned limits are independent, as condition A assumes: but 
on the contrary if the forces be repulsive, they tend to 
move on average in the same direction ; and therefore that in 
dealing with their encounters, or (in the case of rigid bodies) 
their collisions, we cannot legitimately base our conclusions on 
condition A, except in the limiting case when the density — 
i.e. the ratio which the aggregate volume of all the molecules 
in unit of volume bears to unit of volume, — is infinitely small, 
though conclusions based on condition A are approximately 
accurate for gases under ordinary conditions. 

In the present chapter I shall deal only with the case of 
finite intermolecular forces. 

Let m, &c. be the masses of a number of molecules in 
stationary motion. Let x, y, z denote the coordinates of m, 
x\ y\ z those of any other molecule as m'. Let r be the 
distance from m to m\ Let / be a function of r which is 
everywhere positive, finite and continuous, which is of negative 
degree except for very small values of r, and such that 

^iTT^fdr is finite. 



/ 



Then for any m, 2/ is the sum of the values of /, referred 
to the position of m as centre, for all the other molecules ; 
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and if p be the number of molecules per unit of volume, 
S/ = I ^irr^pfdr. Also 2r -^ has corresponding meaning. 

Extension of Clausius' Theorem, 

54. We can now put the Clausian equation in a modified 
form as follows, 



dt 



H't^yt-'tiy]-" «■ 



for the function within brackets has a determinate value at 
every instant, and cannot on average increase or diminish with 
the time, the motion being stationary. 

Again for the same reason, 

d 
dt 



2m|sA' + 2^|W-HS^Ss/.j = 0...(II). 



Here the first 2 denotes that every molecule successively is 
to play the part of m, and the second S includes all the 
molecules except m each with its own x\ y\ z\ and each with 
its own/ referred to the position of m for the time being. 

Subtracting I. from II. we obtain 
^ dx^Ax' ^ dy^Ay' , ^ dz ^M' 

— %m 



{[^i-m-m^f 



+2-S^(^'-^)f=« <™)- 

A • df dfdr 



48 KINETIC THEORY OF GASES. [CHAP. 

And if X, yu,, V be the direction cosines of mm'y or r, drawn from 
m to m\ od — x = \t &c., and 

dr _ dec' dy' dz' 

dt~^irt^^~di'^''~di 

f^ dx dy dz\ 

and the last line of III. becomes 

s- s^ df f^ dx ^ dy ^ dz\ / dx dy' dz\ 

>^ ^ df f^ dx dy dzV 
-^'>'^'dr[^dt+^dt-^''di)- 

That is, since X/jl = \v = fiv = on average, and 

o 
on average, the last line of III. becomes. 

Again, -l^rf =-|J^ ^-rrr^prfdr 

= - i iirr'pf T" " + U-7rr''p/dr + (^ ^'^f% '^^ 
= If by the conditions for /, if -^ = 0.* 

Again, m-^-X,k,(i. 

if X, Y, Z be the components of all the forces, whether external 
or intermolecular, acting on m. 

* If the molecules be scattered uniformly 3^=0 on average. If not 

uniformly, it will be generally negative, because we take a molecule for centre, 
and there are more molecules where p is greater. But (see Art. 106) it will be 
small. If negative it helps my present argument, but I shall treat it as zero. 
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Therefore m -i— S/^' — x — X%f{oc - x). 

Making these substitutions, and arranging the terms, we find 

„ „j.(dxdx' dy dy' dz dz'\ 
'^^'^^J Kdt'dt'^ Jtlii^ didt) 
+ 'ZXS.fiaf -x) + ^Ytfiy' -y) + tZtfinf - z) 

K^ -c df /^ dx , dy , dz\ ,' dx' . dy' , dz'\ 






+ 2Z2/(«' - a;) + 2 F2/(2/' -y)+ XZ%f(z' - z) 
<• ■s' df [^ dx dy dz\ /^ dx dy' dz'\ 

■^^^^'i[^Tt+''i'^''dt){^-di-^^Tt'-''dt) 

= (IV). 

55. Now if there are no intermolecular forces, 

2Z2/(^' -x) + ^ Y%f{y' -y) + %Zlf{^ - z) 

is for an infinite system negligible, and equation IV. is satisfied, 
whatever /may be, by making 

dx dx' dy dy dz dz' _ 

J l^dx ^ dy dz\ (^ dx' dy' dz'\ 

on average for each value of r, and therefore for each value 
of /. This expresses condition A, which thus appears as a 
consequence of there being no intermolecular forces. 

B. i 
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56. But if there are intermolecular forces, we shall find 
that condition A cannot exist. 

In the limiting case when the encounters are binary, or no 
molecule is within the sphere of action of more than one other 
molecule at the same time, we have 

ZS/(*' -x)+ F2/(2/' - y) + Zlfiz -z) = - Rrf. 

In the general case consider a spherical shell of radii 
r ... r + dr described about m. Let jB be the intermolecular 
force between m and a molecule ra within the shell, 6 the angle 
made with x by the line mm drawn from m to m'. Let fr be 
the value of / for the radius r, and let 2^ denote summation 
for all molecules in within the shell. Then if Xy &c. relate 
only to the forces exerted on m by molecules within the shell, 

Xy{x -x)= -frRr {tr COS Of, 
the mean value of which, if the molecules be distributed through 
the shell at haphazard, is —^%.Br. 

Again X includes, not only the forces acting on m due to 
the molecules in the shell r . . . r + dr, but also all other inter- 
molecular forces acting on m, e.g. forces due to molecules 
within the shell r ... r -\- dr\ which forces w^e will denote by 
R'. The part of X^f{x — x) due to these forces is 

— frRr^r cos ^2r' COS 6\ 

where 6' corresponds to 9 above. If we assume 

2,.cos0S,.cose' = O (1), 

on average, we shall have 

But whether X^f{x—x) be equal to — ^^^Rrf ov not, it is 
nevertheless negative, and not zero. What has been proved 
for X is equally true for Y and for Z. So that 

XS/(^' - X) + FS/(2/' -y) + Ztf{z' - z) 
is in any case negative, and in the limiting case of binary 
encounters, or wherever we can assume (1), is equal to — ^Rrf, 

Since, as we have seen, in equation IV. the term 
tXtf{x' - ^) + &c. 
is negative and not zero, it follows that the sum of the remaining 
two terms, that is, 
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^ ^ ^ (dx dx^ dy dy' dz dz^\ 
^"^^J Kdi'dt^Ttli^dtTt) 

^ ^ df [^ dx dy dz\ / dx dy dz\ 

is positive and not zero. But this is inconsistent with condition 
A. It follows that we must discard condition A when inter- 
molecular forces exist. 

57. First solution of equation IV. Let for any value of r 
(dx^ dx' dy^ dy' dz ^ dz'\ 

"^[Tt^^dt-^di^^lt^dt ^'" dt) = ^^" ^^ ^^ ^^^^^^^^ 

where <^ is a function of r such that 47rr^0, and therefore also 

4!'7rr^(l>fy vanishes at either limit, when r is zero or infinite. Also 

dx dii 
assume all terms of the form -^ -^ , &c. to be zero on average. 

Then as in Art. 54, 

-:, df (^ dx dy dz\ [^ dx dy' dz'\ 



= ^t^ -j- , because V = ^^ = ^ = ^ 

= -j 4!rrr'p(f>fdr-^j iirr'p -£fdr, 
and using this in IV. we find 

2ZS/ (x -x) + &c. - iS j"%,rr> ^fdr = 0, 

which is satisfied by assuming for every r 

lr{X (x' - x)+ Y(y' -y) + Z (z -z)}==i4',rr^p'^; 

also I 4!7rr^p(f)dr = — ^ I 4f7rr^p -^ dr, 

Jo Jo ct"^ 

whence 4f7rr^p(bdr = ^Rr, 

• 

in all cases in which we can assume (1). 

Second solution of equation IV. Since / is arbitrary we 

may equate separately to zero the term in IV. involving r -^ , 

4—2 
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and the sum of the two terms involving f\ That is, we make 
for each value of r, 

_^ (^ dx dy dz\ /^ daf dy' dz'\ ^ .^xr x 
, [dx ^ dx' dii ^ dy' dz ^ dz 

^^^ '''W^^ll^<u^^'-dt-^dt^^'di 

+ ZS,. (^' - ^) + F2., (3/' - 2/) + ^S,. (/ - ^) = u. . .(IV b), 

in which X includes the intermolecular forces acting on m due 
to molecules at whatever distance, but x — x, &c. relate to a 
particular distance ?*. 

The result (IV b) and the first solution show that the 
velocities of m' molecules distant r from m are correlated 
with that of m, or 

(dx ^ dx dy ^ dy' dz .. dz'\ 

^ U ^" Jt + Tt -■■ di + di -'■ di) > ^' "" ■^'''^^'- 

And the correlation is a function of r. 

In this solution 47rr^p<^, and in the former one — ^4i7rr^p ~ , 

is a determinate part of the Virial of the intermolecular forces 
acting on m, though not necessarily equal to Rr, except in 
cases where (1) can be assumed. 

When r exceeds the " radius of action " of the molecule m, 
and so iJ = 0, X nevertheless includes the x force due to 
molecules distant r from m, for which E is not zero. But 
ultimately as r increases, S^. cos 6 S^' cos 6' — for all values of 
r' for which R is sensible. So that as r increases, 

m -T7 2r jT + &c. becomes in any case zero. 

But the correlation may possibly retain a sensible value at 
distances from m greater than its radius of action. 

Equation (IV a) expresses that although the velocity of m 
is correlated with that of m', yet it is not on average correlated 
with the velocity of m' resolved in direction 7nm', that is not 

correlated with X ^ + y^ -^ + i^ 4 • For (IV a) with (IV b) the 
dt dt dt 

simplest, though not perhaps the best, solution is the following. 
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If mP represent in magnitude and direction the velocity of 
m, and if it be taken for polar axis of a sphere of radius r 
described about vi as centre, then the average motion of 
molecules on the surface of that sphere is on the meridian 
towards P, and proportional to the sine of the colatitude. 

58. If the force R be attractive instead of being repulsive, 

doc dx 

'11 'T. -\- &c. will be negative instead of positive. The velocities 

of m and m' may be said in this case to be contrar elated. If 
there be two sets of molecules, A and By and every A repels 
every other A, every B repels every other 5, but A and B 
attract each other, the system will tend to form at every point 
two streams, the J.'s moving in one direction, and the J5's in the 
opposite, which we may compare, 

' si parva licet componere magnis,' 

with the mutually attracting and repelling electricities supposed 
to constitute an electric current. 

It should be noted that for any molecule m, 

fdx dx' 



, -. (dx dx o \ 



the summation including all values of r, is if (1) be assumed 
equal to the Virial, taken as positive when the forces are re- 
pulsive, of the intermolecular forces acting on m. (See Art. 85.) 



d v dx 
I have treated -j- —r- as zero in the absence of intermolecular 
dt dt 

forces. Strictly, n the number of molecules in the system being 

finite, and the centre of inertia at rest, it must be negative, but 

it may be neglected when n is great. 

59. I have thus proved that if there be finite intermolecular 
forces, condition A cannot prevail, and 

dx dx' dy dy dz dz' 
dt dt dt dt dt dt 
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has on average a finite value, a function of r, and is positive if 
the forces are repulsive. The investigation does not apply to 
elastic spheres, in which the intermolecular forces are infinite 
during collision. I shall therefore devote the two following 
chapters to proving the same thing for equal elastic spheres by 
an independent method. It will be, I regret to say, a much 
more laborious method. If however the reader is now satisfied 
that I have established my case for finite intermolecular forces, 
he will have little difficulty in provisionally taking for granted 
that the same thing can by some means or other be established 
for the limiting case of elastic spheres, in which the forces 
become infinite for very small values of r. And omitting 
Chapter V. proceed at once to Chapter vi. 



CHAPTER V. 

CORRELATION OF VELOCITIES. ELASTIC SPHERES. 

60. I PROPOSE in this chapter to prove for elastic spheres a 
proposition corresponding to that proved in the last chapter for 
molecules between which finite forces act. The elastic spheres 
shall be each of unit mass and diameter c. The same pro- 
position for spheres of unequal masses or diameters is much 
more laborious, and when proved does not bring us much 
nearer to any probable form of molecule. 



Equal Elastic Spheres, 

Assume the chance that a group of n spheres shall simul- 
taneously have component velocities 



«! . . . ai + dui \ 



7i ... 7i +dyi ) 



. for the first 



otn • • • ttn + da^ ' 

I3n . . . /3n + d/3n ' for the nth 

Jn ... yn + djn 

to be Ce-^'^da^.^dyn. 

If we assume condition A, we virtually assume that 
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Now it will be found that when the diameters are finite 
the motion is not stationary when 

but becomes stationary when 

Q = S (a^ + ^^ + 7^) + 2S6 (aa' + /3y8' -f 77'). 
Here a, a, &c., are component velocities at the instant con- 
sidered of any pair of spheres, and 6 is a function of the distance 
between them at that instant which, except for very small 
distances, is evanescent. I admit of course that given con- 
dition A always existing, the distribution of velocities among 
the spheres denoted by g-^^ (<i^+^^+y-) is not disturbed by collisions. 
That has been proved many times. But the point is that 
condition A itself cannot continue to exist when collisions 
take place. In order that it may continue to exist, it is 
necessary that the system be continually receiving disturb- 
ances, or else that the diameters of the spheres shall be 
infinitely small, so that no collisions occur. 

Definitions. 

61, Let P be any point, r the distance of any other point 
from P. Let / be a function of r satisfying the following con- 
ditions, namely, 

(1) /is equal to unity at all points within a small sphere 
of radius a described about P as centre. 

(2) /is finite, positive and continuous everywhere. 

(3) / is of negative degree at all points outside of the 'a 
sphere and such that the series denoted by 2/ or I Wr^fdr, 

J a 

is ultimately convergent, but contains no high powers of - . 

We may take for 'a' the radius of a sphere which on 
average contains one molecule. 

Let now ^= ^^ , 'n=^~yj^ ^^^f' summation in- 

cluding all the molecules of the system, and / having for each 
molecule the value proper to its distance from P. Then f, 7?, f 
are continuous functions of x^ y, z, the coordinates of P. 
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LetM.,=/J/dxdyd.a<y(f + 2) , in which (| + |)is the 

value of that function at a given instant within the element of 
volume dxdydz, aj is the value of the product ay for the 
molecule (if any) which is within that element, or the sum of 
those products for all molecules within the element if there be 
more than one, and the integration includes all space. Then, 
for any element of volume, aj changes from instant to instant 
by the passage of molecules into or out of the element. 

There are two analogous functions, namely, 

and also three others of the form 

M^^ = lljdwdydz a" ^ &c. 

It will be sufficient in this chapter to deal with Ma;z, and we 
may for the present drop the suffix, and write M for M^z- 



Summary of the Argument 

62. It is necessary first to treat of the case in which the 
molecules, though possessing finite mass, are of infinitely small 
dimensions, so that collisions between them will not occur. 
That being the case, the system will be in stationary motion if 
the velocities are distributed according to Maxwell's law, Q 
being a sum of squares only. This we may admit whatever 
views we may hold as to the possibility of any other law. I 
shall assume Maxwell's distribution to exist. It follows from 
Q being a sum of squares only, that the fact of any molecule, or 
any group of molecules having any given velocities, affords, if 
the whole number of molecules be infinite, no presumption 
whatever with regard to the velocities of any other molecules. 
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We then find for this system the mean values of f ^ [~T') > ^^• 

In this system of material points -7- = 0, which gives rise to 

equations A and B below. It is then proved that if, with the 
same masses, the spheres have finite diameter c, so that collisions 

occur, -f- , which was zero, becomes, if ^^ &c. retain the same 

mean values as before, positive, and the motion not stationary. 
But the motion must become stationary. Therefore in the 
system of finite spheres ^^ &c. cannot retain the same mean 
values as before. It follows that in this system Q in e"'^^ cannot 
be the sum of squares of the velocities only. 



A System of Material Points, 

63. The chance that the molecules shall at any instant 
have for their component velocities in x, Ui ... ai-\- dai, &c., is 
Ae'^^dui ••• dum where ul is a constant and in the system now 
considered 

Q = ai' + a/ + &c. 

But from Art. 61, 

ar/n = |S/- tti/i - . . . - dn-ifn-i- 

It conduces to simplicity to assume fn = 1, but, as will be 
seen, it in no way affects our conclusions. That being assumed, 

^n = ^2/ - ai/i - tta/a ~ &C. 

Substituting this value for a^ in Q, we obtain 

Q = f (2/)^ + (1 +/0 ai^ + (1 +/2O «2^ + &c. 
-2/;a,P/-2/;a,^2/-&c. 

+ S/i/saiOs + S/i/sttias + &c. 

If we now integrate 6~^^ for aj . . . an-i successively with limits 
± 00 , we get the result next stated. 

Let D be the determinant of the coefficients in Q, each of 
the coefficients of products being divided by 2, that is 
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D = 






Let Dpq be the minor formed by omitting the pih row and 
9'th column. Similarly let Dpgrs be the minor formed by 
omitting rows p and r, and columns q and s. Then by inte- 
grating for tti ... a^_i, Q is reduced to ^r- ^2, whence ^^ = ^t^. 

See Appendix ((/). Also by integrating for Wa • • • ofn-i only, Q is 
reduced to 

-'-'1221 -^1221 -^1221 

Again we find (see Appendix (g)) 

Ai = 1 +/i^ +A' +... f\-^ = SiV^ since A^ = 1, x 
A22l = S/^-/l^ I 

Du^=fJ^, \ (19). 

A2 = i>33 = &C.=i) = (2/)^ 

A2 = -/i2/: 

We have then for the chance that Wj and ^ shall lie within 
assigned limits irrespective of the values of Wa • • • «n-i the 
expression Ae'^^'da^d^ in which 



and 






=i> 






=p' 



and 



(s/y 



■ r. 



We see now that the assumption /^ = 1 has not affected the 
result, because writing ^ for/^ does not alter jp, q, or r. 

64. Let now ^ be the value of that function at P the 
origin, ^ + 3^ its value at a neighbouring point P\ Then at 

P', e"'^^' becomes e"'^^' + S^zt^ e"^^. And if ^ be very small, this 



is equal to 



e-^^{1^2hpa,d^). 
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But if at P, ^ is not negligible, we must use instead of 
a, /8, 7 the following functions, namely, 

u = a — ^, v = ^ - 7], w = ry — ^, whence a = u + ^ &c., 
and Q' = q (^ + ^y - 2p u-\-^ ^ + rp, 



or Q' = qu^ + 2q -p^u + (q-2p-\-r) ^\ 

and let e"''^' = <f>{u ^). 

Then if </> {u ^) be the value of this function at P, 0i {it ^) its 
value at P', 

<^i(w|)-0(«l) = 9?^</'(«a 

Again, let </> (i;) dli* be the chance that at P ^^ shall lie 
between u and u + c?t*, whatever f may be, that is 

J —00 

and the same chance at P' is 

Also |</,(^|) = (|)^(«^)-|-^(„|), 

if f^j denotes differentiation with u constant. Therefore 

But OK^)*^^"^^"^' 

because </> (i^ ^) vanishes at either limit. Therefore 

And this is proportional to 
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if h'-h ^^f^' 

whence <^i {u) = ^ (i^) (1 -f 2h'ud^). 

65. Now with the above definition of M we may write 

-Tr = 0, that IS 
at 

the summation being for all the elements of volume, each with 

dxj' 



its own a^ and its own ( "j^ + ~^ ' ■ 



Now we shall find that on average -^ a7 is of the opposite 

sign ^^ ( Tp + 3^ ) ) ^^d 771 ( j^ + 77^) is of the same sign with 

a7. So that —r- consists of two series neither of which is zero, 
at 

but which are equal and opposite. 

66. Let us consider the plane of wz, and | as varying 
in direction z, f as varying in direction 00, Let us first deal 

with -~ only. Consider the two infinite planes z = and 

z = dz, and suppose for a moment -~ constant between them. 

The number of molecules at or near the plane z—dz for which 
I*, or a — f, lies between ii and u-\-du exceeds the correspond- 
ing number at or near the plane ^ = by the quantity 

Ae'^'^^'dn^h'u^fdz. 
az 

That follows from Art. 64. To fix the ideas let u be positive. 
Let us separate from this class of molecules those whose z 
velocity is 7, and suppose 7 negative. Let us then call the 
separated class the class (^^7). Then the number of molecules 
of the class (uy) which pass through the plane z = dz per unit 
of area and time exceeds the number of the same class which 
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pass through the plane z = per unit of area and time by the 
quantity Ae"^'^''^ e-^y du2h'a<y-~dz. And the number of mole- 
cules of the class (1^7) which enter S the space between the 
planes per unit of time exceeds the number of the same class 
which pass out of S per unit of time by 

which is positive because 2^7 is negative. 

If SNdud^ be the number of the class {u^) within >S, we have 

^ dudy = ~ Ae-^'''' e-^y' dud<y2h'iirf ^ (20). 

But if uy be the mean value of 1^7 for all the molecules within 
S, u and 7 now taking all possible values, 

uy = 1 1 Nuydudy i I I Ndtcdy, 
and therefore 

because -7- I \Ndudy = on average, 

lie-h'^'e-^y'2h'uYdudy ^ by (20) 

2h dz ' 

I have assumed in this proposition that the space S within 

which -~ is constant, is the continuous space between two 

infinite parallel planes. That will not occur in fact. But it is 

not essential to the argument that 8 should have that shape, 

or that it should consist of one closed surface. The proposition 

holds equally if S consists of all those elements of volume for 

dP 
which at a given instant ~ has the same value. And it is true 

on average for each element of volume. 



= -^ 
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67. Again, _«^ = _(a-^)T, 

Now a molecule moving within S with positive 7 moves 
from less to greater ^, -j- being positive; if with negative 
7, it moves from greater to less ^. For all molecules moving 
within 8, -j 7^ is positive and equal to rf -^ , US be divided 
into parallelopipeds of base dwdy and altitude PP' 

^p and ^p' being the values of ^ at the upper and lower surfaces 
of the parallelepiped. By the motion of all these molecules 
y^ is increased. 

But at the upper end where ^ is greater, molecules with 
positive 7 pass out of the space, at the lower end they pass in. 
For negative 7 the converse is the case. By these entries and 

exits 7^ is therefore diminished at the rate 1 172 (^p - ^p) dxdy, 

or by the same quantity by which it is increased by the motion 
of the molecules within S. On the whole therefore if ^ alone 
vary 

., d— c?— d-r. 

d — 

We have then from Art. 66, 

d — __ 1 d| 
dt ""^ ~ 2hdz' 

Similarly, if ? only vary, ^^ ^ = - ^ g , 

and therefore if f and f both vary, 

d ,-, 1 fd^ , d^\ 
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and (A) becomes 



llh 



,dz dx) dt 2h \dz dx 



It follows that «7 ;7i ( ;; + j^ ) niust be on average positive, 



and equal to 2^2 + |y. 



68. We must however consider the value of -y- {~ + ~ 

dt \dz dx. 

more carefully thus : 

Let P, P' be two points on a sphere S whose coordinates 

referred to the centre C are respectively x, y, z, and Xy y, — z, so 

that PP^ is an edge of the parallelopiped, whose base is dxdy 

at X, y. Let ^p, ^p' be the values of ^ at P, P' respectively. 

dP 
Then if S be the volume of the sphere, and ~ the mean value 

^ dz 

d^ 
of -^ within the sphere, 

the integration dxdydz being throughout the space S. Let the 

dP 
radius of S be ' a ' (Art. 61). Then we may take -^ as equal to 

fata 

dz 

Suppose a molecule at G whose velocities are a, y8, 7. The 

dP 
change with the time of 8-^ due to the motion of that molecule 

^^ Vdx'^^dAj'^^dz)]] ^^^y ^^^ ~ ^^'^* "^^^ change due to 
the motion of all external molecules is zero on average. 

Let now / be the value of / referred to P as centre. Then 

^p = -^ . And let /' be the value of / referred to P' as 

centre, so that ^p' = ^^> . Then due to the motion of the 
molecule at C 

d ,^ ^ X ^ i d . o^ . d\(df df 
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df 



because on average -^S/=0 &c. The minus sign is here 



required because / is referred to P, and therefore -^ has the 

opposite sign to that which it would have if referred to C. 

Thatis, f=-?f. 
dz r dr 

Now to form equation (B) we shall have to multiply by aj. 
And that will make the terms in a, ^ disappear, because 
OL^oL<y = a/3a7 = 0, on average and therefore on average 

And therefore 

;^-ii-«vi.//«,(i-f), 

the integration being over so much of the plane of xy as is 
mthin S, that is " «t|, |=«y^^f, | being the value 

of -^ when r = a, 
dr 

By symmetry 

' dt dx ' ^fa da' 

It follows from (A) and (B) that in stationary motion of 
the medium of material points on average 

-^'vi if- 4(14^' <°)- 



Case of Spheres with Finite Diameters, 

69. We now pass to the case in which our molecules, 
instead of being material points incapable of colliding with each 
other, become equal elastic spheres of finite diameter c. It will 
be shown that as the result of their collisions with each other, 

—^ , which in the system of material points was zero, becomes, 
so long at least ^s f^ ( ~ j , &c. retain the same values as before, 

B. 5 
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positive, and the motion not stationary. It is worth while to 
show first in a general way by the aid of a diagram how this 
effect is produced by the collisions. 





Let us take the plane of the paper for the plane of xz and 
suppose f alone to vary in z, -^ being at a given instant 
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positive near a certain point. That being the case, there is at 
the instant in question near that point a mean motion, or 
stream, of the molecules in direction x increasing as z increases. 
We may suppose this instantaneous stream to be positive for 
positive values of z, negative for negative values of z. It is 
indicated in the diagram by the long arrows. 

Let us consider two kinds of collisions, (1) and (2), the 
directions of motion of either sphere before and after collision 
being shown by the small arrows. For simplicity these directions 
before collision are here drawn parallel to x. 

If we examine collision (1), we see that the molecule which 
before collision was moving towards positive x is deflected 
upwards, i.e, towards the positive stream, and the molecule 
which before collision was moving towards negative x is deflected 
towards the negative stream. In either case the effect is to 
increase the stream 'pro tanto. 

If we examine collision (2) we see that the reverse is the 
case. The effect of the collision is to diminish the stream. But 

now, -^ being positive, there are more collisions of the kind 

(1) than of the kind (2) per unit of volume and time. And 
therefore collisions tend on the whole to increase the stream. 

For the same reason, it will be found, they make —rr positive. 

It is true that in a rare medium the effect here indicated is 
very small. But the question is only of its sign, not of its 
magnitude. 

70. Since every collision changes the direction of motion 
of the colliding spheres, the expression in Art. 66 derived from 

/d^ ^ d^\ djay) 
\dz dxj dt 

contains, in addition to the result calculated in that article, 
a new term derived from collisions which we have now to 
calculate. 

5—2 
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Let V be the half relative velocity of two molecules, \, fju, v 
its direction cosines before collision, X', fi , v after collision. 
Let A,F= Vx, vV — Vz. The number of collisions per unit of 
volume and time for given V is 2ttc^p V. Hence, for the change 
with the time of VxVz due to collisions with given V, we have 

27rc^pV,V^{X'v'-Xp), 

the bar denoting mean values. We have now to find the 
values of X^u^ and \v for all collisions, given V. 

At the instant of collision let the centres of the two colliding 
molecules be ^, -S; so that AB is the line of centres. Let BD 
be the relative velocity. Or, if DE = BE, one molecule has 
velocity D£, and the other BE, in addition to the velocity of 
their common centre of inertia, whatever that may be. Let 
the angle DBA = 0, If A0 = BO, is the point of contact. 

Let BX, BZ be the directions of the axes of w and Zj and let 




the angle between the planes DBA and DBZ be </>, and the 
angle between the planes DBA and DBX be (j>\ 

Then, if oo, y, z be the coordinates of A referred to 
as origin, and if the direction DE be that of X, //,, v, 



f* n . 

z-= — v-^ COS Q -\--^n\ — v^^\wQ cos <^, 



iT = — X ^ cos ^ + s "^1 — ^^sin cos ^' 
The complete definition of <^ shall be this : — 



.(22). 
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When the plane DBA, turning round DB, contains the axis 
of ^, ^ = for that position of BA which makes the least angle 
with the axis of positive z. And 0' has a corresponding defi- 
nition for the axis of x (see Art. 72, post). 



71. If ^0, ?o be the values of ^ and f at 0, then at A 
^ - ^0 = ^ ( - 2^ cos ^ + Vl -z/2 gin cos ^) -1^ 



.f r (23). 



f~f,=-(-Xcos^ + Vl-Vsin^ COS </>'), 

Then by the same reasoning as employed in Art. 67 the 
number per unit of volume of pairs of molecules whose com- 
ponents of half relative velocity are 

is at A 

- ^JL^-^v- e-a.(V+r/) dV^dVydV, {1+4^'F»^- f„}. 

TT y 77" 

Here owing to the variation of ^ in Zj Vx is involved asym- 
metrically with respect to Vy and F^. Let us suppose now F 
constant, but the angles 6 and cf) to vary. The general mean 
value of VxVz is of course zero. But the mean value of V^Vg, 
given F, in the asymmetrical system is to be found by integra- 
ting the above expression according to 6 and <^, and is not zero. 
Let us denote the mean value of VxVz, or \vV'^, for all values 
of 6 and ^, F being constant, by (VxVz\ the mean value of the 
same function when Falso varies being denoted by V^Vz- We 
have then, treating ^ alone as varying, 

(VxVz) = — f*" 2 sin cos d0 f^dcp (1 -h 4A'F^ (|- ?o)) 

TT Jo Jo 

in which ^ — ^o '^^s ^'^^ value given in (23), and the general 
means of functions of V^ and F^, or \V and vV, are taken 
after integration according to and ^. That is, writing Vx/V 
for \ &c. 
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Jo 



+ 



VF^-4- F'' 
_j ^_ — ^ sm cos (f) 



Wc VJV: d^ 



(24), 



S V dz 

because the general mean value of VxVz, or \vV^, is zero. 

We have thus proved that, on the average of all colliding 
pairs of spheres, (VxVz) is before collision of the opposite 

7^ TO 

sign to -^ or ~- . This does not depend on the result of 

collision. It is therefore true if, instead of being elastic spheres, 
the molecules be centres of force, or whatever be their form. 

72. We have next to calculate the mean value for all 
collisions of XV, or the value of \v after collision. This will 
depend on the form of the molecule, and we shall now deal only 
with elastic spheres. It will be found that \V is of the same 
fd^ d^^ 



sign with ^^ + ^ 

We use the same notation as before. When the plane DBA 
turns round LB until it contains the axis of z, let Az, A^, which 
are on opposite sides of DB, be the two positions in which that 
plane is cut by the circle which A describes. Let Az be nearer 
to the positive axis of z than J./. Then, for Az, = 0; for 
Az, </> = TT. Similarly, when the plane turning round BB con- 
tains the axis of x, A^, A^ are the two positions of A in that 
plane, and, if J. a; be nearer than A^' to the positive axis of ^, <^' = 
for Ax, (j>' — TT for Ax^ 

Let 6 be the angle between the plane of DB and z and the 
plane of DB and oo, so that c^" = (^ — e, <^ = </>' + e. It will be 
seen that, if Xv is negative, cos e is positive ; if \v is positive, 
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COS € is negative. In either case 



cos e = 



— \v 



sine = 



f^ 



Vi-xWi-i/^' 



Vi - x^ Vi - 1/2 ' 

we have then, to find X and v\ 

V =-v cos 29 + Vl -1/2 sin 26 cos 0, 



;:} 



.(25). 



\'=-Xcos20 + V'l-Vsin2^cos^' 

Whence, substituting V^ for X'F, and F/ for v'F, we derive 
the two symmetrical systems 

F; = - F. cos 26* + VfTTIV sin 2^ cos ^ ' 



V,V, 



V^ = - F^ cos 26 - / ' sin 10 cos <^ 

V K ~ ~t~ V ,, 



F/ = -Kcos20- 



a; T^ ' y 

F.F, 

vf/ + v: 



V...(26). 



sin 2^ cos <^' 



F/=--F^cos2(9 + \/F/+ F/sin2^cos0'; 
The term involving sin € will disappear in the subsequent inte- 
gration, and is omitted. The above systems hold for all values 
of Fa; and F^. 

73. We now find that for given F, taking into considera- 
tion only the variation of ^, or -7^ , 

{V^V^) = - l^de rdct> 2 sin e cos 0(1 + ^h'V^f^,)V^V^, 

TTJo Jo 

F' F' 

But with the above values of X', i;', or -^ , ^ , 



f dl9rdI</>2sinl9cos(9F,'F/ = 0, 

Jo JO 



as is easily seen. Therefore 
1 (^'" 



(VafV;) = - rd0rdcl> 2 sine cos 4^h'Va,^-^oV^'V; 

TTJo Jo 

= - 4^h'c I dO rd^ sin cos 0V^ 

TT Jo Jo 

X |- ^^ cos g + ^ ^-' + ^"^ sin e cos J F; F/ g , 
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by (23). In this expression we have to substitute from (26) 
the values of V^f and V^, and integrate according to 6 and </>. 

To effect the integration according to <^, we reject odd 
powers of cf), and for cos^ (f> write ^ . We have then 

( V^' F/) = - ^ 4A'o ^de S^d6 sin 6 cos (9 K, 

-sin^COS<^l l^a;'^/ 



V, VF'4- F' . 



F 



: - ^ 4/^'c ['^6^^ f " c^d) sin cos dV^FOH, 
nrdz Jo Jo 



where 



n.... . VF,^ + F/ 



i^= j- ycos(9+^ y " sindcoscf) 

G = {-V, cos 2(9 + VF/+F/ sin 2(9 cos </>), 
H=(-V^ cos 2(9 - , ^^^^ -^ sin 2(9 cos cb) , 

( K,' F/) = - g 4/i,'c -^-y-^ .1 sin (9 cos^ 6> cos^ 2(9 d0 
+ ^ 4A'c -^V^ 1 sin (9 cos^ 6 sin^ 26^ (i(9 

- ^ 4/i'c ^- ^'^^"^"^^^ i sin2 e cos (9 sin 2(9cos 2(9d(9 

+ -^ 4h'c -^p' 1 sin^ e cos e sin 2(9 cos 2(9 d(9 
c?^ F Mo 

= -^4;i'c-V^ sin6>cos2 0cos2 2(9d(9 
c?^ F Jo 

+ ^ 4/i'c -^' 1 sin (9 cos^ (9 sin^ 2(9 dO 
dz F Mo 

- ^ Wc V'V\ f '%in2 6> cos (9 sin 2(9 cos 2(9 dd 

dz "^ ''Jo 

+ ^ 4/i'c -^-^ i sin^ (9 cos sin 2(9 cos 2^ d0 
dz F Mo 
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+ ^ Wc -^^ h sin^ COS e sin 2(9 cos 20 d0 
dz F ^Jo 

~rf^^^' F 3.5.7 +^^^^"^3.5.7 ^^ ^' 



but 



F.^ 1 v:v: 



V ""3' F4 3.5' 



or F"/F=5.-^-^. 

Hence the expression becomes 

"^^U^n ^"^^ -11 + 12-4 + 10 

1 d? 72^2 
""3.5-d5*^^~F~- 

By symmetry if we treat -^ as varying instead of -— , we 
obtain the expression 



1 ^u^y^^v: 



3.5 rfa; V ' 

and therefore, given V, 

And therefore comparing this with (24) 

(..■r,,-(nF.,=«'om(i.3L)(i.i) m. 

74. To find the rate of change with the time of Va^Vz 
given F, we multiply the last expression by 27rc2p F ( F being 
the half relative velocity), that is by the number of collisions 
which take place with given F per unit of volume and time. 

Let us denote the change due to collisions by ^ . 

Vv 

The result is, given F, 

|(F.r,.s*Wp(1.3lj)K,r.-(f.g). 
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And then we introduce the factor 

IT IT * u zi 

and integrate for all values of F^,, Vy, V^ and so we obtain 



/' 



75. This being the value of ^ VxVz, what is the value 

of|(«7)? 

If Ux, Uy, Uz be the components of the common velocity of 
the two spheres, a, ^, 7 the velocities of one, a', /8', 7' those of 
the other, 

Therefore ay + a '7' = 2 ( [^^ C7^ + F^ F^). 
Hence since Ux, Uy, Uz do not vary by collision, 

|(«7 + «V) = 2|Frn (28). 

It follows that on average of all collisions 

and therefore on average 

76. We have next to consider the change in the term 

a7-7- f-j^ + -7^ J due to the variation of a, or 7 by collisions. 

As before let us treat ^ only as varying. Let us use the 
construction of Art. 68. Consider two spheres in collision 
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their point of contact being at G the centre of the a sphere. 
The z ordinates of their centres are 
c 



^ (- z/ COS ^ + Vl — v^ sin 6 cos <^), 

and ^ (— i; cos ^ + Vl — v^ sin 6 cos <^). 

pi 
To find the value of ^ ^p due to the collision, let us refer/ 

to the point P. Its value at the point of contact being /, 
its value at the centres of the two spheres respectively is 

/+ ^ (- 1/ cos ^ + Vl - v^ sin 6 cos <\>) -~ , 
and /— ^ (— z' cos ^ + Vl — z/2 sin Q cos 0) -j- , 

with corresponding terms in -^ and ~ which will disappear 

on integration. 

Now for the first sphere a is increased by the collision by 
^x — ^xy for the second it is diminished by the same quantity. 
Therefore as the result of this collision 



9|p = 9 27 = - 2?^ j^(- '' c<^s + \/l-v^ sin 6 cos <^)(F^'-F,). 

For the minus sign here introduced see Art. 68. 

Substitute for V^ — V^ its value from (26), multiply by 
2 sin 6 cos 6, and integrate according to 6 and <^. That gives 

yd^p = - ^ ^ dJ y\^^de{^d<i> 2 sin e cos OFG, 

where F=^{- v cos 9-\- \/\ — v^ sin ^ cos <^}, 

a=[-\2 cos2 ^ - ^"^ ^ sin 2(9 cos <^) . 

That is tfd^p = " S ^ £^^^^- 

Now multiply by ^ird'pV, the number of collisions given V 
per unit of volume and time. That gives 

dt 15 >" ^^^"^fdz- 



dz' 
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Similarly t'= - ^^o^pF.F.^f , 

and |(^.-r.) = -^.o^.F.n^^(f-f). 

and ljjda:dy(^p-^p,) = -f^'^c^pV,V,jjdxdy~(^f^-f^ 

----TTcr'oFF — i^;S 
- u'^''P^''*''2fada'^' 

where the integrations are throughout the section of the ' a ' 
sphere made by the plane of xy, S is the volume of the a sphere, 

and - -T^ is the value of - -j^ when r = a. 
a da r dr 

But -^jj dxdy {^p - ^r) = 

as in Art. 68, therefore 

dtdz~ U'"^('^''^'lfada- 

, d (dP d^\ 32 „ „ ,^ 1 1 df 

This expresses the change oi -j--{- ~ due to collision between 

two spheres. Let a, /3, 7 be the velocities of one of the two 
colliding spheres, a', y8', 7 those of the other. 

Then ,— contains the two terms 
at 

a ^ (^ ^\ , , a /^ d^ 
dt\dz dx) ' dt\dz dx 

But oL^U^^V^ a'=U^^V^ 

as in Art. 75. 

Therefore (a7 + a'7') ^^Vz = 2F^^F/ on average, 
because UxUzV^Vz = on average. 
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Therefore for each sphere separately 

because ay = [7,=[7; + V^V^ + CT/F; + U^V; 
= 4i F/F/ on average. 

We have then on average, remembering Art. 68 (C), 

d (d^ d^\ 4 ,1 (d^ dn^ ,„, 

Adding this to Cj, we find for the term in -^ due to collisions, 



dt 



-i-v/JH^Ad+D" <«• 



We interpret this as follows. As the result of collisions, 
molecules moving in a given direction, as that of ^, are on 
average deflected towards that side where the average motion 
of the molecules, that is the stream, is for the time being in x. 
And so collisions tend to increase the stream. Compare Art. (57). 

77. Now referring to equation (A), we see that -^ , which 

in the medium of material points was zero, and so the motion 
stationary, has, when the spheres have finite diameter, become 
positive, namely 

2h\dz^ dx) 



, 2 . 
+ 3,rc. 



'4(1-0} i-i 



of which the first two terms are together zero, as in (A), and the 
third is positive. Therefore if in the system of spheres of finite 
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diameter f -7- j &c. have the same values as in the system of 
material points, —j- is positive. 

But in stationary motion —r- must be zero. Therefore in 

stationary motion of the system of finite spheres ( ~ 1 &c., 

cannot have the same mean values which they have in the 
system of material points. 

78. It might be suggested that the explanation of the 
difficulty consists in the fact that when the molecules have 
finite diameters the quantity of momentum transferred across 
any plane per unit of area and time is increased in a certain 
ratio, 1 : 1 -f ^/c, where, in the case of elastic spheres, k = fTrc^p. 
But this increases both terms in (A) in the same ratio, and 

therefore, so far as this property is concerned, -7— remains 

zero. This effect arises from the usual convention concerning 
elastic bodies, according to which there is ou every collision 
an instantaneous transfer of momentum through a certain space. 
It has nothing to do with the change of direction of the relative 
velocity, which also occurs on collision, and which we now have 
considered. The first effect would take place in precisely the 
same way, if after every collision the direction of the relative 
velocity were restored by Maxwell's corps of demons to what it 
was before collision. 



CHAPTER VI. 

CORRELATION OF VELOCITIES. ELASTIC SPHERES. 

79. The solution of the problem of the last chapter is 
this. In the system of spheres of finite diameter ^ has become 

That gives by arranging the terms in (D) 

^ = jjfd.dyd. [ay ^ (f + |) - ^(f + i)^ 

-^S P 2h LU^ dx) ^ U^ dx) J 2h [ dz-^ dx) J ' 

We may assume that «7 j^ ( j^ + j^ I is not altered by the 

substitution of f + r for f , or ay |(g + J) = 0. 

The first two terms of this equation are together zero 
by (A) Art. 65. The remaining two terms can be made zero by 
suitably choosing the ratio of the means 

\dz dx ) ' \dz dx) 
We have, namely, writing k for f Trc^p, 

X^dz dx) \dz dx))\dz dx) ' 
whence 

if fc^, fc^ &c, are negligible. 
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80. Instead of using M^^ = ijl dxdydzci'y (3^ + ;p) j ^^i' 
the analogous terms Mr,cy or Myz, we may make use of 

Mxx = I / Idxdydza? -^ , 

or the analogous terms Myy or M^z, and we obtain the same 
result. I have set out the mathematical work in the Appendix 
by reference to Arts. 66, 68. We may therefore now use the 
complete expression for M, namely 

\dy dx) ' ydz dxj 

+^^(I4J)} «. 

and we then reduce -^ to zero by making 

\dxj 1—K\dxj ' \dy) 1 — K\dy) 

Now so long as the distribution of velocities is represented 
by e~^^, and Q is the sum of squares of the velocities, ^ &c. and 

( 77 ) ^^* ^^^^ have either the same or less values than those 

found by the method of Art. 63, for material points. In 
stationary motion therefore in the system of finite spheres Q 
cannot be the sum of squares only. It must contain products 
of the form b (aot + /3/3' + 77'). And the question now is, what 
must be the form of the coefficients b, or as I shall call them 

coefficients of correlation, in order that the ratio I j ) ' i^) 
may have the required value. 

81. By way of illustration let us suppose n molecules within 
a sphere of radius R, their velocities being distributed according 
to Maxwell's law e~'^^, where Q is the sum of squares of the 
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velocities, and let us impart to each molecule, in addition to the 
velocity which it has in that motion, the small velocity X' in x. 

This being done, the chance that one of these molecules 
shall have velocity represented by lines drawn from the origin 
to some point or other within an element of volume dfl, whose 
centre is distant p from the origin in direction making the angle 
i|r with the axis of x, is 

Ae-^P' (1 + 2hpX' cos >/r) dn. 

The chance that another molecule shall have velocity repre- 
sented by lines drawn from the origin to some point or other 
with an element of volume dfl' distant q from the origin, in 
direction making the angle <^ with jp, is 

Ae-^^^dD.' [1 + 2hqX' (cos i/r cos </> + sin i/r sin </> cos e)}, 

where e is the angle between the plane of jp and q and the plane 
of p and X, The chance that both molecules shall have the 
velocities aforesaid is 

A^e-^ (P^+9^) dndD.' (1 + '2hpX' cos i/r) {1 + '^hqX' (cos y\r cos 

+ sin -^ sin ^ sine)}. 
Multiply this by \ sin i/r di/r, and integrate for all values of 
6 and y\r, and we obtain 

^2^-Mp2+^^ c^XlrfXl' ( 1 + ^' Xyq cos <^ , 

expressing now the chance that the velocities p and q, of two 
molecules, in whatever direction either is, shall make the angle 
<f) with each other. Now let u, v, w be the components of p ; 
u\ v\ w' those of q. Then 

dO^dDf = dudvdwdu' dv' dw\ 
pq cos ^ = uu' 4- vv' + ww\ 

Therefore, since, with these substitutions, the expression does 
not contain 0, it gives the chance that two molecules, both 
within the R sphere, shall respectively have velocities u.,M-\-du, 
&c., in the form 
^^^-H^+^+w^+u-^-vv'^^-w^) du.^.dw'U-v'^ X'^ {uu -h vv' + ww')\ , 
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that is, assuming X'^, &c., to be small, 

in which 

Q' = u^i-v^-^w'' + u^ + v"" +w''^-~^ X'-' (mi + vv ■+ ww'). 

Evidently, if we take Y' and Z' into account, we shall have to 

substitute X'^ + Y^ + Z''^ for X'^ in the last expression. 

We see then that in the case now treated the coefficient of 
correlation 6 has the form 

3 

or, since Z^ + Y' + Z' = ^ y . 

X, F, Z denoting the velocities of the centre of inertia of the n 
molecules, 

b = --(X''-^Y"' + Z'')l(X~'-\-T' + Z') (31). 

82. Next let us employ the converse method, namely, the 
coefficients b being supposed given, to find the ratio ^'- : f -. 
And first let us consider a finite number of molecules, and all 
the b coefficients equal, and very small, so that b\ ¥, &c , may 
be neglected. 

We have generally, referring to the determinant of Art. 63, 

P = ^l^ , and therefore 
zhlJ 

8p being the variation of ^ due to the introduction of the b 
coefficients. Now it is proved in the Appendix (n) that in this 
case dD = 0. Therefore 

^ ~ 2hl)~ 2hD • 
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Also A,-= 1+/,% /.A fj.,... 



AU /.A 1+//,- 

and (see Art. 83) d (1 +ff) = - 6/„ d (1 -irf^) = - 6/2, &c., 
dAf, = k(h-hf,-hf,),&Lo. 

If we were now to make every /= 1, we find that every 

constituent in Dn as it originally stood is in the varied form of 

2—6 
Dii multiplied by —^ . Hence 

But on the same assumption that every /= 1, 

2hD 2nh' 

1 n— 1 

And r = aA,/2/.i)=-2^-2-^' 



whence 


r^ 


+ v" 


'+r= 


3 « 




n being large, 


or 











as we obtained in a different way in Art. 81. 

83. Next, let us assume 6, the coefficient of correlation 
between the velocities of two molecules, to be a function of the 
distance r between the molecules to which it relates. And let 
it be required to find, on this hypothesis, the ratio 



r:Bo,(f)%(DV 



For the system of finite spheres we shall have 



6—2 
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where bpq is the value of b for the two molecules whose velocities 
are a^, Uq, &c., whence, substituting for «„ as in Art. 63, 

Q = (1 +,/;^ - i,„/0 «.^ + (1 +/? - b,nA) OLi + &C. 

-2a,^(l-J6,a/)-&c. 

+ 2 (//, + i&i, - -li.^/, - i6»/) «,«, + &c. . . .(.32). 

Then f + r^ = ^'L+|_^" = r + 26 "^l"/ 2A A 

r = S6f^/2AA 

every 6 being small, so that its square and higher powers may 
be neglected, and, for the same reason, the variation of D may 
be neglected, as will be easily found. Appendix (n). 
Now in this case [see Appendix (o)] 

in which p and q respectively take all values from 1 to n, there 
being n molecules in the system. Also if we assume the 
molecules to be scattered through space without regard to the 
position of p, with density p, XXhpgfpf^ may be represented by 
the intesfral 

Ifp 4^7rr'pdrbpgfq] 

J c 

and if bpq is very small except for very small vahies of Vpq (see 
Appendix {p)\ we find 

.CO .00 

If J ^irr^pbpqfqdr = lf; ^nrr'^pbpqdr (33), 

J c ■ c 

= 2//26, if 26 = 'i-jrr'pbdr, 

J C 

and therefore f ''^ = 9p 

= -i2/;26/2Ai) = -j26.p, 
because S// = Dn ; and if 6 be so chosen that y%b = — /c, 

Y' = Kf (34). 

84. Again, it can be shown that under the same conditions 

(f)=«(i' <->■ 

and this is the required ratio. 
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For since | = «! ^+ a2 ^+ &c., 

dP d fi d fo r. 

if ^.4|*«- 

1 d^ 0, e^ , 

whence ^n = -K- - ^ «! tt — <^2 ^ — &^c., 

Un ^Z "n "n 

and Q = ai' + a/+ ..- + «/ 

n 

Now the determinant of this expression differs from that of 
Art. 63 only in the form of the coefficients. Therefore 

W/ ~2hl)~2fi ■ 

Again, when the b coefficients are introduced, we have, 
substituting for 0^ in 6i»aia„, &c., 

( c\ ^1^2 7 1 ^2. 1 0\\ o 

and by the same process as before 

\dz) " \dz) '^hD'" ^h'^Kd'z) ' 
It follows from the above that any form of b which makes 

^S6 = — «:, or ^ / ^irr'^phdr = — ac, 
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and satisfies the equation (33), gives a solution of 






\dzj \dz, 

and the solution is independent of the form of the function/, 
so long as that function satisfies the conditions of Art. 61. 

The value of S6 is thus determinate. But the particular 
form of h which will make | / '^^irr'^phdr = — /c is not determi- 
nate. For instance, let 

That is a solution of (33) and (35) for all values of fi greater 
than 3. 

/// — 3 c^ 
Asfain, let 6 = - -„ when r < a, 

LL-3 c'af" , 
6 = — - — when r > a. 

That also is a solution for all values of fi greater than 3. 
Here a may be the radius of a sphere which on average contains 
one molecule, as in Art. 61. 

85. The law of distribution of velocities being e~^^, where 
Q = a^a^ -j- hi^cL-^a^ + a.^oii + &c., 
let a, /3, 7 be the velocities of a molecule of mass m, a', ^', 7' 
those of any other molecule. Then, h being small, it is easy 
to prove that 

mt (^ + ^^' + tV) = - S ^^, 

= - TS6, 

that is, ^m2 (oa' + W$' + 77 ) is equal to the Virial (taken as 
positive) of the intermolecular forces, represented in case of 
elastic spheres by collisions. In Art. 58, dealing with finite 
intermolecular forces, we advanced a step further, proving name- 
ly that the value of Jm (a? + ^' + 77') for any given r is equal 
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to a particular part of the Virial specially related to that value 
of r. If we are to make a similar step in advance in the case of 
elastic spheres^ it must be by dividing the Virial, if that be 
possible, into distinct portions, each specially related to a 
particular r. Instead of the Virial of the definite force which 
m' exerts on m at distance r, let us calculate the chance that 
m\ so situated, shall be the next sphere to collide with m, 
multiplied by the mean Virial of that collision if it occurs. 
That gives us a part of the whole Virial which stands in 
special relation to r. 

Consider a molecule m. The number of collisions between 
m and m with relative velocity V in time dt is proportional to 
irc^pVdt. Suppose the whole system to have attributed to it a 
velocity equal and opposite to that of m for the time being, so 
as to reduce m to rest. Then the sphere m' which collides with 
m with relative velocity V must have described a certain free 
path, which must have commenced at some distance r ..,r-\-dr 
from m. 

The chance that it shall have commenced at distance 

r.,.r + dr is of the form qe'^^dr, where q is independent 

of r, but a function of V, But w! starting from collision at 

distance r from m, is as likely to be moving in any direction as 

in any other; we must therefore divide e"^*" by 47rrl Again, 

2 Fc 
the Virial of the collision if it takes place is - m -r- by Art. 28. 

Let us then introducing the factor — ^ , assume for given V 
„ irc'V c _2 Vc 

and therefore, given V, 

/"* 2 mV^ /"^ 1 

^26 == I 4}7rr^pbdr = - irc^p ~-j- I -e-'t'dr 

2 , mV^ 

and taking means for all values of V 

i'^b = - ^TTc'pT = - kT. 
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To express 6 as a function of r, we must perform the 
integration according to F, and not according to r. That is 

2 3 47rr^ 4Jo 
in which g is a function of V, and G is the usual constant. 

That appears to satisfy all necessary conditions for b. But 
it is not a complete solution of the problem, because the intro- 

duction of the factor — ^ is not proved to be necessary, except 

to satisfy (33). It is possible that, as in Art. 57, the corre- 
lation, that is, b, may depend on the angles which the velocities 
of two molecules make with r, the line joining their positions, 
as well as on 7\ 



CHAPTER VII. 

FINITE INTERMOLECULAR FORCES. BINARY ENCOUNTERS. 

86. It is worth while to apply the method of Chapter v. to 
molecules between which finite forces act. Nothing is assumed 
concerning the force acting between two molecules, except (1) 
that it acts in the line joining their centres, (2) is a continuous 
function of the distance between their centres, (3) becomes 
repulsive and infinite as that distance is indefinitely diminished, 
(4) is evanescent at all distances greater than a certain line 
c which is very small compared with the dimensions of the 
space in which the system of molecules is moving. A sphere of 
radius c described about a molecule shall be called its sphere of 
action. 

It seems necessary to assume (3) that for sufficiently small 
distances the force becomes repulsive and infinite, if we are to 
hold that there is a very small distance within which two 
molecules cannot approach each other. 

87. With this definition of a molecule, let us assume that 
the number of molecules scattered through a given space is so 
small, or the aggregate volume of their spheres of action bears 
so small a ratio to the space, that no molecule is ever within 
the sphere of action of more than one other molecule at the 
same instant. In other words, the encounters are binary. Or 
that the number of such complex encounters which take place 
in unit of volume and time bears so small a ratio to the 
number of binary encounters, that we may without appreciable 
error assume all the encounters to be binary. 
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88. Suppose first a molecule P fixed in space. Let PX 
be any line through P. 

Let Q be a point just outside of the sphere of action of P 
and distant I from PX. And let another particle at Q of mass 
m move with velocity q parallel to XP. An encounter takes 




place, and since the encounters are assumed to be binary, 
m describes a symmetrical curved path QRQ', R being the 
apse. Let the angle XPR = i|r. The whole effect of the 
encounter, so far as regards the change of direction of the 
motion of m, is the same as if m received an impulse 2mg cos y^r 
in direction PR. Then also P experiences an impulse 2mgcos'^ 
in direction RP. 

Let the number per unit of volume of molecules which at 
or near Q have velocity q ...q ^-dqm direction XP be /(g) dq. 
Then the sum of the impulses which P experiences in unit of 
time from all the molecules which so move is 2mf{q) q^dq cos ^\r, 
and the mean force per unit of time on P due to all molecules 
moving from Q in direction XP is 



2m I /(g) q^ cos -^ dq. 



yjr being a function of g, and the integration including all values 
of q. If I also varies ^ is a function of q and I. 

Now suppose a circle described about X with radius Z in a 
plane at right angles to XP. Let m, moving parallel to XP, 
be before encounter anywhere on that circle, and let the angle 
made by the plane QXP with the plane of the paper be e. R 
will evidently lie on a circle parallel to the circle described 
about X. 
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If PZhe drawn, at right angles to PX in the plane of the 
paper, the mean force on P resolved in direction PZ due to the 
molecules moving from Q in direction parallel to XP is 

2m cos e I f(q) (f eos -^ sin '>^dq. 

And if now I also vary, the mean force on P in direction PZ 
due to all molecules moving in direction XP, whatever the 
value of /, is 

2m I cos edel ^irldl jf(q) q^ cos yfr sin ylrdq, 

i|r being now a function of I and q. It is of course zero. 

89. If P, instead of being at rest, be moving with velocity 
u in direction PX, we must write w + g for g in the above 
expressions. The force on P for any given value of e will 
then be in direction BP, and equal to 

2m / 27rldl I dqf{q) (q + uf cos yjr. 

And the mean force on P in direction PZ will be 

2m / cos ede 1 27rldl I f{q) (q + iif cos i|r sin y^dq. 

Now let the plane of the paper be the plane of XZ, XP the 

direction of x. Let the function f have at P the value Jo, and 

dP 
at any point whose z ordinate is dz the value Jo + x^^-s^- Then 

dP 
the velocity of P relative to Q will be g + ^ +;7^ ^ cos e. 

diZ 

And the mean force on P in direction PZ is now 
47^w^ I cos 6C?e 1 Idl \ dqf{q) ( g + w 4- i -^^ cos e j cos -^/r sin -^/r. 

That is, if u, v, w be the component velocities of P, 

dw /*"■ f^ f I dP \2 

-T. = — 47rm I cos ede I ZdZ I dqf{q) {q-\-u + l~-QOS€\ cosi/rsini/r 

= - 47rm I WZ I dqfiq) {q-\-u)l-^ cos i|r sin y^. 
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And u being negative let us write - u for u. Then 
dw 



""dt 



= + 47rm I Idl \ dqf{q) (qu + tt^) l-^coB^f sin i|r 
= + 47rm J Mljdqfiq) v? cos i/r sin -^ ^ , 
on average of all values of w ; or writing ^ for mii^ 



2h 



dw _^ d^ 1 
dt ~ dz2h 



where k = 4<7r I Ml I dqf(q) cos -i/r sin ^Ir. 

By symmetry 'Z/;^ has the same value mutatis mutandis, 

and therefore -y^uw = k \ ^ + ^] , 

dt \dz doo) 

(S-i)j,«"='(i-S)" «^)- 

90. This agrees in form with the result obtained for 
elastic spheres in Art. 76. But k, instead of being a known 
constant, now depends on the unknown law of force and is a 
function of cf". Further, in the expression for k, cos-v/r is 

TT 

necessarily positive, because i/r lies between zero and - , and 

Li 

also cos i/r sin -v^r is necessarily positive, and k is necessarily 

TT 

positive. Now if i/r > — , cos ^ sin -^ diminishes as -^ increases ; 

that is as q increases, because -~- is positive ; if -^/r < - , cos -v/r 
sin -^/r increases as q increases. 

TT 

Also the condition that -vlr shall be less than — is that I be 
^ 4 

less than a certain magnitude, which we may call l^, and the 

greater q is the smaller must Zq be. It follows that as g'^ or 

the mean kinetic energy, increases, it becomes less probable 

TT 

that ->^ shall be less than - , and therefore more probable that 
COS -^/r sin -v/r shall diminish as q increases. 
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Therefore if c^ be great enough, cos i|r sin i|r taken on average 
of all values of I and g', will ultimately diminish as (^ increases. 
And therefore k must for sufficiently high values of 5^ that is 
of the mean kinetic energy, or temperature, if the two are 
identical, either diminish as the temperature increases, or in- 
crease less rapidly than the temperature. 

91. We have now to consider the other term in -y^ , 

at 

namely, 

' dt\dz dxj* 

Let us again assume the encounters to be binary. As in 
Art. 68 describe a sphere of radius a about any point 0, and let 
PP' be a double ordinate parallel to z, P and P' being on the 
surface of the sphere. 



Then g = |j"ci*rf2/(^p-f^)^ 






Suppose two molecules to encounter one another, the point 
of contact of their spheres of action being at the centre of the 
sphere. Let F^;, Vy, V^ be the components of their relative 
velocity before encounter. We may suppose their common centre 
of inertia at rest, since its motion will not affect |p on average. 

d§ 
Let us consider the change of -^ due to this encounter 

only, that is 



//"Kf-t)*"?' 



d^ 
First to find -y^ , let / be referred to P as centre, and let 
at 

/o denote its value at the centre of inertia of the two spheres of 

action. Let / be the value of / at the centre of either sphere 

when the encounter begins. After the encounter has ceased / 

will have become/', corresponding to the final relative position 

of the two spheres of action. Let V.^, V^ have become F^;', F/. 
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If 8fp be the change in ^p due to the encounter, d^p is, 
omitting the factor =-. which is on average constant, propor- 
tional to, and has the same sign as 

VJf'-V,f, 
that is, if we write Vx = \V, VJ =X'V, &c., the same sign as 

x'f-V- 

Now let be the angle between the relative velocity, whose 
direction cosines are X, fi, v, and the line of centres at the 
commencement of encounter. 

Then 2yfr — 6 is the corresponding angle at the end of the 
encounter, yfr being the angle so denoted in Art. 88. 

Then we have, ^ being the radius of the sphere of action, 
/=/,+ |(-^cos0+Vr:^sin^cos(^)£, 

y'==y; -f- 1 (- i/ cos 2f^e-h Vf-^sin 2yfr^^0 cos (f>) -£] 

also as shown in Art. 72, 

\v . ^ 

X' = — \ cos 2i|r Sm Sl/r COS <^, 

Vl — V^ 

and therefore, omitting cos ^ and writing ^tt for cos^ </> as the 
result of integrating according to </>, 

xy ' = X\f, H- 1 Xi. . (cos 2^^r cos 2^-^~ ~ 4 sin 2^;r sin 2^J^^) -^ 

:=xy; + |xi^. ^{cos0(4cos^^-4cos^^ + l) 
— cos ^ 2 cos^ 'v/r sin^ -^ 

+ sin 2 sin -v/r cos i|r (2 COS^ 'i|r ~ 1) 

+ sin 6 sin -x/r cos '\/r (2 cos^ -^ — 1) j. 
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Also X/= X/o 

-\-^^(—v COS 6 + \/i — v^ sin 6 cos <f>) -f- 

= X/i — Xi' s cos ^ ~ , 
on average. Also for the pair of spheres (X' — X)/o = 0, and 
^ ^"2 cos (9 sin (9c?(9 (X/ - X/) 

'0 



+ 



f 

Jo 
Xi; I J^ ['" ^ *^^^'^ ^^^ ^ ^^ ^^ ^^^' '^~^ ^^^' ^ + 2) 
Xi/|^p2sin2(9cos^d(93sin>|rcosA/r(2cos2'i/r-l)...(36), 
which can be easily shown to be positive, being less than yjr, 
and yjr less than ^tt. Therefore °^7 ;7: ( 3^ + ^ ) i^ positive. 

92. If instead of being repulsive as hitherto assumed, the 
force between the molecules P and m be attractive, it is possible 
that the approach of the two may result in a permanent union, 
each molecule revolving in a closed curve round the common 
centre of inertia. Leaving this case out of consideration, the 
form of the path described by m with reference to P regarded 
as fixed will for attractive force be as in this figure. 

Q 




The angle yjr is now between ^ and tt instead of between 

zero and ^ as in the case of the repulsive force. The impulse 

on m is represented as before by 2mq cos yjr, cos yfr being now 
negative. The mean force on P due to a succession of en- 
counters is in direction PR instead of RP. The reasoning of 
Art. 89 will apply equally well to this case, except that in the 
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result 1^-77 will be of the opposite sign to -^ instead of the 



same 



T (dP d^\ d 
sign, and I -y- + ;/~. ) ;7t ^'^ will be negative. 



That leads in stationary motion to a positive sign for the 
coefficients h. It is conceivable that the force between two 
molecules, being a function of the distance r between them, 
may change sign at some value of r. If how^ever there be a 
limit of distance within which the centres of inertia of two 
molecules cannot approach each other, the force must ulti- 
mately, as r is indefinitely diminished, become repulsive. And 
as the density increases the near approaches becoming more 
frequent, the repulsive forces must ultimately predominate. 
The effect of the attractive forces, if such exist, will be, not to 
make the 6's positive, but to make them have less negative 
value than they otherwise would have. 

93. So far we have assumed the encounters to be binary, 
so that for each encounter either molecule describes a complete 
symmetrical orbit w^ith reference to the other, complete that is 
in the sense that after the two have separated to the very small 
distance c from each other, the motion of either is sensibly 
rectilinear in the asymptote. If we remove the restriction to 
binary encounters, we may reason as follows. 

d^ 
Assume -^ to be positive at the origin P. Then a molecule 

at P moving with ii negative has greater velocity relative to 
molecules on the side of positive z than to those on the side of 
negative z. On the average therefore the potential of the 
mutual action of P and the other molecules will be greater 
with z positive than with z negative. Let p be that potential. 

Then ~ is positive. 

Therefore if u, v, w be the component velocities of P, 
—- = — -r^ is on average negative. And therefore u ~, is on 
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average positive, or has the same sign with ~ . Also its sign 

has no relation to -^ . 
ax 

Similarly ^^ -^ is on average of the same sign as -y^ , and 

d^ 
independent of the sign of -p . And therefore 

[d^,d^\d — _ /d| d^\/dw du\ 
\dj2 dx) dt \dz dx)\ dt dt) 

is necessarily positive. 

And therefore in a system of molecules of the kind now 

considered -^ is positive, and the motion not stationary if 

^^, rf, ^ of Art, 61 have the same mean values as they would 
have in the medium of material points. To effect stationary 
motion, f 2, rj^y i^ must be increased as we found to be the case 
with elastic spheres. 



CHAPTEE VIII. 

GENERAL THEORY OF THE STATIONARY MOTION. 

94. Having established the results contained in Chapters 
IV. and v., I shall now assume for the law of distribution of the 
coordinates and velocities among an infinite number of mole- 
cules in stationary motion the following, namely : The chance 
that at any instant the coordinates of the molecules, n in 
number, shall lie respectively between the limits 

^1 and Xi + dx^, 

Vi „ Vi + dy^. 



Oury^ an CI Xiy^ "f" (XlJi"Jl , 

Vn » Vn + dyn, 

and their component velocities between the limits 
«! and «! + doLi > 

7i » 7i + djiy 



On and an -^ dan, 

In » "^n-^d^n, 



is 



(7 V 2 
- j VD^ e-'^^ dx^ . . . dznda^ . . . ^7^ . 
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Here Q = a, (a,' + A^ + 71^) + a, (a,' + ^i + 72O + &c. 

+ 612 (aitta + /3i/32 + 7172) + &c. 

The ' a ' coefficients are independent of the positions of the 
molecules. Every &, as hpq, is a function of the distance, r^q, at 
the instant in question between the molecules to whose 
velocities the suffixes pq relate, such functions satisfying the 
conditions of Chap. v. B is the determinant of the coefficients 
inQ. 

Again, A may be explicitly a function of the coordinates 
X, y, z, &c., or may be a constant, VZ) is a function of the 
coordinates, only as contained in the coefficients h. Then 
evidently 

n 

Instead of using the word chance^ we may say that the 
above expression represents the time during which on the 
average of any very long time the coordinates and component 
velocities will be found respectively within the limits aforesaid. 

95. The function Q may conveniently be put in the form 

Q = a^ai + i6i2aia2 + khzOL^c^^ + &c. 
+ a^oL^ + \h1a2OL1 + Ihc^a^a^ + &c. 
+ &c., 
in which each line is appropriated to a particular molecule. 

96. According to the results obtained in Chapters iv. 
— VI. it is necessary for stationary motion that ^'^ + rj^^ + ^'^ 
(Art. 79) shall have a certain value. That necessary value it 
will have if suitable values be given to the coefficients h as 
functions of r. The coefficients h being so determined, the 
motion is stationary, so far as not to be disturbed by the 
collisions or mutual forces of the molecules. I assume now 
that in Q the b coefficients have the values so determined. 

7—2 
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97. We can now deduce the fundamental property of this 
distribution, namely, that for any given set of positions of the 
molecules, or as we shall call it, for any given configuration, 

dQ dQ . 

ofj 7— = tta -J— = <^c., on average. 



For -j^ = 2(Xiai + h^^a^ + tisOts 4- &c. 



dQ^ 
doL^ 

If, with «! constant, ttg, «3 &c., assume all possible values 
consistent with the conservation of energy, we find (see 
Appendix Qi)) that the mean value of 

-7^ , or 2aiai + h^^a^ 4- &c., 

is J- ai, D being the determinant of the coefficients a, b, &c. 

in Q, and Du its coaxial minor obtained by striking out the row 
and column containing ai. Therefore also 
dQ^D^ 

But it is also shown (Appendix (g)) that 





<- 


~hD' 




whence on average 








dQ 


D J9„ 
DnhD 


1 
'"h 


Similarly 


dQ 
"'da,'- 


n 




and 




dQ_ 


&( 



98. Inasmuch as the system is infinite, and there is no direct 
action between any molecule and any other molecule except at 
distances very small compared with the dimensions of the 
system, not only is the sum of the potential and kinetic 
energies constant, but the kinetic energy of the whole system 
is itself separately constant. Let the kinetic energy be 
T = J m, (a,' + /5,2 + y,') + i 771, (a,' + ^,' + 72') + &c. 
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Here if ai be velocity of translation m is mass, if a be 
velocity of rotation m denotes moment of inertia, and so on. 
Then we have by the constancy of f, 

mifli "^ + '^202-^ +&C. =0 11. 

For the same reason Q is separately constant, or 

f«^ + ^^ + &c. = III. 

dui dt d(X2 dt 

Now II. and III. suggest the inference 

§j '""§.h'"-'^ '"■ 

An obvious solution of this is Q = ^^m {a^ 4-/3^ + 7^). But, 
as proved in Chapters iv. and v., this does not give stationary 
motion when the molecules have finite dimensions, except in 
the limiting case of infinite rarity. 

99. Concerning the Maxwell-Boltzmann Law 
rriiOLi = m^OL^ = &c. 

This follows at once from Art. 97, and IV., if IV. be accepted. 
It must however be noted that IV., although a consistent, is not 
a necessary, consequence of II. and III. For while ai, ota, &c. are 
unchanged, let all the molecules undergo small displacements. 

Then -^ , -~ &c., depending as they do on the intermolecular 

Cut Ojv 

forces, will receive variations, which may be arbitrary, as the 
displacements are arbitrary. If under these circumstances 

3-^ , &c. remained unchanged, IV. would be a necessary conse- 

quence of II. and III. But -j^ &c. do not remain unchanged. 

So IV. is not proved by II. and III. 

If we assume condition A, we assume in effect that the 
h coefficients are all zero, and Q = %md?. It follows then at 
once from the law of distribution e~^^ that m^a^ = 'm^^^ = &c. 
The law is therefore true whenever condition A can legitimately 
be assumed. It is true, that is, in the limiting case of infinite 
rarity. 
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I am not aware that it has ever been proved in any other 
case, or by any assumption not equivalent to condition A. 

If it were true in the general case, we should have, comparing 
11. with III, 

dQ j dQI . 

and therefore by Art. 97, 

miDu mJD^ 
or ^1 Ai = '^2^2 = &c. 

But Dii B22, &c. are functions, not only of the masses of the 
molecules rrii mg, &c., but (if they be elastic spheres) also of 
their diameters, (if they be centres of force) of their effective 
volumes. It seems therefore to follow that the law 

mitti^ = m2a2^ &c. 

cannot hold universally. It can be accepted only on the 
authority of the great physicists by whose name it is known. 

100. If therefore we write 
dQ / 

fiiy /JL2, &c. will in general have different values for different 
kinds of molecules in our system. Only in the case where the 

molecules are all of the same kind, may we write ~ I m = fia, 

where fju is constant. 

In all cases however in which the b coefficients are very 
small fii, /ii2, &c. will be very approximately equal to each other. 
For in the limiting case when the 6's are all zero 

Q = diOCi^ + agCta^ + &c., 
and we know that ai = mj, ag = ma, &c. When the 6's are not 
zero, every coaxial minor, as i)n, consists of the product of all 
the axial constituents mgrng . . . , and of other terms each of which 
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contains the product of not less than two 6's. The quantities 
m^aiy m^cL^y &c. will therefore, as long as the 6's are small, differ 
from each other only by small quantities of the second order. 

Field of external force. 

101. We may now equate to zero the time differential of 
A^De~^^y and so obtain, applying Boltzm ami's general method, 



\dt dx dt dy dt dz) 



i- s=^. §- 

Again let the forces acting on the system have a potential x- 
Then 

dai __ dx 

doi2^__dx n 
^ dt dx^ 



And the equation becomes, arranging the terms, 

dA^hA 
dx m 

V dz m dz dr^J 



dx da J \ dy m dy dp/ 



d^ J 

^ 1 ^^ dDdb . 
2 VD db dx 



Let us first consider the terms which involve differentiation 
of the Vs, Since every 6 is a function of the corresponding 
r, and of that alone, 

db ^db dr _db X — x^ „ 
dx dr dx^ dr r *' 

where r^ = (a; - xy + (y - y'f -\-{z- z'f. 
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mi P ^o dDdb ^^ dD dh dr 

iheretore Sia -sr-j- = SSa-yr- -r -y- 

db ax do dr dx 

_ dD dbi2 X1—X2 dD db^^ x^ — x^ ^ 
^db^^dn^ ^2 ^dbisdr^s r^s 

dD dbi2 x-^ — X2 r. 

(AjOi2 (X^12 '?^12 

because -7-^^ = — r^ &c. 

dxi dx^ 

The above expression can by taking the terms in pairs be 
resolved into a series of terms of the form 

dD db (oL — a') (x — x) 
db dr r 

and is therefore zero on average. 
Similarly 

v-v* / / r^r^/ fx ( db ^db db\ 

can be resolved into a series of terms of the form 

and is therefore also zero on averaofe. 
102. We have then 

V dx m dx doLj 

+ &c. = 0, 
or dividing by A, 

(dlogA h dxdQ\ ^ f ^l^g^. Ai%rfQU&e. = 0. 
\ ^ dx^ mi dxi daj \ ^ dx^ m^ dx^ daj 

And taking mean values of -— from Art. 97, 

(dlogA ,dxD\ f dlogA . dy D \ , 
\ axi dxi niiDiJ \ dx2 dx2 m^D^J 



VIII.] GENERAL THEORY OF THE STATIONARY MOTION. 105 

And since — ^:r— = — ^^ = &c. = ii, 

V ax dxi J \ dx^ dx^ '^ ) 

This is satisfied by assuming, if the molecules be all of 
the same kind, 

J.=€-^'^x (37). 

We have then for the law of distribution of coordinates and 
velocities in a field of external force whose potential is y^ 

n 
7 V 2 

"'^ VD 6-Vx 6-*«rf«i ... dZnd0L^...d^n (38). 



a 



Now -^^f-'^-T, 

and 7ir=iSm(a^ + ^ + 70 = ^, 

the system consisting of n molecules of the same kind. 
But Q = ttittg^ + ^laaitta + &c. 

_ aiD„ + 6]2 Ji2 + ^2^2 + &c . 

Therefore Q = fjinT (39). 

103. With regard to the inclusion of the inter molecular 
forces under the potential % the same considerations apply as 
in Chap. iii. I shall assume that we may use the form 

As an example let us now again consider the system of 
equal elastic spheres in vertical column, Arts. 29 and 45. In 
that case % is the potential of the constant vertical force. Let 
us suppose Ky and therefore all the b coefficients, very small. If 
the X velocity of any molecule, as Wi, be given, the mean value 
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of U2 is — y~ Uj, where B is the determinant 



i) = 



612 613... 
2 623... 



for in this case the a coefficients may each be replaced by unity. 

Now i)i2 contains the term —612 multiplied by axial coefficients 
only, and other terms each of which contains the product of at 
least two 6's. Therefore if each b be small enough, the mean 
value of U2 is ultimately — ^b^^Ui. 

Similarly itg = — 1 b^^u^ , &c.. 



and therefore -^ = /xmiWj = (2 — ^X¥) 11^ , 



dQ 
dui 

or m,fjL = (2-^t¥). 

It thus appears that e"^/^^ differs from e~^^^, when k is small, 
only by small quantities of the second order. Now the factor 
g-(«-Ko)^ or e~*, in the expression for the density in Art. 29 
corresponds, not to e~^'^^, but to e"^'*'^. That is, if yfr be the 
potential of the intermolecular forces, e"" is what €~^^'^ or e~'^'^^ 
becomes in the limit when the forces become infinite for the 
distance r = c. 

And K being small the solution agrees with that of Art. 29. 
It does not however agree with that of Art. 29 when ^-, &c. 
cannot be neglected. Which result then are we to accept 
in this case ? Evidently the result e"'^'*^, because Art. 29 is 
based on the tacit assumption of condition A, and therefore 
holds only for small values of k. To sum up the discussions of 
this problem. 

8 

If K be neglected, the law of density is €~^^^, with ^ = xi • 

If k\ &c., but not K, be neglected, the law is 6~^^^6~'', with 
T— — 

If K^j &c., are not to be neglected, the law is e'^'^^e"", with 

4/xA' 
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Boltzmann's Minimum function. 

104. The following theorem is not to be compared in 
elegance with the theorem founded on condition A, Art. 38, 
but has the advantage that condition A is not assumed. 

Let /{utl . . . UnOOi . . . Zn), or /, be any positive function of the 
coordinates Xi,..Zn, and velocities Vr^..,Uny of a system of 
Qiolecules, and of certain parameters a, 6, Szc. Let 9/ be a 
small variation oi f, a, b, &c., being constant. Let us write 

Let the variation 9/, or qf, be made subject to the following 
conditions, viz. 

JL.,fl-\-qdxi,., dzndui , . , dun 

= 11" 'fdxi . . . dzndui . . . dllny 
or writing dsda- for the continued product of differentials, 

[j...fqdsda- = (1). 

Also /J. . ./TTq Gdsda = jj, . .fGdsd<T, 

or jL.fqCdsdcr^O (2), 

where G is any function the constancy of which is prescribed. 

Again, let 

H = jj...fi\ogf-l)dsd<T, 

H + dH = jj...fTTq {log(/lT?) - 1} dsda. 
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We have then 

H+dH = jj,..f{logf-l)dsda 

-f- 1 1 . . ./log 1 + q dsda 

+ jj,,,fqlogfdsd(T 

+ I / . . .fq log i + g dsdo- 

— I \...fqdsda 



= H + jj,.J^dsda- 

+ jj"'f{t-^i)dsda 

-\-^y..fq\ogfdsda 
= H-\-l{i..fqHsda 

^ jj..,fqlogfdsda, 

because q^ is to be neglected. 

In order that dH may be necessarily positive, or H may be 
least possible, subject to (1) and (2), it is necessary, and it 

is sufficient that jj fq log f dsdo- = 0, that is by (2) that log / 

shall be a function of C, that is that/= J.e"'^^, where 

jj..,fGdsda 
is constant. 

n 

We may write ( - As/D for A and Q for G, The actual 
value of H when minimum is 

|log^ + log^+|logD-| (40). 
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Physical effect of the b coefficients. 

105. From the law of distribution at which we have now 
arrived we can find the mean value of any of the products of 
the form aa', &c. For instance, a^^^D^^jhD, D being the 
determinant of the coefficients. And all the 6's being negative, 
A2 is necessarily positive. If we have only two variables x 
and 2/, and if 

Q — aoc^-\- hxy -\- cy^, 
h being negative, we find A2 = - 6, which is positive. 
So with three variables, x, y, z, and 

Q = aiX^ + tta^^ + a^z^ 
+ &i2^y + bisooz + b2sy2!, 
we find i)i2 = - 2aJ)j2 + 613623 , 

which, every b being negative, is positive. (See Appendix (6) 
as to the sign of these minors.) 

The conclusion to be drawn from this is that if two mole- 
cules are so near each other that b is not negligible, they are 
on average moving in the same direction. There is at every 
point on average a determinate energy of stream motion. It 
follows further from this uniformity of direction that if a0, 07, 
^7, be the mean values of those functions for all molecules 
within any small finite space, (a/Sy, &c. have greater value 
than they would have if the molecules were material points. 

Another consequence of the b coefficients. 

106. But we might also, treating the velocities as con- 
stants, proceed as follows. 

Let for a moment the chance that two molecules distant r 
from each other shall have velocities in x, u ,., ii + du and 
u' ,,,u' -\- du' be proportional to 

in which b may be g — , and w > 3. 
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It is convenient here to treat b as positive, and prefix to it 
the minus sign in the index. That being the case, if b be not 
negligible, that is if r be small, uu is more likely to be positive 
than negative. If r be great, and therefore b negligible, it is 
indifferent. It follows by the theory of inverse probabilities 
that, it being given that uu' is positive, r is more likely to be 
small than if nothing were known of the sign of uu\ If on the 
other hand it be given that uu' is negative, r is less likely to be 
small. 

Suppose a molecule at P with velocity u, and an element of 
volume dco at P', distant r from P. Let the chance that there 
shall be within do) a molecule moving with velocity u' be 'yfrdco, 
when PF\ or r, is great, and 6 = 0. Then when b=\^0, the 
chance becomes 

which is greater or less than ylrdco, according as uu' is positive 
or negative. 

Further, the chance that there shall be a molecule within 
dco with velocity either it or — u is, since a priori -f u' and — u 
are equally probable, and the system is supposed infinite, 

uu' being now treated as essentially positive. That expression 
is equal to y^day if 6^ and higher powers of b are neglected. But 
if 6^, &c. be taken into account, it is greater than -y^dw, 

107. It follows from this result that under the law of 
distribution of coordinates and velocities at which we have 
arrived in this chapter, if at any instant we calculated for each 
of the molecules the value of S6, or the sum of the b coefficients 
connecting the u velocities of that molecule with the velocities 
of each of the others, the average value of S6 would be 
greater than it would be were all the molecules scattered 
through space at random. That signifies that the molecules of 
the entire system would be collected into denser and rarer 
masses — the effect of which evidently is to increase the mean 
value of %b for all molecules. Further the ratio ^'^ : p would 
have a greater value for the denser masses than for the rarer 
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ones, or the stream motion would be greater for the denser 
masses than for the rarer ones. So that the system would tend 
more and more, with increasing number of molecules in a given 
space, to assume the form of a number of denser aggregates, say 
clouds, moving through a comparatively rare medium. 

108. It is not essential to my argument, but I will here 
give what appears to me the best possible definition of density 
at a given point at a given instant, for molecules of finite 
dimensions. 

I would define namely as follows. The density at any 
point P at any instant is — S6, where 26 is the sum of the 
h coefficients, all being negative, of all the molecules of the 
system referred to the point P — that is, h for any molecule is 
a function of the distance of that molecule from P. 

109. It is assumed throughout this investigation that the 
quadratic function Q is positive, that is, that the coefficients 
a, b, are such that Q cannot be made negative, whatever values 
we assign to the velocities. It is assumed that the distribution 
of the molecules in space, and therefore the coefficients 6, are 
on the average such as would be calculated by the method above 
indicated, the whole number of molecules of the system in a 
given finite space, i.e. p, being known. Should that condition 
ever fail, the motion of the system can no longer be represented 
by the law considered in this chapter. It is reasonable to 
expect that the change of the mathematical formute will 
coincide with a change in the physical state of the system. 

110. The condition that Q shall be necessarily positive, 
whatever values be assigned to the velocities, is that the deter- 
minant D, of the coefficients in Q and all its coaxial minors must 
be positive. See Appendix {p). As the h coefficients increase in 
absolute magnitude, Z), being positive, diminishes. But as the 
density p increases cceteris paribuSy the b coefficients do generally 
increase in absolute magnitude, and D, with all its coaxial 
minors, generally diminishes. The b coefficients are generally 
also functions of the mean kinetic energy, or temperature, 
as well as of the density, and if the temperature be high 
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enough, probably diminish as that increases (Art. 90). As the 
density increases, or temperature diminishes, or both, the 
system may ultimately arrive at a state, in which D or some of 
its coaxial minors becomes zero or negative, and then the change 
of physical state ensues, whatever its true description be. Now 
to any group of molecules we choose to select, and therefore to 
all the molecules forming one of the denser masses or clouds 
above mentioned, belongs a particular coaxial minor, that 
namely obtained by striking out all the axial constituents of D 
except those which belong to molecules of the group. Gener- 
ally therefore the condition may be expected to fail for some of 
these coaxial minors before it fails for the complete determinant 
D. So that as density and temperature vary continuously, the 
change of state, or liquefaction if such it be, will take place 
partially and gradually, some portions of the system having 
passed into the changed state, while other portions remain in 
the original state. Further, the higher the temperature, at all 
events after a certain point is reached, the greater the conden- 
sation necessary to make D, or any coaxial minor vanish, and 
therefore to produce the physical change in question. 

Such is the process which our analysis leads us to expect. 
Physicists may consider how far it corresponds with what is 
known to take place in gases under condensation, or on what 
(if any) farther hypothesis it may be made to correspond 
with it. 



CHAPTER IX. 



OF MOLECULES AS CARRIERS. 



111. The molecules being elastic spheres, and condition A 
being assumed, to find the mean free path for a sphere issuing 
from collision with velocity w. As the mean free path has 
received more than one definition, I define it as follows. Sup- 
pose a sphere m to start from a collision with velocity g), and to 
undergo n successive collisions. After each collision let each of 
the colliding molecules have restored to it the velocity which it 
had before collision. Then our sphere will describe n paths 
lj,l2... In between its collisions, each with the same velocity «. 
The mean firee path for velocity co is then 



^w — 



ii + 4+...+^n 



n 

when n becomes indefinitely great. 

If yfr be the velocity before collision of a sphere mf which 
collides with m, E the angle between the directions of motion 




of the two spheres before collision as in the figure, the relative 
velocity of approach is 

R = Vft)^ + -^/r^ - 2(oylr cos E, 

Let c be the sum of the radii of m and mf. Let/('^) dyfr be 
the number of m' spheres in unit of volume. Then the number 
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of collisions which m undergoes with the m' spheres per unit of 
time is 

h I I /('^) dyfrTTC^R sin EdE^ird'N^ suppose, 

and the mean free path for an m moving with velocity «, if no 

collisions occur except with spheres of the m! class, is \„ = ^rp- . 

If there be two classes m^ and m^, and m can collide with another 
sphere of either class, we shall have a corresponding number of 
collisions per unit of time between m and m^, which we will 
denote by irc^N^' Then the mean free path for velocity co is 

\^ = — —r^ ^TrT and so on. To calculate the mean free path 

TTC^ (iVi + iVa) 

is to calculate Ni + i\^2) ^c- 

112. Another method of obtaining this result is as follows. 

Let S be the whole space in which our spheres are moving. 
Consider a sphere m moving with velocity o>, and another sphere 
mi whose velocity and direction of motion are comprised re- 
spectively within the limits i/r . . . y^ -\- dyfr and E ... iJ -\- dE^ as 
in the figure. 

Now TUi may be anywhere within S. If its centre be at this 
instant within the cylinder whose base is a circular area of centre 
m and radius c, the sum of the radii of in and mj, at right 
angles to iJ, and whose height is Rdty then, but not otherwise, 
a collision will occur between m and nij within the time dt after 
this instant. 

The chance that such collision shall not occur is then 

TTC^Rdt 



1- 



S 



Now the number per unit of volume of spheres tUi whose 
velocity and direction of motion are comprised within the limits 
aforesaid i^ f(yfr)d'\lr^ sin EdE. And therefore the number 
within S is >S/('>|r) dslr ^ sin EdE. Let us now assume that the 
aggregate volume of all the vii spheres within S is negligible 
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compared with 8. Then tjie chance that during dt m shall not 
collide with any one of the class mi is 

raised to the power of 

^/W^f Jsin^d^, 

that is, in the limit when dt becomes infinitely small, 

1 - 7rc^Edtf(yfr) df ^ sin EdL\ 

and the chance that m shall not during dt collide with any mi, 
whatever be the values of yfr or E, is by the same reasoning 

1 - ircHt I / fi^fr) di/r J iJ sin EdE = 1 - irc^dtN^ suppose. 

Similarly if there be other classes of spheres, mg, m^ &c., 
with which m may collide, the chance that it shall undergo no 
collision in the time dt is 

1 - ircHt {Ni + iVg 4- &c.) = 1 - ird'dtN, 

By the same reasoning the chance that m shall undergo no 
collision in n successive intervals of time dt is 

1 - TTC^Nndt = 1 - TTC^m, 
if t = ndt 

Also if X be the space described in time t 

\ 

\ = ct)^, or t = - . 

(O 

Therefore the chance that m shall survive for a distance X 
without undergoing any collision is 1 — ttcW - . 

Let l-7rc^N- = 6 (X) = <f). 

CO f ^ ^ r 

Then <f> is the proportion of the whole number of spheres m 
which, starting with velocity ©, traverse the space \ without 
collision. Also <^ (\ + d\) = (f>(\)(j> (d\). It follows that 

dd> ^ m^N , J .. 7 'rrc^N 
-r = _ A = _ ]cA if ]c = , 

aX ^ ft) ^ o) 

8—2 
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and therefore <^ = Ce'^K Also </> = 1 if X = 0. Therefore C = 1 
and = e~^^. Whence if \ be the mean value of X given co 



Jo 



e-ic^XdX=l= "" 



which is the same result as we obtained above (HI). 
The expression 

TTC^ I dyfr I \in EdEf{yfr) d'>\r Vo^ + i/r^ - 20)-^ COS E 
Jo Jo 

can be calculated numerically in terms of ttc^ and h if 

TAe fundamental equation. 

113. Let (? be any quality, e.g. mass, colour, charge of 
electricity, momentum in given direction, &c. which a molecule 
may possess in greater or less degree, and which it can carry 
with it unchanged by its own free motion, and which, as 
regards the aggregate possessed by two colliding molecules 
is not changed by collision. Suppose a quantity of gas con- 
tained in a vertical cylinder of which the height is very small 
compared with the diameter of the base. Every molecule that 
enters the cylinder through its upper flat surface shall enter 
charged with the quantity Gi on average of our supposed quality 
G. Every molecule that enters the cylinder through the base 
shall enter charged on average with the quantity Gq. When 
the motion has become stationary, every molecule within the 
cylinder at a height z above the base will be on average charged 
with a quantity of the thing in question between Gq and (?i. 
Let it be called G(f). G(z) in stationary motion will be a 
function of z. 

114. Now consider a layer Pi of molecules in the cylinder 
between the planes Zi and z^ + dzi. If the molecules within 
that layer were enclosed between two elastic planes Zi and 
Zi -\- dzi for a finite time, the distribution of velocities among 
the molecules within it would become, under the circumstances 
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supposed, completely independent of the quantity of Cr which 
they happened to carry. It would become in the strictest sense 
of the term " molecular ungeordnet." In other words condition 
A would prevail. 

Collisions would occur within that layer, and any molecule 
emerging from such collision would be as likely to be moving in 
any one direction as in any other. The chance of its having 
any particular direction would be independent of its charge G. 

Consequently if any such molecule, having undergone col- 
lision within the layer Pj arrives without further collision at 
the layer P, the mean quantity of G which it brings with it is 
the mean quantity of G per molecule for all the molecules 
within the layer at Pj. This is the fundamental assumption. 

115. Now let Gz be the mean value of that quantity for 
all the molecules within the layer z ,,, z -\- dz Bit P, Suppose a 
molecule to arrive at that layer from above with velocity 
between v and v + dv, and in direction making angle between 
6 and 6 + d6 with the vertical. Let X be the distance which 
on its arrival at P it has traversed since its last collision. Then 
such last collision must have occurred at a height A, cos above 
P. Therefore the quantity of G with which the molecule is 
charged is the mean quantity of G for the height ^ + X cos 0, 
that is 

Giz + Xco80)= Gz -hXcos 0^ , 

because \ cos is small. 

The number of molecules which, having velocity and direction 
between those limits, pass through the horizontal plane P per 
unit of area and time is vf(v) J sin cos 0d0dv, 

Hence the quantity of G carried through the layer P per 
unit of area and time by molecules coming from above is 



r2 (-2 fin. 

vf(y) dv\ sin cos 0G {z)d0 + \\ \ cos^ sin 0d0 ^ f{v) vdv 
Jo Jo cL^ 
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if all values c 



and taking the mean of all values of X and v it is 

dG 

dz ' 



Here \—\ f(v) vXdv, 



Jo 



By the same reasoning the quantity of G carried through 
the layer P per unit of area and time by molecules coming 

from below hlG{z) — ^\~^ . 

And the quantity passing through P downwards exceeds 
the quantity passing through P upwards per unit of area and 

time by ^\ -T- . 

116. This is the fundamental equation from which Boltz- 
mann calculates the rate of diffusion, viscosity &c. in gases as- 
suming condition A to prevail. It would represent with complete 
accuracy the initial motion of the system, if the quality G were 
suddenly attributed to the molecules in the degree proper to z, 
but without regard to their velocities. 

It takes no account, as he points out, of any variation in the 
distribution of velocities which may be caused by the quality 
G itself, for if for instance G be momentum in given direction 
at right angles to the axis of our cylinder, such a change in 
distribution must exist. It is assumed that the effect of this 
change in distribution on the mean free path X is so small as to 
be negligible, a condition which must be secured if the given 
terminal conditions, e.g. Gq and (?i, do not differ very widely 
per unit of distance. 

117. The method also tr.kes no account of the stream. The 
effect is to set up a stream, so to speak, of G from the upper to 
the lower parts of the cylinder. And that, as will be found, 
causes a small error in the result. The error may indeed be 
negligible so long as the stream velocity is very small compared 
with molecular velocities, which will be the case if the terminal 
conditions do not differ widely per unit of distance. 



IX.] OF MOLECULES AS CARRIERS. 119 

It therefore does not sensibly affect the results calculated 
from Boltzmann's formula in any case to which as I understand 
he proposes to apply that formula. It is worth while however 
to consider the sign and order of magnitude of this deviation 
from complete accuracy. 

We are asked to suppose that in stationary motion molecules 
carrying positive 0, have a mean velocity u downwards ; mole- 
cules carrpng negative G have a mean velocity u\ u and u 
being functions of z, Hydrodynamical considerations make Gu 
and Gu^ sensibly constant throughout the cylinder in stationary 
motion. 

But when collisions take place in the layer Pi, condition A 
is supposed to prevail so that the molecules issuing from 
collision are as likely to be moving in any one direction as in 
any other, and to carry with them, in whatever direction they 
do move, the quantity of G due to the height 0. 

In fact molecules with positive (?, having before collision 
the mean velocity u downwards, will issue from collision with 
some mean velocity downwards, less indeed than u but com- 
parable with it, and therefore will not be moving in all directions 
indifferently, but on the average downwards. And molecules 
moving downwards will possess on average a greater quantity 
of G than that due to the layer in which their last collision 
took place. For let o) be the absolute velocity of a molecule 
m before collision. Let ^jr be the velocity before collision 
of the other molecule m concerned, E the angle between 
their directions as in the figure (Art. 111). Then the velocity of 
m after collision is the resultant of (1) the common velocity of 
m and m', (2) their relative velocity, in whatever direction it be, 
after collision. But for the relative velocity after collision all 
directions are equally probable. The relative velocity therefore 
contributes on average nothing to the resultant. 

The resultant velocity after collision is then on average 
the common velocity. And the velocity after collision re- 
solved in its direction before collision is the common velo- 
city resolved in direction «, that is if m = m' a> + i|r cos E. 
We have then for the mean value of the velocity after 
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collision resolved in the direction before collision the ex- 
pression 

A I dyjre-^'f'' -^^ f Vw^ -f- i/r^ - 2G)i/r cos E (« 4- -v/r cos E) sin EdE. 

Jo Jo 

That can be calculated as a function of co. As however we 
intend to treat it as negligible, it is not worth while to calculate 
it here. But we can easily prove that it is necessarily positive. 
We may therefore conclude that the mean velocity after col- 
lision of a molecule carrying positive G resolved in its direction 
before collision is/o) on average, where/is some positive quantity 
less than unity. 

It follows that the mean velocity downwards of the mole- 
cules carrying positive (?, which before collision was u, is after 
collision fu. 

We have thus proved that the molecules carrying positive 
G issue from collision in the layer Pj with mean velocity fu 
down the cylinder. For the same reason the molecules carrying 
negative G issue from collision with mean velocity /i^o upwards. 
Therefore of the molecules which reach the layer P, having had 
their last collision in the layer Pi, the proportion which carry 
positive G is not the same as the proportion of all the molecules 
within the layer Pi which carry positive G, as it was assumed 
to be in forming the fundamental equation, but is greater. 

118. The effect of the introduction of the h coefficients in 
this case. 

It has been shown in Chapters iv. — vi. that molecules near to 
one another are on average moving in the same direction. The 
energy of their relative motion is therefore less than it would 
be were condition A assumed, and therefore their encounters 
pro tanto less frequent. That tends to increase the mean free 
path. 

On the other hand under assumption B the molecules tend 
to become aggregated into denser and rarer masses. That 
tends to increase the frequency of their encounters and there- 
fore to diminish the mean free path. 
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It is probable that for ordinary gases under ordinary con- 
ditions the b coefficients will be very small. And if so the 
effect of the b coefiScients on the phenomena of diffusion &c. 
will be of the same order of magnitude as those sources of error 
mentioned in Arts. 116, 117, which as we have said Boltzmann 
and others have agreed to neglect in calculations concerning 
these phenomena. I shall not therefore here further consider 
the effect of the use of assumption B instead of A on such 
calculations. 



CHAPTER X. 

ON THE HYPOTHESIS THAT HEAT OR TEMPERATURE IS REPRE- 
SENTED BY THE KINETIC ENERGY OF MOLECULAR MOTION. 

Natanson's Theorem, 

119. M. Ladislas Natanson (Interpretation cinetique de la 
fonction de dissipation — Bulletin de VAcademie des Sciences de 
Cracovie, Deeembre 1893) considers a mediutn composed of 
molecules in motion. The components of the molecular velocity 
of a molecule are u, Vy w, those of the apparent or stream 
velocity of an element of volume are ^, rj, f. I have here inter- 
changed Natanson's symbols to make them agree as far as 
possible with my own, and so avoid confusion. Evidently, taking 
mean values, 

iiz=v = w = 0. 

Let if = i jjj pi^' + v'-^ ?0 dxdydz, 

p denoting density, and 

E=^ iij piu' + v' + w^) dxdydz. 

The integrations are throughout all space occupied by the 
system. 

120. M. Natanson gives no definition of ^, ?;, f, but is, as I 
understand, willing to accept Boltzmann's definition — that f at 
P is the momentum in x of the centre of inertia of all the 
molecules contained in a sphere of radius r described about P 
as centre divided by the volume of that sphere, in the limit 
when r becomes infinitely small. 

I should propose myself to define it as in Art. 61 by the 
function /, so that f, 77, f, Natanson's components of stream 
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velocity, agree nearly with the ^, 97, f, of Chapter v. But 
M. Natanson does not on his own assumption require any 
further definition. 

The normal pressures per unit of surface are 

Pxx = pU^ Pyy = pV\ Pzz = py^' 

The tangential pressures are 

P^y^pUV, PxZ^PUWy Pyz = pVW, 

121. He then gives the following fundamental equation, 
viz. q being any function of -z^ + ^, v + rj, w-\- ^, the time varia- 
tion of q is found from 






Here X, F, Z are the components of external force, and —- 



dq 



dq 



dt 
the change of 5' due to encounters, ^ the total change of q 

with the time. 



Writing u + ^ for q, he finds 

And again, putting q = (u + ^y + {v + r^f + (w + ^y, he obtains, 
neglecting terms of the third order, 

P J^{^^ + V^ + ^ + u^ + ^ + w') 
+ ^ (2fK + 2^p^ + 2r,o^) 
+ ^ (2^puv + 2rjpiP + 2^pvw) 
-\--r{2^puw + 2r]pvw + 2^pvf) 



9 



+ Ps?(^' + ''' + ^' + ^' + ^' + ^') (5)- 
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Combining (4) and (5) we get 

= P^f(^' + v' + ^' + u'' + v' + w') (6). 

Let f = a, ^/ = b, §=c, 

dx ' dy ' dz 

dz dy ' dx^ dz~ ' dy^dx~ • 
Add to the first member of (6) 

i (m^ + w^ + «i;^) J + |p (m2 + ,r^ + w2) (a + 6 + c), 

which is zero. Also write 

d 
dt' 



-T y (u^ + v^ + w') 



And then in the equation so obtained integrate for x, y, z, 
and make l^ + mrj 4- ?zf = at the bounding surface, I, m, n 
being direction cosines of the normal to that surface. Assume 
further that X = F = ^ = 0, or no external forces act. 

The result is 
,— + / 1 Upu^a -h pv'^b + pvfc + pvwA + puwB + puvG) dxdydz = 0. 

Treating equation (5) in the same way, we obtain 
-7- — ill {pv?a + pv"^!) + pw'^c + pvwA + puwB + puvG) dxdydz = 0. 

And therefore, there being no external forces, -7- = ~' ~Ji > 

or the kinetic energy of stream motion can increase or diminish 
only at the expense of molecular energy. 
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122. Let us now write 

Sp = pi? + pv^ + pw\ 

F = (p - pu") a + {p- pv^) b + {p- puP) c 

— pvwA — puwB — puvG, 

By a further use of the fundamental equation, putting 
q = (u-^ 1)2, and then q^v+rj x w + ^, we find 

+ 2 (pvwf + 2 (pmuf + 2 (pwi;)'}. 
And finally 

dE dK 



dt 



= -~^\\\{F-pe)dxdydz, 



where ^ = a + 6 + c = -^ + -^ + :/. 

aa? a^/ dz 

123. It appears then that E, the molecular or heat energy, 
gains or loses at the expense of the energy of stream motion, 
according to the law of encounters " loi des chocs moleculaires." 
If, that is, as the result of encounters the tangential pressures 
are increasing, E diminishes, and K increases. Now this is 
precisely what happens according to Chap. V. ante, when, 
for instance, instead of being material points, the molecules 
become elastic spheres of finite diameter. As the result of 
collisions in that case ^^ in my notation becomes ^^ + f '^ &c. 
And as we have seen. Art. 105, (a7)2,(a/3)^ and (^87)^ increase, 
but only up to a certain point, at which namely the ratio f 'Vf ^ 
has the proper value for stationary motion. The relation 
between the energy of visible motion K and that of molecular 
motion E is in Natanson's theorem the same as that between 
K and heat in Thermodynamics. 

Messrs Bryan and Boltzmann's method, 

124. In order to confirm the hypothesis that the temperature 
of a gas or other substance is proportional to the mean kinetic 
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energy of its molecules, Messrs Bryan and Boltzmann (Vienna 
Sitzungsherichte math, naturiu. Classe, Band III., Abtheil II. 
Dec. 1894, also referred to by Boltzmann, Gastheorie, p. 136, et 
seq.) treat the following case. X and Y are two infinite parallel 
planes. In the space Si to the left of JT is a gas A, In the 
space ^2 to the right of Y another gas B. In the space S 
between X and Y there acts on molecules of gas A a force 
directed from Y whose potential becomes infinite on Y, and 
there acts on molecules of gas B a force directed from X whose 
potential becomes infinite on X. These conditions insure that 
no molecule of gas A will be found in S2, and no molecule of gas 
B will be found in Si. But in the space S there will be a 
mixture of A and B, and encounters will take place between 
their respective molecules. 

Messrs Bryan and Boltzmann use generalised coordinates, 
but it wdll be sufficient in this notice to deal only with velocities 
of translation, as though both A and B molecules were elastic 
spheres. 

It is assumed that the mean kinetic energy of A molecules 
at any point within S is the same as their mean kinetic energy 
within Si, and that the mean kinetic energy of B molecules 
within S is the same as their mean kinetic energy within S2. 

It is assumed further that as the result of encounters taking 
place within S, between A and B, the molecules of A acquire, if 
they have it not already, the same mean kinetic energy as those 
of B. That result seems to me to be questionable or to be only 
approximately true, if the molecules have finite diameters or 
spheres of action. But I state this with diffidence when such 
authorities are against me. 

Making however these two assumptions, the conclusion 
follows that the mean kinetic energy for the A molecules in Si 
must become equal to that of the jB molecules in fi^g- ^^^ 
therefore that if the mean kinetic energy be originally different 
for A and for B, a finite quantity of energy will pass from Si to 
S2 , or vice versa. 

The condition that no energy shall pass is that the mean 
kinetic energy shall be the same for A as for B. But we may 
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regard the two gases as two bodies in contact, which condition 
is more nearly realised as S is diminished indefinitely. The 
known experimental condition that no energy shall pass between 
two bodies in contact is that both shall be at the same tempera- 
ture. Hence in the case of the two gases at all events kinetic 
energy plays the same part as temperature. 

Professor J, J, Thomsons method. 

125. The same hypothesis is tested in a different way by 
Professor J. J. Thomson (Application of Dynamics to Physics and 
Chemistry, p. 91). He puts the following case. Suppose a tube in 
which is moveable a piston of very great mass M. On one side 
of it is a set A of material particles each of mass mi in motion. 
On the other side a set B of material particles each of mass mg in 
motion. Considering the two systems of particles as two sub- 
stances, and the piston as a conductor of heat, the condition 
that no heat shall pass through the piston from one substance 
to the other is that the two substances shall be at the same 
temperature. If temperature be represented by kinetic energy 
of translation, then the condition that neither set of particles 
shall on average either lose or gain energy by collisions with 
the piston ought to be that the mean kinetic energy of transla- 
tion is the same for the A as for the B set. 

Let U be the velocity of the piston, v^ the velocity 
normal to the surface of the piston of a particle m^ before 
impact. It is assumed that particles striking the piston re- 
bound from it as perfectly elastic bodies. Let U\ w/ denote 
the values which U, u^ assume after impact. Then we have 

m^Ui -f- MU = miV + MU\ 
mi^2 ^ jf [72 ^ ^^^'2 ^ jjf ;7'2 

From which deduce 

, „ ,, 4ilf2mi CT^ - 4ifmi (if - m^ Un^ - ^Mm^^u^^ , . 
m, {u, - V) = TW^^^ <^)- 

Now add together the equations of this form for all collisions 
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which take place in unit of time. Let Ni be the number of 
such collisions. Also let 6 = \MU'^ be the kinetic energy of the 
piston, and 6i the mean kinetic energy of the particles of the 
set A which strike the piston. As I understand Professor 
Thomson, the mean is taken for all collisions, not for all 
particles. 

It is assumed that the velocity U of the piston is zero 
on average of time, because the piston is being struck on the 
other side by particles of the B set and the pressure is supposed 
to be the same on either side. Hence we may neglect in the 
expression {a) the term which contains U in the first degree. 

The expression then takes the form 

N.^e. = N.i2e-m^-M-^. (1). 

Note. This should, I think, be f^j instead of J^j because as M 11^—26-^ 
so ^1^1^= 1^1, ^1 denoting the whole mean kinetic energy of all molecules 
which collide with the piston or with other portions of the elastic surface 
bounding the A molecules, i.e. 

But this does not affect Professor Thomson's argument. 



Similarly for the B set we have 

' (ilf+m,)* 
and since dO = — N^dOi — N^d2 = on average 



^-^#, = Jf,(2«-l«,)-j^-=, (2), 






{M + m^f ' (if + nhT ■ 



This gives id -16,^^1 {6, - 6,) 






^e-i6,=---Aie.-e.). 
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Substituting in (1) and (2) 

and therefore d0i = dO^ = 0, if 0^ = 02. 

It follows that if 02 > 0i the effect of all the collisions is 
to increase 0^ at the expense of 02, and vice versa. And the 
condition that no energy shall pass from A to B or from B to A 
is that 02 = 01, or the mean kinetic energy of translation for all 
particles colliding with the piston is the same for A as for J?. 

126. Now let f(u)du be the number per unit of volume of 
particles of the A set whose velocities normal to the piston lie 
between u and u-\-dii. Then the number of such particles 

which strike the piston per unit of area and time is I f{u) udu. 

Jo 
The sum of the kinetic energies of their motion normal to the 

/•OO 

piston is ^nii I f(u)u^du. The mean kinetic energy of motion 

Jo 

normal to the piston for all A particles irrespective of their 

/•QO 

striking the piston is ^nii I f(u) v?du. Let f{u) be so chosen 

Jo 
that the ratio 

/»00 /«00 

I f{u)ii^du I f(u)u^du 

Jo_ . J O 

/•OO • /.QO i 

\ f{u)udu \ f(u)du 

Jo Jo 

is an absolute constant, and therefore the same for A as for B, 
That being the case, when 0^ = 0^^ or the mean kinetic energy 
per collision is the same for A as for B, the mean kinetic energy 
per particle is also the same for A as for B. And therefore the 
condition that on the whole there shall be no transfer of energy 
through the piston from ^ to 5 or vice versa is that the mean 
kinetic energy per particle is the same for A as for B. 

B. 9 
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127. If we assume for A particles 

and for B particles fi^^) — CV"^^^^*^^, 

we satisfy the condition. 

Possibly there may be other forms of f(u) which satisfy it. 
But however this may be, some assumption with regard to the 
relative frequency of different values ofu appears to be necessary. 

Further, f{'u) being properly chosen, we satisfy at the same 
time the two conditions (1) that the pressure shall be the same 
on either side of the piston, the number of particles per unit of 
volume being the same for A as for B, and (2) that there shall 
be no transfer of energy. 

The law of distribution of velocities thus appears to be 
of the essence of the matter, if the molecular motion is to 
represent heat. 



The Second Laiv of Thermodynamics. 

128. Our gas system being in stationary motion according 
to the laws above investigated, we might impart to it a certain 
quantity of heat hE. That being done, it would assume a new 
form of stationary motion, in which the mean kinetic energy T, 
the volume «, and the mean potential energy x^ will be gener- 
ally different from what they were in the original motion. If 
P be the external pressure which in the original motion was 
just sufficient to prevent expansion, the change of volume dw 
can be effected only by doing an amount of work P^ft) against 
the external force P. A part of the energy dE supplied to the 
system is spent in doing this work. The remainder is spent in 
increasing either (1) the mean kinetic energy T, or (2) the 
mean potential energy, ^, of the intermolecular forces, so that 
dE = XdT + 9% + Pdcoy X being the number of molecules. 

129. If the external forces P also have a potential, there is 
in general a determinate relation between dco and dT, and a 
determinate relation between do^ and d^- So that in this case 
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there is only one independent variable, and no question can 

dE . 
arise whether or not 7^- is a complete differential. Any theorem 

proved on the hypothesis of P having a potential is not the 
Second law of Thermodynamics as usually understood. 

But the external force P is in Professor Thomson's language 
a "controllable" force. It is possible therefore by suitably 
varying P to maintain the system in stationary motion with 
any arbitrarily assigned values of T and w. But dT and 
9ei) being given, 9% is determinate. There are then, P being 
controllable, two independent variables. It is convenient to 
use P and m for independent variables. 

130. If in rectangular coordinates the abscissa represents 
ft), and the ordinate represents P, any point in the diagram 
represents a determinate state of the system, in which it is in 
stationary motion with the values of P and w corresponding to 
the point, it being understood that the controllable force P is to 
have such value given to it as will be necessary to maintain 
that motion stationary. So also any curve drawn on the plane 
of the diagram represents a series of states through which the 
system might be made to pass successively, P being made to 
vary as required. But it is understood that the change of 
state, that is the passage from one point in the curve to 
another, is effected so slowly that stationary motion is always 
attained, and therefore the kinetic energy corresponding to the 
controllable coordinate co is for our purpose to be taken as zero. 

If the curve be a closed curve, the system having passed 
through all the states denoted by the curve, has at the end 
of the process the same values of P and o), and is therefore 
for all experimental purposes in the same condition as it was 
at the beginning. Contemplating the system itself alone, 
nothing that we can observe by experiment has happened to it. 
Now if 6 denote the absolute temperature, the second law of 

Thermodynamics asserts that in this complete cycle 1 — = 0, 

dE . 
or -^ is a complete differential of some function of 6 and co, 

9—2 
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That function was called by some writers the Thermodynamic 
function, by others the Entropy of the system. 

dE 

131. As I have said, -^ being a complete differential, no 

change takes place in the system itself as the result of any 
complete cycle of operations of the kind described. But wE^^O, 

and 1P8g)=|=0. Generally a quantity of heat denoted by \dE, 

drawn from some external source, has been converted into 

mechanical work \Pd(o outside of the system, and a further 

quantity has been transferred from a hotter to a colder external 

body, as shown in treatises on Thermodynamics. But if our 

system be enclosed in a non-conducting envelope, so that no 

heat or other form of energy is allowed to pass into it from 

without, or out of it into external space, it may by suitably 

varying the controllable P be made to expand or be compressed, 

and in so doing it does external work, positive in case of 

expansion, negative in case of compression. This work is done 

at the expense either of the absolute temperature 6 or of the 

potential ')(^ of the system itself The line traced on the 

diagram in such a process, a line that is for which dE = 0, is 

called an adiabatic line. If the cycle be A BCD on the diagram, 

AB denoting expansion at constant absolute temperature 6, BG 

adiabatic, CD compression at constant absolute temperature 6\ 

and DA adiabatic, it is found that the heat which must be 

supplied from without to effect the expansion AB is to the 

heat withdrawn during the compression OD, both measured in 

mechanical units, as 6 : 6\ If temperature is represented by 

dE , 
the kinetic energy, we ought to find T : T' :: : 0\ or ^^ is 

a complete differential. 

132. In the Philosophical Magazine, January, 1876, a 
proof of this proposition was given by the present writer. 
This was employed afterwards in an improved form by 
Dr Watson in his Kinetic Theory of Gases, and has been 
accepted by Professor Bryan. It was founded, as will be seen, 
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on Boltzmann's theorem, that the chance of a group of \ 
molecules having their coordinates between the limits 

a?i . . . iUi 4- dxi , 



(in which position their potential is x) is proportional to 

e"'^^ dxidyi . . . dzf, . 

This is established only on the assumption of condition A. 
I will here give that proof of the second law in order that we 
may subsequently see what modifications are required in it 
when we abandon condition A, and employ the more general 
method of this work. I here use h instead of 2/i of Art. 44. 



We have 9^' = \aT + 9% + Pd(o. 



3 

Let T = ^= the mean kinetic energy of a molecule. 



Then T ^ "3" ' 

, BE .2A 2A- 2A 
and Y^^S^^-^'^'S^^^'S W* 

Now substitute for Pdco from Clausius' equation, Art. 24, 

f Pw = SJm {u' + v' + w') + ^ttR?, 

or Pdco = lXm(u'-\rv'-{-w')—-{-^Xl^ — . 

"^ CO CO 

Now if there be a general cubical expansion, co becoming 
(o 4- day, every line r in the system becomes 

dco 
r + i-r, 

or 9r = 4 — r. 

^ CO 

Therefore the last equation becomes 

"^ ^ ' CO 
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Again, for each configuration 

^ da) 

Hence taking mean values 

aw 

Here V^ is the mean for all confie^urations of ~ , that is 
dco "^ dco 

the mean of a variation, while dx is the variation of a mean. 

The distinction between 9y and v^ dco is of the essence of the 

^ dco 

matter. 

Again, 

.V- / o o ox 3^ \do) 
^l^m (u'' + v' + W-) — = T — , 
"^ ^ ^ CO h CO 

if X, be the number of molecules of the system. 

We have then 

„^ Xdco dv ^ 
k CO dco 



3 ^ CO ^^ da 

And therefore 



and ^ P9ct) = |X ^r :? 9^- 



y = A-ar-f^.^ + ix-^--3^aco (II). 

133. So far our equation is derived from physical con- 
siderations. 

Now consider the auxiliary function 

u = log I 1 1 e~^^ dx^. . .dzK • 

Then du = ^r 9/i + -7— 9«, 

du , 



and "^ 9A = ~- ^- 9 f [ [ e'^^ x^^'i . . . d^x 



X.] THERMODYNAMICAL RELATIONS. 135 

because, since the chance of the configuration Wi,,. Zk is 

e~^^dxi ... dzxi 
evidently 

X = jjj- ^'^"^ %^^i "^dz^-r jjj^-^'' dx,, . . dz^. 
Again, -j-dco = - ~r~. i If e'^^ h -^ dx^,.. dzj, 

/ 1 1 e-^^ d (dx^ . . . d^A), 



I €-^^dxj..,dzK 

the last term being necessary because every element of length 
dxi ... is increased in the ratio 1 : 1 + J — by the general 
cubical expansion above mentioned. 

That is -r-d(o = — h -j^ day -h - dco, 

do) do) 03 

because 9 (dx^ . . . dz,,) = X — , 

CO 

and du = — ydh — A -^ H , 

"^ dco (o 

or fx — -p^aft) = iaz^ + fYaA (iii). 

ft) (i(0 

Substitute this in (II), and we have 

which is a complete differential. 

134. Let us now consider what modification has to be made 
in the above proof when Q, instead of being the sum of squares 
of the velocities, has the value given to it in this treatise. 

With the form of Q containing the h coefficients, the 
relative frequency of the several configurations is different 
from what it was when Q was a sum of squares only, but the 
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values of ;)^ and -~ are not altered for any given configuration. 

Hence % and -j^ have different values from those which they 

had in the former case, y and ^ are the same as in the former 

case. Also y = ^ - = -y-, . 
Zfjih 2ti 

Nevertheless the physical equation 

dE = XdT-{-dx + Pd(o (I) 

remains unchanged in form with the altered value of ^x, 
Further when we substitute Clausius' expression (Art. 24) for 
Pdco, the equation 

^ = Xc)logT + -ax + |X^- 3^30) (II) 

remains unchanged in form with the altered values of h, x> 

and ^. 
dco 

Now let us consider the auxiliary function ti. We might 
in the proof above given have put it in the form 



,t = log 1 1 1 ( -J e-'^ (,+S(aH^Hy^)) a^^ , , , ^^^^^^ ^^^ 



In our present problem it would then take the form 

A. 

u = log 1 1 1 e-^'^ (-j VSe-'^^ dw^... dz^doi^ ...dry, 



A. 



= log [[ e-'^'x (h' Vi> e-^^dsda, 



where ds stands for dx^...dzx and da for dai...dyx' Then 
h^ and co are the two independent variables. Let 

Then du = ^j-. dh' + -^— dco. 

art dco 
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Now ^,dh' 
ah 



■^h' ff _,./^ d ,^, , dF' 



\k'^'i^m^''^^-i¥^'^ 



jje-'^''cFdsda^^ ^ "^^^ 

= - X^h' + J. ^ lle-'^'m' %dsda (1). 

jje-'^'>cFdsda^^ "^^ 

. . du ^ 
Again -v- 0(0 

jje->''>cFdsda^^ '^"^ 

the last term referring to the variations of all the elements of 
length dx, dy, dz, consequent on the expansion denoted by dco. 

That is ^- d(o 
d(o 

jje-'''^FdsdJ^ ^ ^ '^"'^ 



^ dco 
+ \ — 

(O 



ft) dco 

+ TT-^ \\^-''-fdsda (2). 

jje-'^'^Fd.dJ^ '^ 
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Adding together (1) and (2) we find 

"^ dco CO 

+ , , ^ W e-'^'x idh' '{^, + 30, '^^ dsdcr. 

r r Q^ 

But now / / Fda = — ^ is constant whatever be the values 

of h and « or h' and a>. It is therefore not affected by the 
variations of h' and &>. 

These variations affect neither the quantities ai...7A, nor 
the limits of the da integration. 

We have then 

^^'=\aiogr+|9(A'x) + §9«, 

a result of precisely the same form as in the former case. 
135. In any case in which 

the proof of the second law for the two independent variables T 
and ft) can be presented in a simple form. 

For we have in that case from II 

BE ^ _ . „^ ^^ do) 

-y =xaiogr+fx-. 

Now let i=z -^^ that is i is the time in which a particle 

moving with the velocity of mean square would describe a 
distance proportional to the linear dimensions of the system. 

Then §3 log ft) = 29 log i + d log T, 

and therefore 

f = xaiogT+|x^^ 

= 2Xd log (iT). 
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We might define \2iT as the Action of the system par excel- 
lence. With that definition the Action is the same for all 
states of the system defined by dE = 0, that is for all states of 
the system on the same adiabatic line. And we may say that 

the thermodynamic function <^, ( ^ = 9^ ) either is, or is propor- 
tional to, the Action. 

136. But dx = -r^dco cannot exist for any system such as 

we have treated as a gas, except for a so-called perfect gas in 
which ;\; = 0. If there be a potential ^ ^^ ^1^> ^^^ relative 
frequency of different configurations of the molecules must 
generally vary when T varies. And therefore although for any 
configuration the change of % as o) varies is independent of T, 
yet the mean of all values of x will not be independent of 
Ty because the variation of T causes some configurations, i.e. 
some values, to be represented with greater, and others with 
less frequency. We may for instance conceive a case in which 
the molecules placed at equal distances from each other would 
exert no sensible forces on each other by reason of their 
distances, and so % = 0. If we set them in motion they will 
in some cases approach so near to each other that % =|= 0. 

137. Professor J. J. Thomson in the work above referred 
to gives a proof, not, as he says, of the second law of Thermo- 
dynamics, but of a proposition analogous thereto. He does not 
confine himself to the two variables T and co, but employs T 
with any number of controllable coordinates cf). Then he 
assumes that the mean potential energy ^y or in his notation 
F, is completely fixed by the controllable coordinates, that is, 
is not a function of T. That restriction prevents Professor 
Thomson's proof from being applicable to gases, as a proof that 

dE . 

^ is a function of temperature and volume, except in the 

limiting case of a perfect gas in which V is non-existent. Nor 
does Professor Thomson propose to apply it to such cases. 
It may well be that in the more general class of problems 
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treated of in his work the condition that V shall be a function 
of the controllable coordinates is not a serious restriction. 

138. The following is another form of solution of this 
problem. Let there be, instead of the one controllable coordi- 
nate ft), any number of them as ft>i ... con- We will now use co 
as a type controllable coordinate. It is assumed that the rate 
of change of any co with the time is so slow that we may leave 
the velocities Wi . . . cb^ out of account. Let Fi . . . V^ be the 
generalised components of momentum corresponding to ft)i . . . con. 
Then it is a characteristic of the stationary motion that every 

dV 
_^ is zero. And therefore ii pi ... p^ be the components of 

external force required to maintain the motion stationary 

p = ~^ T^ for each ft), 

aft) d(d 

by Lagrange's equations. See Art. 24. 

Generally let F be the function denoting frequency, that is 
jF is a function of the coordinates Xi.., Zn, and the velocities 
ai-.-7n, ot the molecules of the system such that the time 
during which on average of any very long time, they lie between 
the limits 



«! . . . ofi -h da^ , 



is Fdx^ ... dzndai ... d^^n, ^^ as we may write for brevity Fdsda. 
Then the mean value of any function as % is 

X=ll X^d^d^ / / 1 Fdsda-, 
^ = f f ^ FdsdcT j jj Fdsda, &c. 

And dx = [[ ix^F + F^X) ^^^^ / // ^^^^^- 
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In the same case we have, if r represent the kinetic energy 
in any particular state, T the mean value of r, 

dT= [j rdFdsda -h jj Fdrdsdo- 

= [jrdFdsda +jJF^ daydsdcT. 

139. Let us at this point assume the function ^ to be a 
product of two functions, F=ff\ where/ is a function of the 
coordinates only. In that case we have also 

X=\xfds jjfds, 

or dx-^day^jxdfdsljfds. 

And now 

^dT + dx-^dco + ^day 
^ d(o d(o 

= 2ar H- f x¥ds - [[ rdFdsda, 

because -7- dco= 11 F -r- dcodsda 

d(o JJ d(D 

= dT- [jrdFdsda; 
therefore 

^-^ = 2dlogT + 2Jpfds-jj'^dFdsda. 

dE 
And now we make -^ a complete differential by assummg 

where <^ and y^ are any functional symbols. 
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A solution of these equations is 

where (7, C are the known constants. 

This sohition is the one we should take under the assumption 
of condition A. Another solution is 

if Q — ^T, This is the solution corresponding to Chap. VIII. 

The subject of the Second Law of Thermodynamics has 
been very fully investigated by Professor G. H. Bryan in his 
Report on the present state of our knowledge of Thermo- 
dynamics, Part II., British Association, 1894 
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MATHEMATICAL PROPOSITIONS. 



In the preceding chapters I have assumed without proof 
certain propositions relating chiefly to determinants. They are 
for the most part very simple, but as some readers may think 
they require proof, I will devote this chapter to proving them. 
I shall require to make use of the following elementary 
property of determinants. If each of the constituents in one 
row is equal to the corresponding constituent in a second row 
multiplied by a constant factor, the determinant is zero. The 
two rows for which this is the case are said to be similar. 
Whatever general proposition is proved of rows is of course 
equally true of columns. 

(a) Let us take for a type determinant of n^ constituents 



D: 



^21 > ^22 • • • ^271) 



^ni> ^n 



Cn 



in which Gpq is the constituent (whatever its value) in the pth 
row and qth column, and C^^ is not necessarily equal to G^iy &c. 
This can be made to represent any given determinant by 
assigning proper values to the constituents (7. 

Let Dpq denote the minor of D formed by omitting the pth 
row and qth. column, that is by striking out Cpq and the row 
and column in which it is ; Dqp the minor formed by omitting 
Cgp. In the determinant or any of its minors the diagonal 
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from the first constituent in the top row to the last in the 
bottom row shall be called the afdsy and its constituents axial 
constituents. 

(b) The sign of Dpq shall be such that 

D = GiiDu -\- C12D12 + . . . + CmDm 
= 601-D21 + C22D22 + . . . + ConJ-^^n 

= &c. 

In order that this may be the case, we must in forming 
any minor Dpq apply the following rule of signs, namely to the 
product of all the axial constituents of that minor the positive 
or negative sign is to be attributed according as j9 + ^ is even 
or odd, and we deduce therefrom the sign of the product of any 
other set of constituents of the minor according to the usual 
rule in expanding determinants. In like manner Dprqs shall 
denote the minor formed by omitting Cpq and (7^,^ and the rows 
and columns in which they are. It follows from the above 
convention concerning signs that 

p=^nq = n 



pq -^pq • 



Then evidently ^^^ = D^,,^, whatever values p and q may 
have. 



dD 



^pq 



(c) Now consider the new determinant 

O21 — "2^1} ^22 t/2^2 • • • ^2n "^P^m 
^ni — ^n*^i y ^■^n2 f^'n/^2 • • • ^nn ^n^n • 

If all the (7s in any two rows, as the pth and ^th, were zero, 
D' would be zero, because the pth and qth rows would then be 
similar. It follows that every term in D contains at least 
n-1 C factors, and therefore that D' contains no products of 
any two or more of the new partial constituents 6a;. For such 
products necessarily enter in pairs equal and of opposite sign, 
as for instance d^x^O^x^ — O^x^Jd^^. 



Therefore 



U = D-tX0pxM 



pvOq-i^pq . 
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(d) Now form a new determinant of {n-{- ly constituents 
by adding to D an (n H- l)th row and {n + l)th column. The 
(n -h l)th row shall be 

P (71+1) 1 > P (n+i) 2 • • • P in+i) n y ^Ctn+i > 

and the (n + l)th column shall be symmetrical with it, contain- 
ing y8i(«+i) ii^ tbe first row, ^2{n+i) i^ the second, and so on. 
Also let ^pq = ^qp, for all values of p and q. 

Let A be that new determinant, namely 

^11 > ^12 ••• ^m> Pi (n+i) J 
^21 y ^22 • • • ^2n » P2 (n+1) > 



A = 



(e) Next in D' write 

/) ^ __ A(n+i) ^ /? /« _ Pi (n+i) ^2 (n+i) 

/J ^ __ P2 (71+1) /3 ^ _ Pp (n+l) Pg (71+1) o^ 



2Cl^7l+l 

Then we have 



2a^ 



Pj? (71+1) Pg (n+1) 



J^po) 



and by the ordinary expansion 



A = 2a^+i D - SS^^ (n+i) y8< 



Therefore 



i)': 



2a, 



*^ (7i+l) Pq (71+1) -^j05 



A. 



'71+1 



(/) In this last deduction we have attributed no particular 
values to. the C constituents. 

It would therefore hold for all the minors of D' not contain- 
ing the suffix n + 1 as well as for D' itself, namely ' 



Dn' = 



D' = 



1 



■A„, 



2a, 



A12; 



71+1 



and therefore 



£>' A ' 



10 
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{g) Now consider the function e'^Q in which h is any con- 
stant, and Q is a quadratic function of ii variables u,, u,,., Un, 
namely 

Q = a{ti^^ + 6i2«^i^2 + «2^2^ -f &c. 
Let us effect the integration 

j e-^^du^ between limits ± oo . 

The terms in Q which contain Un are 

which may be put in the form 
Let us for a moment write 

im% -I- h^U^ 4- &c. 



= ^. 



Then 



and 



2Va« 

J 00 

= e-^Q' multiplied by a numerical factor, 



h 2 



^2*' + &C. 



Comparing these results with (e\ (/) we see that if 
^12, 2^2, h<^,..h^^ 



D containing (>i+ 1)^ constituents, 



and U : 



2a ^ ^lAn 



^1^ — 



2a^ 
^12 — ?;^ — , la^ 



2an 



2an 



'12 ^ , ^^2 — ^r— , . . . , 



D' containing ti^ constituents ; 
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,._ A/ Ai A/ A. 

then W = ^' -W^T- 

That is, these ratios are not altered by the change in the 
coefficients aj, 612 » &c., consequent on the integration accord- 
ing to Un. Evidently the same law holds for every subsequent 
integration. Let us then perform the integrations 



///■ 



€~^^dUs ... dUn, 

JJJ 

and let the result be 



Ill' 



If e-^ (Aiu,''+B,^u,u^+A^u^-i) duidu^ is proportional to the chance 
that two variables Ui and u^ shall lie between the limits 
% . . . Ui-Y dui for the one, and U2...U2 + du^ for the other, we 
easily find the mean of u-^, thus 

the limits of integration being for each u±^ 

1 2J.2 dn 



h 4>A,A^ - A/ ^ci ' 


if d be the determinant 




22I,, A2, 




A2, 2^,: 


But ^ = "^' by the last article. 


Therefore '^'^hD' 


A2 
Similarly UiV^ = j^, 




and so on. 





10—2 
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(h) Under the same circumstances to find the mean value 

dO 
of -~ , or 2aiU-^ + 6l2^^2 + ^13^3 + &c., iij being given (see Art. 97). 

That is to effect the integration 

. ,.dii2 . . . dun€~^^ (2a-iUj + b^^^io + &c.). 



///• 



Let 



2aiUj + 612^2 + ^i3^^3 + &c. = P. 



The coefficient of u^ in P is twice the coefficient of Ui^ in Q. 
The coefficient of every other u, as tio, in P is equal to the 
coefficient of U1U2 in Q. 

After integrating according to u^ the coefficient of Ui in P 
J 2 J 2 

is 2ai — ^ while that of u^^ in Q is a^— —^ . 



^in^2n 



2a^ 



and the coefficient 



The coefficient of U2 in P is 612 ■ 
of ^1^2 ill Q is the same. 

The same relation between the coefficients in P and those 
in Q which existed before integration according to Un exists 
also after that integration. The same rule holds for every 
subsequent integration to U2 inclusive. But when all the 
integrations have been effected to Wa inclusive, Q has become 

r-fr-Ui^, Therefore P has become ^ u^ which is the mean 

2Dii Un 

value sought. 

(i) Let us now consider the determinant 

D= S\ -f,S, -M -f,S,... 

-AS, AA, A^A^> (1+/A-- 

treated of in Chapter v, Art. 63. Here S is written for 2/ 
of that Article. 
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Comparing it with the type determinant, we find 

G^^^8^ o,, = l+/,^ (?33=l+/2^&c., 

C, = C, = ^M &c., G^ = (732 =/i/2, &c. 

Let it be required to find the values of D and its first 
minors. 

(j) We observe that in D every row except the top row eon- 
tains one constituent of the form 1 +f^, e.g. the {p+ l)th row 
contains 1 +//. If it were fp^ instead of 1 +//, the (p-t l)th 
row would be similar to the top row and therefore D would be 
zero. 

If we strike out any column other than the first, say the qth. 
column, we deprive the qth row of its constituent 1 -\-fq-.^. 
Hence the gth row is now similar to the top row. 

And if, to form a minor, we strike out any other row, neither 
the first nor the qth, the determinant, i.e, the minor so formed, 
is zero. We see then that ii p>l, and g > 1, and p^q, Dpq = 0. 
But if we expand D, the constituent Gpq appears only in the 
form DpqGpqy and Dpq = 0, Therefore I) when expanded does 
not contain any anaxial constituent except those in the first 
row or first column. 

From (i) and (j) we see that D = S\ For the same reason 
D^=S^, and so on for every coaxial minor except Du, but as we 
shall see later Ai = S/' + 1. 

(k) We have next to find A2, As, &c. Striking out from 
D the top row and second column to form A2, we find that in 
A2 every row except the top row has one constituent of the 
form 1 +f\ Therefore by the same reasoning as in (j) and by 
the rule as to the sign of the product of the axial constituents, 
A2 = Ai = +/i'S. When we form As? the product of the axial 
constituents is by our rule to be positive. But in order to 
reduce this minor to one in which the top row does not 
contain a constituent of the form 1 -h/^ we must interchange 
the two top rows, and thereby change the sign, of the minor. 
Therefore As = -f./s'S' and so on. In forming any minor D^p we 
have if p is odd, an odd number of changes of sign. 
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(l) Next let us deal with Ai- 

Let A=Ai= l+/x^ AA, /;/;,... 
/1/3, /^/s, l+/3^... 



[appendix 



or 



(7ii=i+/?, tv,==/i/;,&c. 



All products of the form /^^/^^^ evidently disappear in 
expansion, and therefore all terms involving f^ except in the 
first degree disappear, and therefore 

n-l 

A: = 1 + 2 /^ 

1 

= |/^if/„=l. 



(m) Again, if we strike out any column except the first, say 
the second, we thereby deprive the second row of its constituent 
1 4-/2^. And therefore if we now strike out the top row we 
shall have 

/l/2> /2/3J 72/4 > ••• 

A A, (1+//), AA,... 
yi/4, /a/., (1+//). •• 



= —f^ f^ by the rule as to the sign of the axial constituents. 
Similarly A13 = -/1/3. 

Also A,,==S/^-/,^ 

A,,=2/^-/,^&c. 

{n) In the type determinant let now every G become 
G + dG, and let all the 9C's be small, so that we may neglect 
products of them. That being the case 

2? = i q=zi il^pq 
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Let US apply this to the determinant of (i) in which 
On = s% 0., = -M C,, = -/A &C., 

Here as we have seen -j^ = 0, unless either p=i, or q = l, 
or p = q. 
Therefore 

dB = DndGn +ASdG,, +/^(7„ + &c. 
+fiSdC,, +/AdG,, + &c. 
+ S''dG^ + S^G^ + &,c. 
Now as in Art. 83 let us make 

dGn = 0, aOi2 = ^b^„S, dG,, = ^b^S, &c. 

dGa = ibinS, dGsi = ^b^S, &c., 

9(722 = — bmfi , BCss = — 62„/2, &c. 

Then as products of the 6's are to be neglected, we have 

dD = + hn/A' + b^/A' + hnfsS' - &C. 
-b^nfiS'-b^/A-SiC. 

= 0. 

(o) Let us treat D^, or A, in the same way. Here 

Cu = l+/.'. G,,=f,f„&o. 
And aAi = dA, 

dA _ _, dA - „ „ 

And 8(7n = - hnfi , 9^22 = -h^f<,, &c., 

Therefore, since each anaxial term/1/2, &c. appears twice, 

a A = - (2/ -/,^) hr^f, - (S/^ -/2O 62^/2 - &c. 

-/1/2 (*12 - ^in/2 -- &2n/i) -/1/3 (&13 - bmfs - bsnfi) - &C. 
= "~ ^12/1/2 — Oiajijs — ... — Opqjpjq 

+f^%nfx +/A„/i + . . . +/A„/2 +/A„/2 + &C. 
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That is 8A = - 612/1/2 - 613/1/3 " &c. 

In this result every bpgfpfq is to be counted once only, that 
is, we do not distinguish hpq from bqp. If bpqfpfq and bqpfqfp 
are counted separately, we must write 

and this form is the one required in Art. 83. 

(p) To prove that for certain values of b, Xfp^bpqfq = ^fp^^b 
see (33) of Art. 83. Let /(Art. 61) = 4r2 when r>a. 

Taking for origin, let p be any point Op = p. About p 
describe a sphere of radius r less than p, and let g be a point 

on its surface. If ti + 2 = 1, or /= - , the mean value of fq for 

all points q on the sphere is ^ , ovfp. And equation (33) is 
accurately satisfied, so long as r<p. 

If n 4- 2 > 1, we have in usual notation 
sin ed0 



Jo 



n+2 

(p^ + r^ — 2pr cos 6) 2 



1 (p + ry -{p- ry 
2npr {p^ — r^Y 



n. n — ln— 2 ^ , ^ . 

jipn-i^ + ^r-^ pn-3^ ^ ^c. 

1 Z . o 

— ^p^ (^p2 _ ^.2yi 

=/, -h terms containing /p multiplied by ~2*^c. 

p 

And if now 6 = - ( - j , where /-t > 3, 

fpf 4.7rr'pbfqdr = fj^'j ^.irr'pbdr =fp'lb, 
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because, - being negligible, all the terms in the integral which 
follow the term fp^Sb are negligible compared with fp^%b, 
(q) To find the condition that the quadratic function (q) 
Q = aiiii^ + bi2U^U2 + a^a^ 4- . . . 
shall be always positive whatever values be attributed to 

III i ^2 • • • ^n • 

Evidently since ^2 . . . Un may be all zero, a^ must be positive. 
Similarly ag, ag, &c., must be positive. 

Also, given Un, Q is minimum when 
dQ 



dUn 



= 2anUn + hnUi + h^^U^ + &C. = 0. 



And, when minimum, Q must be positive. Substitute in Q 
for Un the value found for it from 



that is 



"^dnUn + binVn. + b^U^ + &C. = 0, 



tan 
Q then becomes 

Q' = aiUj^ + bi2thU2 + a2U2^ + &c., 

J 2 J 2 

where ai^a^-—^, cia' = ag - ^ , i&c, 

and every a' must be positive. Also 612' = 612 — ^- &c. 

Proceeding in the same way Q is ultimately reduced to 

^^oi 1- . , , D 

V = -^TfT '^*i > and as this must be positive, yp must be positive. 

Similarly ^ must be positive, and so on. Therefore D and all 
-^22 

its first coaxial minors have the same sign. But if we only 

reduce Q to a quadratic function of Ui and U2 
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From which and the corresponding equations we see that 
all the second coaxial minors have the same sign as the first 
coaxial minors, and therefore the same sign as D. 

Similarly we can prove that the third coaxial minors must 
have the same sign as the second and so on, so that D and all 
its coaxial minors must have the same sign. And this sign 
must be positive because aj, ag, &c., are coaxial minors. This 
is the required condition. 

(r) Referring to Chapter v, in which it is proved that with 
(7^) spheres of finite diameter c, 

dt ~ "^^^ 2h\dz'^ dx) ' 
where M^^ or itf = 1 1 jdwdydzay ( j + j j , 

to prove the same thing for 



Ma-,:=jjjdxdydza^^. 



'^XX ' 



Consider the space between two infinite planes ^ = and 

d^ 
X = dx, between which ^ is supposed constant and negative. 

Let u—OL — ^. Then it is proved, as in Art. 64, that the number 
per unit of volume of molecules for which a — | lies between u 
and u -f du, when x = 0, exceeds the corresponding number when 
x^dxhy the quantity 

Ge-^''''2h'u'^ dxdii. 
ax 

Therefore the number which enter the space between the planes 
exceeds the number which leave it per unit of time by 

G Udydze-^'^^du^h'ua ^ . 

Hence as in Art. 66, it follows that, for each element of 
volume dxdydz, 
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ax 

Now u and a are not independent, as u and 7 were in 
Art. 66. Hence 



2/1' " 2/1' 4/^2 2/t'2/i ' 

and 2/?,Va2 = — . 

2/1 

^^''''^ dxdt'dxlti' 2h\dx) ^ ^' 

which corresponds to 

\dz doc) dt 2h \dz dx) ' 

of Art. 66. 

It follows, the motion being stationary, that 

dd?^3^/dfY 
dtdx 2h \dx) ' 

or by the process of Art. 68 

When the material points are replaced by spheres of finite 

da.^ 
diameter c, -77 acquires a new term due to collisions, which we 

will denote by ;=— , and we have 

d^da^^d^dVl 
dx dt dx dt ' 

^x^ Vyy Vzf being the components of the half relative 
velocity of two colliding spheres. 
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Also ^^' = 2^cV F(\'^ - XO F^ 

where \V=V^, fiV=Vy, vV = V„ 

and using the same notation as in Arts. 70 — 71, we find 

77 

^p = - 'Ittc'p V4h ^ I V'X f'de r d(l>2 sin cos 

dt IT '^ dxz Jo Jo 

X (- X cos ^ + Vl - V sin cos </>') jX^ (4 cos^ ^ - 4 cos^ ^) 
+ ( I - A.^) sin^' 2(9 cos=^ ^' - 2X ^l^-X" cos 2(9 sin 20 cos <^'j 

4 

whence introducing the factor Ge-^^'V^dV where t'=-j- (2A.)5 

and integrating from F = to F = oo , 

aF,^ 8 , Id^ 
we get -ar^Is'^'^Md^' 

d« dt'dxdt Vo^^<^th\dx) ^^^• 

We have next to calculate the term 

dtdx' 

due to collisions. That gives by the same method as that 
employed in Art. 76 

^ ^= _ _1 27rc^pF^ \^d0 f #2 sin 6 cos ^ 

ot ax TT j i 

X (— X COS ^ + Vl - X^ sin cos </>') 
X (- 2X cos^ (9 4- Vr^2 sin 2^ cos (^') j-{^= ^)' 
Whence we obtain for a pair of colliding spheres 

dtdx ^ 4i5 rdr^ 
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and for each sphere separately 

,9 <^? . , T7-. 22^,1 df , , 



dt dx 45 r dr 



-S"''>a(S)'''y<^) w- 



Adding this to (3) we find for the change of M^x with the 
time due to collisions 



dt 15 '^^^'>^^^ 



2h \dx) 
Using Mxx for M^^ or ilf in Art. 79, we find 

-"vA{(i)'^(2n-i(f)'-^ 

The first line is zero. The second can be made zero by 
making [-j-] • ( 3^) :: a: : 1 — /c as before. 
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